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SECTION  I 


INTRODUCTION 


The  general  objective  of  this  study  is  to  develop  techniques  and  tools 
necessary  for  rapid  design  of  an  active  control  system  for  aircraft  with 
lightly  damped  structural  modes.  The  synthesis  techniques  provided  here 
are  aimed  at  reducing  the  engineering  man-hours  presently  required  for  flight 
control  system  design,  thus  effecting  a  cost  reduction.  Improvements  in  the 
fatigue  life,  ride  qualities  and/or  handling  qualities  of  military  aircraft 
are  sought  by  controlling  the  lightly  damped  modes  and  thus  improving 
their  mission  performance. 

The  present  scope  of  this  study  is  to  develop  programs  to  interface 
the  level  2.01.00  FLEXSTAB  computer  program  system  with  existing 
Air  Force-owned  optimal  control  computer  programs.  These  programs 
represent  advanced  computational  techniques  required  to  perform  quantita¬ 
tive  analysis  of  multisurface  control  systems.  The  resulting  interface 
program  system  is  called  "KONPACT  -  Computer  Programs  for  Active 
Control  Technology.  "  Working  together  with  FLEXSTAB,  KONPACT 
provides  the  capability  to  model,  synthesize,  analyze,  and  design  automatic 
control  systems.  It  can  also  be  used  as  a  stand-alone  program. 

The  work  performed  under  this  contract  is  reported  in  three  volumes: 

Volume  I.  KONPACT  Theoretical  Description  and  Demonstration 

Volume  II.  KONPACT  Program  Listing 

Volume  III.  KONPACT  User's  Manual 
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This  document  reports  the  user's  information  needed  to  execute  KONPACT. 

It  also  contains  two  C-5A  demonstration  examples  and  results. 

The  purposes  of  including  these  demonstration  examples  are  first  to  document 
the  input  data  and  resulting  output  that  are  required  to  design  the  C-5A 
handling  quality  controller  and  the  C-5A  Active  Load  Distribution  Control 

System  using  the  FLEXSTAB  C-SA  model.  This  model  is  documented  in 
References  8,  9  end  11;  the  designed  controllers  are  discussed  in  Reference  4. 

The  second  purpose  is  to  provide  concrete  examples  oi  the  inputs  that  are 
required  to  execute  KONPACT.  A  majority  of  the  input  options  available 
to  the  user  to  execute  KOMPACT  are  contained  in  these  examples.  In 
addition,  the  resulting  outputs  from  the  KONPACT  input  deck  setups  are 
displayed  to  familiarize  the  user  with  the  output  formats.  It  must  be 
noted  that  in  this  report  the  terms  "  vehicle  model"  and  "  vehicle"  are  used 
interchangeably;  both  terms  refer  to  the  mathematical  models  generated 
from  the  level  2.01.00  FLEXSTAB  linear  systems  analysis  program 
(Reference  1), 

Section  II  presents  briefly  the  description  of  KONPACT  programs.  A  short 
description  of  the  design  process  is  presented  first  for  completeness. 
Subsequently,  a  general  description  of  input  cards  to  execute  KONPACT  is 
given. 

In  Section  III  the  program  execution  procedure  for  the  modeling  program 
(KONPACT-1)  is  described.  The  formats  of  the  input  data  deck  and  printed 
outputs  are  also  described  in  detail. 
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In  Section  IV  the  program  execution  procedure  for  the  design  program 
(KONPACT-2)  is  described.  The  formats  of  the  input  data  deck  and  printed 
outputs  are  also  described. 

Section  V  contains  the  demonstration  example-  for  the  handling  quality 
design  for  the  C-5A  FLEXSTAB  vehicle  model.  The  deck  set-up  and  the 
printed  results  are  displayed,  and  the  results  are  discussed. 

Section  VI  contains  the  demonstration  example  for  the  ALDCS  controller 
design  for  the  C-5A  FLEXSTAB  vehicle  model.  The  deck  set-up  is  displayed, 
and  the  final  results  are  discussed. 

The  analytical  techniques  and  algorithms  used  in  KONPACT  are  described 
in  Volume  I.  Volume  I  also  demonstrates  how  these  techniques  are 
applied  to  flexible  aircraft  control  system  designs. 

Documentation  of  KONPACT  was  beyond  the  scope  of  this  contract.  To  aid 
the  user,  the  listings  of  source  programs  which  implement  the  mathematical 
analysis  and  models  presented  in  Volume  I  are  given  in  Volume  II, 


SECTION  II 


DESCRIPTION  OF  KONPACT  PROGRAM 


INTRODUCTION 

KONPACT  is  a  system  of  computer  programs  developed  by  Honeywell 
under  Air  Force  Contract  No.  F33615-75-C-3046.  KONPACT  uses  the 
state  space  approach  for  modeling  flight  control  systems,  and  it  designs 
the  controllers  using  optimal  control  methodology.  KONPACT  interfaces 
with  the  linear  systems  analysis  (LSA)  program  of  the  Level  2  FLEXSTAB 
Program  system  developed  by  Boeing  under  Air  Force  Contract  No. 
F33615-72-C-1172  (Reference  1).  KONPACT  can  also  be  used  as  a  stand¬ 
alone  program. 

KONPACT  operates  on  CDC  6000  and  CDC  7000  series  computers  and 
can  be  easily  modified  to  operate  on  other  computers.  KONPACT  has 
been  written  in  Extended  Fortran  IV  language. 

In  this  section,  the  design  process  is  reviewed  and  motivation  is  given 
for  developing  KONPACT  program.  A  general  description  of  KONPACT 
program  is  presented  in  terms  of  overlay  organization  and  information 
flow.  The  input  cards  for  KONPACT  are  also  described. 
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REVIEW  OF  DESIGN  PROCESS 


Basically,  optimal  control  design  of  flight  control  systems  involves 
two  major  steps.  The  first  step  is  to  obtain  the  state  space  description 
of  the  flight  control  system.  The  second  step  is  to  obtain  optimal 
feedback  gains  via  the  optimal  control  methods. 

The  flight  control  system  consists  of  the  basic  vehicle,  sensor  equations, 
load  equations,  actuator  dynamics,  controller  specifications,  Wagner- 
Kussner  dynamics,  etc.  Level  2  FLEXSTAB  calculates  the  basic  vehicle 
model,  sensor  equations,  and  load  equations.  The  remaining  dynamics 
of  the  flight  control  system  are  input  by  the  user,  and  the  KONPACT 
program  combines  and  calculates  the  flight  control  system  model. 

The  optimal  state  feedback  gains  are  computed  using  the  DIAK  program, 
and  the  feedback  gains  are  reduced  to  gains  only  on  specified  measure¬ 
ments  using  the  FFOC  program. 

OVERLAY  ORGANIZATION 

KONPACT  consists  of  two  programs;  a  modeling  program  (KONPACT-1) 
and  a  design  program  (KONPACT- 2).  KONPACT-1  interfaces  with 
FLEXSTAB  through  the  LSA  program  to  obtain  the  vehicle  model  and 
augments  the  specified  dynamics  to  obtain  the  state  space  description 
(quadruple  data)  of  the  flight  control  system.  These  data  are  utilized  by 
KONPACT-2  (which  contains  the  subprograms  DIAK  and  FFOC  that  are 
described  in  Reference  2)  in  the  design  of  the  optimal  feedback  gains.  Also 
KONPACT-2  interfaces  with  FLEXSTAB  through  the  LSA  program  to 
evaluate  the  performance  of  the  above-designed  optimal  flight  control  system. 


Table  1  provides  a  brief  description  of  programs  KONPACT-1  and 
KONPACT-2  and  their  subprograms.  The  interface  between  KONPACT 
and  the  LSA  program  is  illustrated  in  Figure  1.*  The  overlay  structure  of  the 
KONPACT-1  program  is  illustrated  in  Figure  2.  It  consists  of  a  main 
overlay  and  five  primary  overlays  (Reference  3).  The  overlay  structure 
of  the  KONPACT-2  program  is  illustrated  in  Figure  3.  It  consists  of  a  main 
overlay  and  three  primary  overlays. 


Table  1.  KONPACT  Program  Descriptions 


PROGRAM 

SUBPROGRAM 

DESCRIPTION 

KONPACT-1 

State  space  modeling  program 

STAMK1 

Obtains  state  space  model  from  LSA 
simulator  deck  data 

STAMK2 

Obtains  state  space  model  from  transfer 
function  data 

STAMK3 

Obtains  state  space  model  from  quadruple 
data  and  interconnection  data 

STAMK4 

Obtains  state  space  model  from  simulation 
equations  (user  written) 

CONDK 

Modifies  the  state  space  model  by  scaling, 
shuffling,  truncating, and  residualizing 
the  system  variables 

KONPACT-2 

Optimal  design  program 

DATAK 

Prepares  data  for  DIAK,  FFOC  and  LSA 
programs 

DIAK 

Designs  full  state  feedback  optimal  controllers 

FFOC 

Designs  reduced  state  (practical)  feedback 
optimal  controllers 

*  Illustrations  start  at  page  No.  75 
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INFORMATION  FLOW 


The  normal  sequence  for  obtaining  an  overall  state  space  model  of  a  flight 
control  system  using  the  modeling  program  (KONPACT-1)  is  as  follows: 

•  The  vehicle  model  is  obtained  by  using  either  the  STAMK1  (for 

LSA  data)  or  STAMK4  (for  other  types  of  vehicle  data)  subprograms. 

•  The  actuator,  sensor,  controller,  implicit  and  explicit  models 
are  obtained  by  using  either  the  STAMK2  (with  transfer  function 
input  data)  or  STAMK3  (with  quadruple  input  data)  subprograms. 

•  The  subsystems  defined  above  are  combined  to  get  an  overall 
system  by  using  the  STAMK3  (interconnection  data)  subprogram. 

•  The  overall  system  model  is  conditioned  (modified)  by  scaling 
and/or  shuffling  and/or  truncating  and/or  residualizing  the 
variables  using  the  CONDK  subprogram.  This  subprogram  also 
develops  the  rate  of  change  of  response  variables  when  required. 

The  normal  sequence  for  designing  optimal  feedback  controllers  and 
evaluating  the  performance  of  the  resulting  system  using  the  design 
program  KONPACT-2  is  as  follows: 

•  Full  state  feedback  control  gains  are  obtained  by  varying  the 
quadratic  weights  and  using  the  DIAK  subprogram. 

•  The  resulting  full  state  feedback  control  gains  are  reduced  to 
gains  only  on  specified  measurements  by  using  the  FFOC 
subprogram. 
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•  The  performance  of  the  resulting  closed  loop  system  is  evaluated 
using  the  LSA  program. 

•  The  above  steps  are  repeated  until  a  satisfactory  design  is 
obtained. 


Table  2  describes  the  data  tapes  used  in  KONPACT-1  and  KONPACT-2 
programs.  The  state  space  model  data  (quadruple  data)  and  the  name 
list  data  are  written  on  tapes  QDATA  and  NDATA,  respectively.  The 
vehicle  data  (simulator  deck  data)  is  written  on  tape  VDATA.  The  feedback 
gain  data  from  DIAK  and  FFOC  are  written  on  tapes  DDATA  and  FDATA, 
respectively.  The  overall  system  data  in  frequency  representation  form 
is  written  on  tape  SDSTP  for  use  by  the  LSA  program.  The  DATAK 
subprogram  is  used  in  preparing  data  tapes  for  DIAK,  FFOC,  and  LSA. 


Table  2.  KONPACT  Data  Tapes 


Tape 

Name 

Contents 

Generating 

Program 

Benefiting 

Program(s) 

VDATA 

Simulator  Interface  data  in  the 
form  of  card  images 

LSA 

KONPACT-1 

QDATA 

Quadruple  (A,B,  C,D)  or  state 
variable  representation  data 

KONPACT-1 

KONPACT-1 

KONPACT-2 

NDATA 

Name  list  data  of  the  state 
variable  representation 

KONPACT-1 

KONPACT-1 

j  DDATA 

Full  state  feedback  gain  data 
in  the  form  of  card  images 

KONPACT-2 

KONPACT- 2 

FDATA 

Reduced  feedback  gain  data 
in  the  form  of  card  images 

KONPACT-2 

KONPACT-2 

SDSTP 

Frequency  domain  representa¬ 
tion  of  quadruple  data 

KONPACT-2 

LSA 
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VARIABLE  DIMENSIONING 


Variable  dimensioning  (Dynamic  Data  Storage)  technique  (Reference  7) 
is  used  for  efficient  data  storage.  This  technique  also  facilitates  changing 
the  amount  of  allocated  (required)  storage  space  by  a  data  card  input.  In 
KONPACT-1,  the  subprognm  arrays  (whose  sizes  depend  on  the  maximum 
system  dimension  inputs)  are  stored  in  scratch  storage  blocks  using 
variable  entry  points.  In  the  subprograms  the  arrays  are  dimensioned 
with  integer  varirbles.  These  "variable  DIMENSION  statements"  remain 
unchanged  although  the  amount  of  required  data  storage  is  altered.  The 
maximum  size  of  the  scratch  storage  blocks  is  specified  in  a  "fixed 
DIMENSION  statement"  in  the  main  program. 

The  size  of  storage  actually  needed  by  the  arrays  varies,  depending  on 
the  maximum  system  dimension  inputs.  Thus,  if  the  maximum  size  a 
user  allows  for  his  problem  changes,  only  the  "fixed  DIMENSION  state¬ 
ments"  in  the  main  program  need  to  be  changed.  The  change  of  the  main 
program  of  KONPACT-1  is  done  by  a  precompiler,  as  discussed  in 
Section  III.  The  user  provides  the  new  maximum  system  dimensions  by  data 
cards.  Updating  and  running  with  updated  main  program  are  done  with 
control  cards  in  a  single  run.  For  more  details  on  variable  dimensioning, 
the  user  is  referred  to  Volume  II  (Reference  5). 

Modularization  and  variable  dimensioning  of  DIAK  and  FFOC  subprograms 
in  KONPACT-2  were  beyond  the  scope  of  this  contract. 
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INPUT  CARDS 
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The  first  logical  record  of  the  input  file  (deck)  contains  the  job  control 
cards,  and  the  subsequent  logical  records  are  the  data  deck  cards.  Figure  4 
illustrates  a  typical  input  deck  structure.  It  consists  of  job  control  cards 
and  one  logical  record  of  data  deck  cards. 

Job  Control  Cards 


Job  control  cards  are  statements  instructing  the  computer  how  to  process 
a  job.  Job  control  cards  are  always  the  first  group  of  cards  in  a  deck 
set-up. 

Data  Deck  Cards 


The  data  deck  input  to  the  KONPACT  programs  consists  mainly  of  conventional 
data  cards  containing  information  arranged  according  to  standard  Fortran 
formats.  In  addition,  the  data  decks  may  contain  three  special  card  types: 
Program  Control  Cards,  Data  Control  Cards,  and  Data  Comment  Cards. 

(For  detailed  format  descriptions  of  data  deck  cards,  see  Sections  III  and  IV. ) 

Conventional  Data  Cards --These  cards  contain  the  data  input  by  the  user  to: 

1)  describe  the  model,  2)  describe  the  interconnections,  3)  select  variables, 
and  4)  describe  the  variables. 

Numerical  data  are  generally  read  under  the  5(212,  El 2. 6)  Fortran  format, 
and  the  alphanumeric  data  under  the  20A4  format.  Some  data  cards  contain 
mixed  alphanumeric  and  numeric  data.  The  actual  format  used  by  each 
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data  card  or  set  of  data  cards  is  specified  in  the  text  and  on  figures  illustrating 
data  card  arrangements.  All  alphanumeric  data  must  be  left- justified  and 
all  numeric  data  must  be  right -justified. 


Program  Control  Cards --These  cards  contain  a  dollar  sign  ($)  in  Column  1 
followed  by  descriptive  words  which  identify  the  program  that  will  be  used 
to  process  the  data.  There  are  nine  program  control  cards  and  they  are 
listed  in  Table  3.  Only  Columns  1  through  4  are  interpreted. 


Table  3.  List  of  Program  Control  Cards 


Program  Control  Card 

$LSA  DATA 
$TRANSFER  FN  DATA 
$QUADRUPLE  DATA 
$INTERCONNECTION  DATA 
$SIMULATION  DA  f  A 
$CONDITIONING  DATA 

$DIAK  DATA 
$FFOC  DATA 


Subprogram 

Program 

STAMK1 

KONPACT-1 

STAMK2 

STAMK3 

STAMK3 

STAMK4 

CONDK 

DIAK 

KONPACT-2 

FFOC 

LSA 

(only  prepares 
data  for  LSA) 
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$LSA  DATA 
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Data  Comment  Cards--These  cards  contain  any  appropriate  comments 
the  user  wishes  to  make.  Column  1  should  contain  a  C  and  Column  2 
should  be  left  blank.  The  comment  starts  from  Column  3  and  can  go  up 
to  Column  80.  Any  number  of  comment  cards,  inserted  anywhere  in  the 
data  deck,  can  be  utilized  to  aid  the  KONPACT  user  in  identifying  important 
data  items. 

Data  and  Flow  Control  Cards --These  are  descriptive  words,  starting  from 
Column  1,  which  identify  key  blocks  of  data.  The  data  and  flow  control 
cards  are  listed  in  Table  4. 

Imbedded  Blanks 

The  user  input  cards  to  KONPACT  (Table  4)  contain  descriptive  words 

with  imbedded  blanks,  e.g.,  PRINT  INPUT.  DATA.  These  imbedded 

A  A 

blanks  are  necessary  for  the  execution  of  KONPACT  program. 
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Table  4.  Data  and  Flow  Control  Cards 
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Table  4.  Data  and  Flow  Control  Cards  (Continued) 


Table  4.  Data  and  Flow  Control  Cards  (Continued) 
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Table  4.  Data/ and  Flow  Control  Cards  (Continued) 


SECTION  III 


STATE  MODELING  PROGRAM  (KONPACT-1) 


INTRODUCTION 


KONPACT-1  program  interfaces  with  Level  2  FLEXSTAB  through  the 
linear  system  analysis  (LSA)  program  to  obtain  the  vehicle  model  in 
state  space  form  (quadruple  data).  KONPACT-1  also  augments  or  combines 
the  vehicle  model  with  the  dynamics  of  actuator  model,  controller  model, 
etc.,  as  specified  by  the  user,  to  obtain  the  state  space  model  of  the 
overall  flight  control  system.  The  user-specified  dynamics  can  be  either 
in  the  form  of  transfer  functions  or  simulation  equations.  KONPACT-1 
program  can  also  be  used  independent  of  FLEXSTAB  as  a  stand-alone 
program,  to  obtain  state  space  models  of  flight  control  systems  and 
other  control  systems. 

This  section  describes  the  execution  procedure  for  the  KONPACT-1 
program.  The  description  of  input  data  and  output  data  is  given  here. 

The  program  timing  and  central  memory  requirement  estimates  are 
also  discussed. 

PROGRAM  EXECUTION  PROCEDURE 

Figure  5  shows  the  block  diagram  of  data  flow  in  KONPACT-1.  The  MAIN 
program  reads  card  input  data,  reorganizes  it,  and  writes  it  on  BINPUT 
file.  The  BINPUT  file  is  used  as  card  input  data  by  the  subprograms. 
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The  job  control  cards,  used  for  executing  KONPACT-1  programs,  are  those 
commonly  used  by  the  SCOPE  3.4. 1  operating  system  for  CDC  6000  and  CDC 
7000  series  computers.  For  a  detailed  description  of  job  control  cards,  the 
user  is  referred  to  the  SCOPE  Manual  (Reference  3).  In  the  following  dis¬ 
cussion,  job  control  card  illustrations  are  based  on  using  permanent  disc 
files;  here  the  disc  files  are  used  to  store  output  data  required  for  succeeding 
runs.  The  user  is  referred  to  the  ASD  Computer  Center's  User's  Guide 
(Reference  6)  if  magnetic  tapes  are  used. 


When  executed,  the  FLEXSTAB/LSA  program  produces  the  simulator 
deck  data  for  the  vehicle  (Reference  1).  Figure  6  illustrates  the  job 
control  cards  for  storing  the  aforementioned  vehicle  data  on  VDATA  file. 
More  than  one  set  of  vehicle  data  can  be  stored  on  VDATA  file.  Each 

set  is  identified  by  the  first  card  in  the  set  which  contains  the  description 
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(i.e. ,  label)  of  the  particular  vehicle  model. 

Since  it  contains  card  images,  the  data  on  file  VDATA  are  formatted. 
KONPACT-1  program  does  not  accept  unformatted  data  on  VDATA  file. 

Job  control  cards  to  create  QDATA  and  NDATA  files  are  illustrated  in 
Figure  7.  There  are  two  logical  records  of  data  following  the  job  control 
cards.  The  first  logical  record  contains  the  specification  of  KONPACT-1 
subprograms  used  and  the  maximum  system  dimensions  needed  for  the 
job.  This  record  is  data  for  the  precompiler  program  PRECOM.  The 

*\ 

second  logical  record  is  the  modeling  data.  The  state  space  model 
(quadruple  data)  for  each  system  specified  by  the  user  is  written  on  QDATA 
file  and  the  name  list  data  for  each  system  is  written  on  NDATA  file. 

The  data  on  QDATA  and  NDATA  files  are  unformatted.  Figure  8  illustrates 
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the  job  control  cards  for  adding  new  system  data  on  QDATA  and  NDATA 
files  using  KONPACT-1  program. 

INPUT  DATA 

The  input  data  for  the  KONPACT-1  program  consists  of  three  groups. 

The  first  group  is  the  simulator  deck  data.  These  data  are  on  file 
VDATA  in  the  form  of  card  images.  The  second  group  is  the  data  for 
the  precompiler  program  PRECOM  in  the  form  of  cards.  The  third 
group  of  data  is  the  additional  modeling  data  also  in  the  form  of  cards. 

Simulator  Deck  Data 

The  simulator  deck  data  are  produced  by  FLEXSTAB/LSA  in  the  form 
of  a  punched  card  deck.  This  is  subsequently  input  on  file  VDATA  by 
the  user  for  KONPACT-1  program.  The  FLEXSTAB/LSA  equations, 
which  represent  the  simulator  deck  data,  are  shown  in  Figure  9.  The 
simulator  deck  data  structure  with  the  formats  is  shown  in  Figure  10. 

The  first  card  describes  the  simulator  deck  data  and  serves  as  a  header 
card.  If  more  than  one  set  of  simulator  deck  data  is  written  on  file 
VDATA,  then  a  particular  simulator  deck  data  can  be  retreived  by  specifying 
the  corresponding  header  card.  Following  the  header  card  is  the  matrix 
name  card  which  contains  one  of  the  standard  matrix  names  (defined  in 
Figure  9),  along  with  the  row  and  column  size  of  the  matrix.  The  matrix 
name  card  is  followed  by  a  deck  of  cards  which  contain  the  elements  of 
the  matrix.  The  formats  for  the  matrix  data  are  given  in  Table  5.  The 
last  card  of  the  simulator  deck  data  is  a  matrix  name  card  with  the  name 
’-FINISHED*  and  size  0x0. 
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Table  5.  Formats  for  the  Matrices  in  the  Simulator  Deck  Data 


CARD  1  (A10,  215) 

Column 

Contents 

Explanation 

1-10 

NAME 

Name  of  the  matrix 

11-15 

NROWS 

i 

Number  of  rows  in  the  matrix 

16-20 

NCOLS 

Number  of  columns  in  the  matrix 

CARD  SET  2  (6E10.0) 

1-10 

Aij 

Elements  of  the  matrix 

11-20 

Ai.  j+1 

Stored  by  rows,  six  elements  per  card 

21-30 

Ai.  j+2 

31-40 

• 

• 

• 

41-50 

Ai+1. 1 

51-60 

Ai+1, 2 

Precompiler  Data 

The  user  defines  the  KONPACT-1  subprograms  and  the  values  of  the  maximum 
system  dimensions  needed  for  executing  his  job.  The  KONPACT-1  basic 
subprograms  are  STAMK1,  STAMK2.  STAMK3,  STAMK4,  and  CONDK. 

If  the  user  does  not  specify  any  subprograms  to  be  used,  the  precompiler 
assumes  that  all  subprograms  will  be  used  in  the  execution  of  the  job. 


20 


* 


The  names  of  the  maximum  system  dimensions  (Table  6)  are  NXM,  NRM, 
NUM,  NYM,  MSB,  AND  MTB  where: 

NXM  =  Maximum  number  of  states 

NRM  =  Maximum  number  of  outputs 

NUM  -  Maximum  number  of  inputs 

NYM  =  Maximum  number  of  summing  point  equations  =  (no.  of  all 
internal  inputs  +  no.  of  all  internal  outputs)  -  (Reference  4). 

MSB  =  Maximum  number  of  systems  to  be  combined  (interconnected) 
at  one  time  s:  9 

MTB  =  Maximum  number  of  transfer  function  blocks  in  a  system. 

Only  nonzero  system  dimensions  need  be  specified.  The  formats  for 
precompiler  data  are  shown  in  Table  6, 


Each  system  (i.e.,  vehicle,  actuator,  controller)  in  KONPACT-1  is 
assigned  a  system  number  by  the  user.  System  number  9  is  reserved 
for  the  Implicit  model.  An  interconnection  of  one  or  more  systems  is 
defined  as  a  new  system.  KONPACT-1  program  obtains  the  state  space 
model  (quadruple  data)  for  each  system  along  with  a  name  list  data  for 
the  system  variables.  The  system  variables  are  the  states,  inputs, 
and  outputs  of  the  system.  When  the  state  space  model  of  a  system  is 
conditioned  (i.e.,  scaled,  shuffled)  by  KONPACT-1,  its  system  number 
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Table  6.  Formats  for  Precompiler  Data 


;  CARD  SET  1  (A6) 

Column 

Contents 

Explanation 

1-6 

NAME 

Name  of  the  subprogram 

CARD  SET  2  (A4,  13) 

1-4 

NAME 

Name  of  the  Maximum  System  Dimen¬ 
sion 

5-7 

MSD 

Value  of  the  Maximum  System  Dimen¬ 
sion 

is  retained.  The  state  space  models  (quadruple  data)  for  all  systems  are 
written  on  file  QDATA  along  with  the  label  cards  to  locate  data.  The 
name  list  data  for  all  systems  are  written  on  file  NDATA  using  the  same 
labels.  There  may  be  more  than  one  system  with  the  same  system 
number.  A  reference  table  associating  the  most  recent  system  description 
with  a  system  number  is  written  on  file  NDATA.  System  numbers  along 
witli  the  reference  table  are  used  for  fetching  quadruple  data  and  name 
list  data  of  systems  from  QDATA  and  NDATA  files. 

Modeling  data  consist  of  several  groups  of  data.  They  are  the  print 
specification  group,  the  system  reference  specification  group,  and  the 
system  description  group.  A  general  picture  of  modeling  data  is  given 
in  Figure  11. 


The  print  specification  group  specifies  the  printed  output  the  user  desires. 
The  various  print  specification  cards  and  their  descriptions  are  given  in 
Table  7.  Any  number  of  print  specification  cards  can  be  used.  If  no 
print  specification  cards  are  used,  a  default  specification  consisting 
of  PRINT  INPUT  DATA  and  PRINT  FINAL  OUTPUT  DATA  will  be  assumed 
by  the  program. 


Table  7.  Print  Specification  Cards 


Print  Specification  Card 

Result 

PRINT  NOTHING 

Nothing  is  printed 
on  the  line  printer 

PRINT  INPUT  DATA 

Input  data  cards 
will  be  printed 

PRINT  OUTPUT  DATA 

Output  data  for  each 
system  modelled  will 
be  printed  in  detail 

PRINT  FINAL  OUTPUT  DATA 

Only  the  state  space 
model  data  for  each 
system  will  be  printed 

PRINT  EVERYTHING 

Debugging  Information 
will  be  printed  in  addition 
to  the  above 

23 


The  system  reference  specification  group  specifies  the  system  description 
that  should  be  used  in  the  reference  table.  Figure  12  illustrates  the  system 
reference  data.  An  example  of  reference  specification  is  given  in  Section  VI. 

There  are  as  many  groups  of  system  description  cards  in  the  model  input 
data  as  there  are  systems  whose  state  space  model  is  being  computed  by 
the  program.  Figure  13  illustrates  the  system  description  data  for  a  system 
described  by  FLEXSTAB/LSA  equations.  The  FLEXSTAB/LSA  simulator 
deck  data  is  on  file  VDATA  and  only  the  header  card  is  needed  to  retrieve 
this  data.  The  name  list  data  arrays  for  the  system  variables  is  shown 
in  Table  8.  If  the  name  list  data  is  not  used  by  the  user,  a  default  name 
list  is  produced  by  the  program  and  is  useful  in  keeping  track  of  the  system 
variables. 

Figure  14  illustrates  the  system  description  input  data  for  transfer  functions. 
This  form  of  system  description  is  most  often  used  for  actuators,  sensors, 
and  controllers.  First  the  transfer  function  data  in  the  form  of  either 
rational  transfer  function  coefficients  or  transport  (time)  delay  parameters 
are  specified.  The  connection  data  for  the  transfer  function  blocks  are 
specified  next.  This  is  followed  by  name  list  data  in  the  end.  The  format 
for  the  transport  time  delay  parameters  are  given  in  Table  9.  An  example 
of  transfer  function  data  and  connection  data  specification  are  given  in 
Section  V. 

Figure  15  illustrates  the  system  description  input  data  when  data  are  in 
state  space  form  (quadruple  data).  First  the  quadruple  data  are  specified 
and  are  followed  by  the  name  list  data. 
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Table  8.  Formats  for  Name  List  Data 


CARD  SET  1  (12, 

6X,  2A4,  4X,  10A4,  4X,  4A4) 

Column 

Contents 

Explanation 

1-2 

NNS  ! 

Number  of  the  state  variable 

9-16 

VNS 

Name  of  the  state  variable 

21-60 

DESS 

Description  of  the  state  variable 

65-80 

UNITS 

Units  of  the  state  variable 

CARD  SET  2  (12, 

6X,  2A4,  4X,  10A4,  4X,  4A4) 

1-2 

NNO 

Number  of  the  output  variable 

9-16 

VNO 

Name  of  the  output  variable 

21-60 

DESO 

Description  of  the  output  variable 

65-80 

UNITO 

Units  of  the  output  variable 

CARD  SET  3  (12, 

6X,  2A4,  4X,  10A4,  4X,  4A4) 

1-2 

NNI 

Number  of  the  Input  variable 

9-16 

VNI 

Name  of  the  Input  variable 

21-60 

DESI 

Description  of  the  Input  variable 

65-80 

UNTI 

Units  of  the  Input  variable 

Table  9.  Formats  for  Transfer  Function  Data 


CARD  SET  1  (6E12.  6,  212) 

Column 

Contents 

Explanation 

1-12 

TD 

Transport  (or  Time)  delay  in  seconds 

13-24 

XX 

X-coordinate  value,  at  which  wind  gust 
is  acting,  in  feet 

25-36 

XR 

Reference  X-coordinate  value  in  feet 

S  37-48 

UU 

Velocity  of  wind  in  feet/sec 

49-60 

OMEGM 

Maximum  frequency  in  radians /sec 

61-72 

DELPHM 

Maximum  phase  error  in  radians 

73-74 

ND 

Number  of  Denominator  terms  in  Pade 
approximation. 

75-76 

NN 

Number  of  Numerator  terms  in  Pade 
approximation. 

Note:  The  rational  transfer  function  data  and  connection  data  are 

specified  as  matrix  data  and  are  always  read  using  the  format 
(5(212,  E12.6)).  Only  the  nonzero  elements  need  be  specified. 


Figure  16  illustrates  the  system  description  input  data  for  a  system  described 
as  an  interconnection  of  several  systems.  The  maximum  number  of  systems 
that  can  be  combined  into  one  system  is  nine.  The  interconnection  data, 
which  specifies  the  system  numbers  of  the  systems  being  interconnected, 
are  specified  as  a  matrix  data.  Only  the  nonzero  elements  need  be  specified. 
The  interconnection  data  are  followed  by  the  name  list  data.  An  example  of 
interconnection  of  several  systems  is  given  in  Section  V. 


Figui  e  17  illustrates  the  system  description  input  data  for  a  system  described 
as  the  conditioning  or  modification  of  a  system  whose  state  space  model 
has  already  been  obtained  by  the  program.  The  conditioning  data  consists 
of  scaling  data,  response  specification  data,  and  reduction  and  shuffling 
data.  The  formats  for  scaling  data  are  given  in  Table  10.  The  conditioned 
or  modified  system  gets  the  same  system  number  as  the  original  system 
and  will  be  the  current  system  in  the  system  reference  table.  An  example 
of  conditioning  a  system  is  given  in  Section  VI. 


Table  10.  Formats  for  Scaling  Data 


Note:  Name  of  the  scaled  variable  such  as  those  listed  in  Table  8 
should  be  of  the  following  form:  X(4),  R(8),  U(7)  etc. 


KONPACT-1  program  writes  the  state  space  model  (quadruple  data)  for  each 
system  on  QDATA  file.  It  uses  the  system  description  card  as  the  label 
to  store  data  for  several  systems  on  the  same  file.  This  label  card  is 
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used  for  retrieving  the  quadruple  data.  Similarly,  the  name  list  d^ta  for 
each  system  are  written  on  NDATA  file  using  the  system  description  card 
as  the  label.  In  addition,  the  system  reference  table  is  written  on  NDATA 
file  for  continuation  runs  (see  Figure  8).  Figure  18  shows  the  way  in  which 
data  is  stored  on  files  QDATA  and  NDATA. 

Print  Output  Data 

The  print  output  depends  on  the  print  specification  of  the  user.  The  following 
summarizes  the  sequence  of  print  output: 

•  Input  data  cards  are  printed  (if  PRINT  INPUT  DATA  is  specified). 

•  For  each  subsystem  defined,  the  input  data  to  the  system  is  printed 
after  the  system  heading  (if  PRINT  OUTPUT  DATA  is  specified). 

•  For  each  subsystem  defined,  the  state  space  data  (quadruple  data) 
and  the  name  list  data  are  printed  after  the  system  heading  (if  either 
PRINT  OUTPUT  DATA  or  PRINT  FINAL  OUTPUT  DATA  is  specified). 

PROGRAM  TIMING  AND  CENTRAL  MEMORY  REQUIREMENT  ESTIMATES 

Program  Timing  Estimate 


An  approximate  formula  for  the  central  processor  time  (TCp)  in  decimal 
seconds  is: 


T„  =  T 


Load  cp 


+  fl&±as  ±  NU)  *  ns 


where: 


^Load  cp  =  ProSram  Loading  Time  (typically  4  seconds) 

NX  =  number  of  states  of  the  system 
NR  =  number  of  outputs  of  the  system 
NU  =  number  of  inputs  of  the  system 
NS  =  number  of  systems  modeled. 

An  approximate  formula  for  Input/Output  Time  in  decimal  seconds  is: 

TI/C  ■  TLoad  I/O  +  (NX  +  NR  +  NU>  *  NS 
where: 

TLoad  I/O  =  Pro£ram  Loading  Time  (typically  20  seconds). 

Central  Memory  Requirement  Estimate 

The  precompiler  program  can  be  used  to  obtain  the  central  memory  required 
for  executing  the  job.  The  precompiler  data  was  discussed  earlier.  An 
example  job  set-up  is  shown  in  Figure  19.  When  this  job  is  executed,  the 
central  memory  required  to  execute  KONPACT-1  program  is  printed  out. 
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SECTION  IV 

DESIGN  PROGRAM  (KONPACT-2) 


INTRODUCTION 

KONPACT-2  program  interfaces  with  the  modeling  program  KONPACT-1 
via  the  data  file  QDATA  for  the  design  model.  KONPACT-2  program  also 
interfaces  with  the  linear  system  analysis  (LSA)  program  of  FLEXSTAB  via  the 
data  file  SDSTP  for  performance  evaluation  of  the  open  or  closed  loop  system. 
KONPACT-2  program  consists  of  three  subprograms  namely  DIAK,  FFOC, 
and  DATAK.  The  DIAK  subprogram  developed  by  Honeywell  (Reference  2) 
computes  the  optimal  state  feedback  gains.  The  FFOC  subprogram  also 
developed  by  Honeywell  (Reference  2)  computes  the  optimal  feedback  gains 
on  specified  measurements  (i.e. ,  Simplified  Controller  Design).  The 
DATAK  subprogram  developed  under  the  present  contract  reads  data  from 
data  files  and  data  cards  and  prepares  the  input  data  to  execute  DIAK  and 
FFOC  programs. 

The  DATAK  subprogram  also  computes  the  frequency  domain  representation 
of  the  open  loop  or  closed  loop  system  for  use  by  the  LSA  program.  The 
KONPACT-2  program  can  also  be  used,  independent  of  KONPACT-1  and 
LSA  programs,  as  a  stand-alone  program  for  computing  optimal  feedback 
gains. 

This  section  describes  the  execution  procedure  for  the  KONPACT-2 
program.  The  description  of  input  data  and  output  data  are  given  here. 


PROGRAM  EXECUTION  PROCEDURE 


Figure  20  shows  a  block  diagram  of  the  data  flow  in  KONPACT-2.  The 
DATAK  subprogram  reads  data  for  DIAK  and  open  loop  quadruple  data 
from  file  QDATA  (obtained  by  executing  KONPACT-1  program)  and  prepares 
the  data  for  DIAK  subprogram  on  a  scratch  file.  The  DIAK  subprogram 
reads  the  data  on  the  scratch  file  and  computes  optimal  state  feedback 
gains  and  writes  them  on  file  DDATA.  The  DATAK  program  then  reads 
FFOC  data,  open  loop  quadruple  data  on  file  QDATA,  and  gains  on  file 
DDATA  and  prepares  the  data  for  FFOC  subprogram  on  the  scratch  file. 

The  FFOC  subprogram  reads  the  data  on  the  scratch  file  and  computes  sub- 
optimal  feedback  gains  on  specified  measurement  and  writes  them  on  file 
FDATA.  Finally  the  DATAK  subprogram  reads  LSA  data,  open  loop 
quadruple  data  on  file  QDATA,  and  gains  on  file  FDATA  (it  can  also  read 
gains  on  file  DDATA  or  gains  on  cards);  it  computes  the  frequency  domain  repre¬ 
sentation  of  the  open  or  closed  loop  system  for  LSA  program  and  writes  it  on 
file  SDSTP.  The  data  on  files  DDATA  and  FDATA  are  formatted  whereas 
the  data  on  file  SDSTP  are  not  formatted.  The  job  control  cards,  used  for 
executing  KONPACT-2  programs,  are  those  commonly  used  by  the  SCOPE 
3.4.1  operating  system  for  CDC  6000  and  CDC  7000  series  computers. 

For  detailed  description  of  job  control  cards,  the  user  is  referred  to  the 
SCOPE  manual  (Reference  3).  The  data  files  of  KONPACT-2,  namely 
DDATA,  FDATA  and  SDSTP,  are  written  on  permanent  disc  files  which 
are  cataloged  at  the  end  of  the  job.  To  use  magnetic  tapes  instead  of  permanent 
disc  files,  the  user  is  referred  to  the  ASD  Computer  Center  User's  Guide 
(Reference  6), 
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Job  control  cards  to  execute  KONPACT-2  program  is  illustrated  in  Figure  21. 
The  input  data  cards  follow  the  job  control  cards.  In  Figure  21  it  is  assumed 
that  the  input  data  consists  of  data  for  DLAK  subprogram,  data  for  FFOC 
subprogram,  and  data  for  obtaining  frequency  domain  representation  of  the 
closed  loop  system  for  LSA  program.  In  practice,  however,  several  iterations 
are  made  with  DIAK  data  alone  and  subsequently  with  FFOC  data  alone  to 
obtain  final  feedback  gains.  DIAK  programs  can  also  be  used  for  covariance 
analysis  alone  or  time  response  alone  for  the  final  closed  loop  system.  The 
user  must  attach  ?.nd  catalog  the  required  data  files. 

INPUT  DATA 

Figure  22  illustrates  the  card  input  data  for  KONPACT-2  program.  It 
consists  of  four  different  groups  of  data.  They  are  the  print  specification 
data,  the  DIAK  input  data,  the  FFOC  input  data,  and  the  LSA  input  data. 

The  print  specification  data  is  described  in  Section  III.  In  this  program, 
the  print  specifications  control  only  the  printed  output  from  DATAK  sub¬ 
program  and  have  no  control  over  the  printed  output  from  the  DIAK  and  the 
FFOC  subprograms.  Modifications  in  DIAK  and  FFOC  were  beyond  the 
scope  of  the  contract. 

The  DIAK  input  data  and  the  FFOC  input  data  are  shown  in  Figure  23  and 
Figure  24,  respectively.  The  sequence  and  meaning  of  data  cards  for  DIAK 
and  FFOC  programs  as  documented  in  Reference  2  are  retained  fully.  Under 
KONPACT-2  program,  comment  cards  can  be  inserted  anywhere  in  the  data; 
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options  are  provided  to  read  certain  matrix  data  from  QDATA,  DDATA 
and  FDATA  files  in  addition  to  reading  them  from  card  input  data.  The 
options  for  various  matrix  data  in  the  KONPACT-2  program  are  shown  in 
Table  11.  These  options  are  selected  by  the  use  of  data  control  cards. 

The  list  of  data  control  cards  is  given  in  Table  12.  The  user  should  consult 
Reference  2  for  the  details  of  the  input  data  for  the  DIAK  and  the  FFOC 
subprograms.  Several  examples  are  given  in  Section  VI. 

The  LSA  input  data  (to  compute  frequency  domain  data  for  LSA  program) 
are  shown  in  Figure  25.  The  open  loop  quadruple  data  are  input  from 
QDATA  file.  Option  is  provided  to  read  the  feedback  gains  from  DDATA 
or  FDATA  file  or  from  card  input  data.  If  feedback  gains  are  provided, 
the  KONPACT-2  program  computes  the  frequency  domain  representation 
for  the  closed  loop  system.  If  no  feedback  gains  are  provided,  the  KONPACT-2 
program  computes  the  frequency  domain  representation  for  the  open  loop 
system. 

OUTPUT  DATA 


File  Output  Data 

Optimal  state  feedback  gain  matrix  data  are  written  on  file  DDATA  when 
DIAK  subprogram  is  executed.  Optimal  reduced  feedback  gain  matrix  data 
are  written  on  file  FDATA  when  FFOC  subprogram  is  executed.  Frequency 
domain  data  are  written  on  file  SDSTP  when  DATAK  subprogram  is  executed. 

The  data  on  files  DDATA  and  FDATA  are  formatted,  and  the  data  on  file 
SDSTP  are  not  formatted. 


33 


'  ■'t’H'vM 


1  ..jy 


'Tin — mmmmm  m  itiri^m'Trnmn- 


Table  11.  Options  for  Obtaining  Matrix  Data  for  the  KONPACT-2  Program 


XXXXXXXXXXXXXXXX 


.a  . . 

b  • a 


4-»  ^ 

m  .3  S 


(0  3 

S  9 


rt 

§  cti 


a  c 
c  o 


1  I  s 

g  5  9 

^  O  tJ 

id  i  s 

&  6  .2 

co  O  Q 


^  2 

«  $ 

I  I 

3  rt 
W  +3 

o  3 


a  .a 

s  § 

J  £ 

0)  r-t 

**  3 

a  > 

K  4) 

■*■>  i-3 


jS  o 
co  U 


n 

t  |  .s 

«  s  E 

s  s  « 

,  8  * 

3  8.  s 

8,  I  e 

m  K  k 

«  *  2 

K  o  3 

0)  iJ  » 

*»  C  cfl 

rt  o  <u 

«  u  2 


.a 

S  -g 

a  S 

.s  s 

3)  3 

a  ■& 

CO  d) 

«  is 

$  a 
ft  5 

.2  ■o 


o 

I  7 

4)  . 

dj 

ft  In 
4)  «2 

3  s 

co  .5 

t-t  ?* 

rt  W 

§T3 
4) 

S'  * 

O  ft 


ST) 
4) 

■a  .a 

O  ft 


,4 

8 

CM  _  J  1-1 


CM  t— I 

y  X  X  X 
X  X  X  2 


ftooxxuwp<pqcy<;<cQ<; 


DELK(AK1(X))  Predictor  gains  for  X 


Table  12.  List  of  Data  Control  Cards  for  the  KONPACT-2  Program 
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Table  12.  List  of  Data  Control  Cards  for  the  KONPACT-2  Program  (Continued) 


Print  Output  Data 


The  following  items  summarize  the  print  output. 

•  Card  input  data  are  printed  (if  PRINT  INPUT  DATA  is  specified). 

•  Output  data  from  the  DIAK  subprogram  are  printed. 

•  Output  data  from  the  FFOC  subprogram  are  printed. 

•  Frequency  domain  data  are  printed. 

PROGRAM  TIMING  AND  CENTRAL  MEMORY  REQUIREMENT* 

ESTIMATES 

It  is  very  difficult  to  estimate  the  time  required  to  execute  KONPACT-2 
program  because  of  the  several  options  provided  in  the  DIAK  and  FFOC 
subprograms. 

A  central  memory  of  165000  octal  is  required  to  execute  KONPACT-2 
program.  This  value  is  based  on  the  maximum  system  dimensions  presently 
used  in  program  KONPACT-2.  The  values  of  the  maximum  system  dimensions 
and  the  main  programs  in  which  they  are  defined  are  shown  in  Table  13. 

For  changing  the  maximum  system  dimensions,  the  user  should  consult 
Reference  2. 
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Table  13.  Maximum  System  Dimensions  Used  in  KONPACT-2 
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SECTION  V 


DEMONSTRATION  EXAMPLE  1 


INTRODUCTION 


The  example  chosen  in  this  section  for  demonstration  is  the  design  of  a 
handling  qualities  controller  for  the  C-5A  basic  vehicle  (cruise  flight 
condition).  This  model  was  generated  from  the  Air  Force  Level  2.01.00 
FLEXSTAB  Computer  Programming  System  (Reference  1). 

In  the  following  subsection,  the  problem  statement  and  the  method  of 
solution  are  discussed  first.  Then  the  deck  set-up  for  solving  the  problem 
using  KONPACT  is  explained.  Finally,  the  computer  printout  and  the 
final  design  results  are  discussed. 

It  is  assumed  that  the  user  is  now  familiar  with  the  optimal  control 
programs  (DLAK  and  FFOC)  and  the  design  procedures  listed  in  Reference  2. 

PROBLEM  STATEMENT 

The  problem  is  to  design  a  controller  for  the  C-5A  basic  vehicle  (cruise 
flight  condition)  so  that  the  controlled  vehicle  will  have  the  desired  handling 
quality  (Reference  10).  The  handling  quality  criterion  in  this  example 
is  prescribed  in  terms  of  the  short  period  poles.  The  desired  short 
period  poles  are  given  as:  <un  =  3  rad/sec  and  £  =  0. 7. 
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METHOD  OF  SOLUTION 


The  C-5A  handling  quality  controller,  designed  via  the  optimal  control 
synthesis  technique,  is  carried  out  in  three  parts:  a)  design  model 
generation,  b)  controller  design,  and  c)  performance  evaluation. 

Design  model  generation  consists  of  three  steps.  The  first  step  is  to 
obtain  the  required  dynamics  (i.  e. ,  vehicle,  actuator,  pilot  model, 
implicit  model)  needed  for  the  C-5A  handling  quality  design  process.  The 
second  step  is  to  augment  the  C-5A  vehicle  model  with  appropriate  pilot 
model,  actuator,  and  implicit  model  (see  Figure  26).  The  third  step  is 
to  select  the  design  responses  and  obtain  a  design  model. 

Controller  design  consists  of  two  steps.  The  first  step  is  to  obtain 
optimal  state  feedback  controller  gains  by  varying  the  quadratic  weights 
on  selected  design  responses  so  that  the  handling  quality  criterion  is  met. 

The  second  step  is  to  obtain  reduced  feedback  gains  so  that  only  a  specified 
set  of  sensors  is  used,  and  the  handling  quality  criterion  is  not  compromised. 

Performance  evaluation  of  the  closed  loop  system  is  obtained  to  check  the 
design. 

Design  Model  Generation 

Figure  26  shows  a  schematic  diagram  of  the  design  model.  In  the 
following  paragraphs,  the  steps  for  obtaining  the  design  model  are 
described. 
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C-5A  Basic  Vehicle  for  Cruise  Flight  Condition  (Static  Elastic  Symmetric 
Dynamic  Math  Model) --The  C-5A  basic  vehicle  for  cruise  flight  condition 
(static  elastic  symmetric  model)  is  provided  by  the  FLEXSTAB  LSA 
program  in  the  form  of  a  simulator  deck  data. 

This  model  is  given  system  number  1  (S-l,  Figure  26)  and  STAMK1 
subprogram  in  KONPACT-1  is  used  to  obtain  the  state  space  model  from 
the  simulator  deck  data.  The  system  variables  are  the  following: 

States:  u,  w,  q,  0 

Outputs:  q  ,  n  ,  n  ,  n 

s  aj  a2  a3 

Inputs:  6a.  6e.  6^.  «e<  Wgj.  wg2>  wg3>  *g2>  wgg. 

WSsl  W*s2'  WSs3'  *«.!■  *Ss2'  *gS3 

Table  14  gives  their  description  and  also  forms  the  basis  for  the  name 
list  data  input. 

C-5A  Reduced  Vehicle  for  Cruise  Flight  Condition  (Reduced  Static 

Elastic  Symmetric  Dynamic  Math  Model)- -The  handling  quality  controller 

design  requires  only  the  short  period  dynamics  (i.  e. ,  u>  and  q),  input, 

and  q  and  n„  outputs.  The  reduced  model  is  obtained  by  truncating 
s  “1 

the  states  u  and  0  and  omitting  the  inputs  and  outputs  not  needed.  Also, 
the  states  w  and  q  are  defined  as  outputs.  This  reduced  model  is  obtained 
by  using  the  CONDK  subprogram  in  KONPACT-1.  The  system  number  1 
is  retained  for  the  reduced  system  (S-l,  Figure  26).  The  system  variables 
in  this  case  are: 

( 
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Table  14.  Description  of  the  System  Variables  for  the  C-5A  Static  Elastic  Vehicle 
(Cruise  Flight  Condition) 
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States:  w,  q 

Outputs:  q  ,  na  ,  w,  q 

®  1 

Inputs:  $e 


Actuator  Model- -The  actuator  is  given  by  the  transfer  function: 
b 


H  (s)  =  -f- 
a  s  +  a 


;  a  =  b  =  7. 5. 
a  a 


(3) 


The  actuator  is  given  a  system  number  2  (S-2,  Figure  26)  and  is  modeled 
using  the  STAMK2  subprogram  in  KONPACT-1. 


Pilot  Model — The  pilot  model  is  a  first  order  lag  filter  given  by  the  transfer 
function: 

bn  * 

HP(S)  *  a  +T  1  ap  =  °'1,  bo  s  0,223  X  10  •  (4) 

P  p  P  P 

The  pilot  input  time  constant  is  chosen  as  10  seconds  since  this  gave  good 

GLA  (Gust  Load  Alleviation)  performance  when  used  with  rate  model 

-3 

following  (Reference  10).  The  gain  b  is  chosen  to  be  0.  223  x  10  in 

P  _  2 

order  to  obtain  an  rms  pilot  input  of  0,  5  x  10  radian  (0. 0287  degree). 

This  is  an  important  factor  in  reduced  (not  full)  control  since  optimal 
gains  depend  on  covariance  analysis. 

The  pilot  model  is  given  a  system  number  3  and  is  also  modeled  using  the 
STAMK2  subprogram  in  KONPACT-1. 


Interconnection  of  Reduced  Vehicle,  Actuator,  and  Pilot  Model  into  a 


System  Called  Plant--The  interconnections  for  the  plant  are  shown  in 
Figure  26.  The  interconnections  are  described  by  the  following  equations: 


UI1/RI1  I  UI1/RI2  |  UI1/RI3 

"l 

UI2/RI1  !  UI2/RI2  I  UI2/RI3 


UI3/RI1  '  UI3/RI2  !  UI3/RI3 


il 


i2 


i3 


+ 


(r)  =  (R/RI1  !  R/RI2 


R/RI3) 


(R/U)  (U) 


(5) 


(6) 


where  U.^,  U.2#  and  Ui3  are  the  inputs  to  systems  1,  2,  and  3,  respectively, 
and  r.j,  ri2,  and  r  are  the  outputs  from  systems  1,  2,  and  3,  respectively. 
U  is  the  input  to  the  plant,  and  r  is  the  output  from  the  plant.  Definitions 
of  the  inputs  and  outputs  in  terms  of  the  system  physical  variables  are 
given  in  Table  15. 


Jr 

I 
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Table  15.  Input/Output  Definitions  for  Obtaining  Static  Elastic  Plant  Model 


System 

No. 

Name 

Input  Definition 

Output  Definition 

1 

Reduced  Vehicle 

Uil  ’  8e 

ril  =  Co1  <Vnal'w*q} 

2 

Actuator 

Ui2  =  6ec 

ri2  =  6e 

3 

Pilot  Model 

Ui3  ’  % 

ri3  =  rp 

4 

Plant 

U  =  Col  {!)  ,u  } 
p  c 

r  =  Col  tqs.nal,  w,q,rp} 

Thus  the  non-null  interconnection  matrices  are  given  in  the  following: 


UI1/RI2  =  (1) 

UI2/U  =  (0  1) 


UI3/U 

R/RI1 


R/RI3 


(7) 
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The  interconnection  data  are  used  by  subprogram  STAMK3  of  KONPACT-2 
to  model  the  plant.  The  plant  is  given  a  system  number  4  (S-4,  Figure  26). 

Implicit  Model- -The  implicit  model  is  used  to  generate  model  following 
error  rate  responses  (Reference  10).  These  model  following  error  rate 
responses  are  weighted  during  the  design  step  to  achieve  handling  quality 
criteria. 

The  implicit  model  construction  is  described  in  the  following. 

The  reduced  C-5A  vehicle  is  described  by  the  equations: 


*  *  AUw+A12q+BU  6e 
q  =  A21w+A22q+B21  6e 


where  the  values  of  the  coefficients  for  cruise  flight  condition  are  given  by: 


-0.67868  8741.2 


B  = 


0. 33079  1 
-1, 6064  j 


\  -0.0001874  -1.1011 j  \ -1.606 

The  implicit  vehicle  is  described  by  the  equations : 


W,/J 

M 

1  AuwM 

+  A12q  M 

+Bn 

U  M 

^  M 

A 

=  A21wM 

A 

+  A22qM 

+  B12 

u  M 

where  wM  and  qM  are  state  variables  of  the  implicit  model,  and  uM 
is  the  input  to  the  implicit  model. 


As  can  be  seen,  the  coefficients  of  the  first  row  of  the  implicit  model 
transition  matrix  are  the  same  as  the  corresponding  elements  of  the 
vehicle  equation.  The  second  row  coefficients  are  obtained  from  the 
handling  quality  criteria,  the  prescribed  short  period  poles.  These  are 

given  as: 

u )  =3  rad/sec  and  C  =  0*7. 

n 


I 

; 

■■ 


1 

f 

Jt 


i 


The  corresponding  characteristic  equation  is: 


2  2 
Ms)  *  s  +  2  C«n  s  +  «»n 


(ID 


since: 


tr(A)  =  -2  Cu>n  -  -4.2 

det(A)  -  uj  2  *  9 
n 


(12) 


and: 


All+A22=  “4-2 


A11  A22  "  A12  A21  *  9 


(13) 


V 

m 


i 


and  the  results  are: 

A 

A  *  -.0007562  and 
Ag2  s  ”3. 5213, 


(14) 


The  implicit  model  is  always  given  the  system  number  9  (S-9,  Figure  26). 
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Interconnection  of  the  Plant  and  the  Implicit  Model  (Overall  System)-- 
The  interconnections  for  the  overall  system  are  shown  in  Figure  26. 
Table  16  gives  input/output  definitions  for  this  case.  The  non-null  inter¬ 
connection  matrices  are  given  by: 


UI9/RI4 


UI4/U 


R/RI4 


R/RI9 


(0  0  0  0  1) 


0  1 


0  l/u0  0 


The  interconnection  data  are  used  by  subprogram  STAMK3  of  KONPACT-1 
to  model  the  overall  system.  The  overall  system  is  given  a  system 
number  5  (S-5,  Figure  26). 


Table  16.  Input /Output  Definitions  for  Obtaining  Overall 
Static  Elastic  Syslem  Model 


System 

No. 

Name 

Input  Definition 

Output  Definition 

4 

Plant 

u  *  Col  {T)  ,  u  } 
i4  p  cJ 

ri4  ’  Co1  <V  nal ’  »•  <»•  rp5 

9 

Implicit 

Model 

u.  n  =  u.. 
i9  M 

ris'CollwM  ’M  5 

5 

Overall 

System 

u  =  Col  {T|  ,  u  ) 
p  c 

r  -  Col  (qs.  nal.  rp.  e^.  eq.) 

Overall  System  (Design  Model) --The  design  response  specifications, 
sensor  specifications,  control  input  specifications,  and  disturbance  input 
specifications  are  used  by  the  CONDK  subprogram  of  KONPACT-1  to  obtain 
the  design  model.  System  number  5  is  retained  for  the  design  model 
(S-5,  Figure  26). 


The  specified  design  responses  for  this  example  are: 

/">■  “•  V  V  V  V 


The  specified  sensor  responses  are: 


n 


al’ 


The  specified  control  inputs  are: 

<uc). 
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Ihe  specified  disturbance  inputs  are 


he  implicit  monel  is  automatically  truncated  by  the  CONDK  subprogram 
a  ter  mode,  following  error  rate  equations  are  obtained  internally.  (For 
an  exp  icit  model,  a  system  number  other  than  9  should  be  assigned. ) 

ZtZT7 needed  by  the  subprograms  DMK- FFoc-  -  *» 

ONPACT-2  are  given  by  the  following  equations  (see  Reference  i), 


X  = 


(F)x  +  (G^u  +  (G  )T) 


r  =  (H)x  +  (D)u 

m  =  (AM)x 


(16) 


where: 

x  is  the  state  variables 
u  is  the  control  variables 
^  is  the  disturbance  variables 
i  is  the  design  responses 
m  is  the  measurements. 


The  quadruple  data  (A.  B.  C.  D)  representing  the  design 
the  previously  mentioned  matrices  as  shown  below: 


model  contain 


A 

C 


(17) 
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Controller  Design 

Figure  27  represents  a  block  diagram  of  the  design  process.  The  starting 
design  response  weights  are  shown  in  Figure  28.  In  the  following  paragraphs, 
the  steps  for  obtaining  the  handling  quality  controller  design  are  described. 

Optimal  State  Feedback  Gains  for  Handling  Quality  Controller  Design- - 
Using  the  overall  system  design  model  described  earlier,  the  optimal  state 
feedback  gains  are  computed  by  the  DLAK  subprogram  of  KONPACT-2  for 
the  specified  quadratic  weights.  The  quadratic  weights  are  varied  and  the 
optimal  state  feedback  gains  are  computed  until  the  handling  quality  criterion 
is  met.  The  variation  of  quadratic  weights  and  the  resulting  system  gains 
are  shown  in  Figure  29. 

Reduced  Optimal  Feedback  Gains — The  optimal  full  state  feedback  control 
law  is  given  by: 

6  =  (K  )q  =  (K  )n  .  +  (K.  )  «  +  (K  >  r  (1C) 

ec  qg  s  nal  al  6e  e  rP  P 

Reduced  optimal  feedback  gains  are  obtained  for  two  cases  by  the  FFOC 
subprogram  (see  Figure  30).  In  the  first  case,  the  gain  on  measurement 
alone  is  reduced  to  zero.  In  the  second  case,  the  gains  on  measurements 
and  ngj  are  reduced  to  zero.  The  effect  of  these  reductions  is  presented 
in  Figure  31.  The  first  reduction  maintains  the  desired  C-5A  handling  qualities. 

Performance  Evaluation 

In  Figures  32  through  34  the  time  responses  of  the  open  and  closed  loop  system 
are  obtained  to  demonstrate  the  desired  handling  quality.  These  time  responses 
are  due  to  the  elevator  input  contained  in  each  figure.  The  short  period  pole 
locations  are  also  checked  out  which  is  demonstrated  in  Figure  30. 
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DECK  SET-UP 

Figure  35  shows  the  static  elastic  simulator  deck  data  which  are  calculated 
by  the  FLEXSTAB/LSA  program.  The  precompiler  data  needed  for  KONPACT-1 
are  shown  in  Figure  36.  It  consists  of  the  subprograms  that  will  be  needed  to 
execute  the  KONPACT-1  input  data.  The  input  data  to  produce  the  design 
model  of  the  overall  system  are  shown  in  Figure  37.  The  KONPACT-2 
input  data  required  to  obtain  optimal  state  feedback  gains  (using  the  DIAK 
subprogram)  are  shown  in  Figure  38.  The  KONPACT-2  input  data  required 
to  compute  the  feedback  gains  only  on  specified  measurements  (using  the 
FFOC  subprogram)  are  shown  in  Figure  39.  The  KONPACT-2  input  data 
required  to  obtain  time  response  (using  the  DIAK  subprogram)  are  shown 
in  Figure  40. 

OUTPUT  DESCRIPTION 

KONPACT-1  output  data  are  shown  in  Figures  41  through  53.  The  corresponding 
input  data  are  shown  in  Figure  37.  The  C-5A  vehicle  simulator  deck  data, 
for  cruise  flight  condition  in  matrix  form,  are  given  in  Figure  41,  and  the 
corresponding  vehicle  quadruple  data,  along  with  the  name  list  data,  are 
given  in  Figure  42,  The  reduced  vehicle  name  list  data  and  quadruple  data 
are  given  in  Figure  43.  The  actuator  transfer  function  data  are  given  in 
Figure  44,  and  the  corresponding  actuator  quadruple  data  along  with  the 
name  list  data  are  given  in  Figure  45.  The  pilot  model  transfer  function 
data  are  given  in  Figure  46  and  the  corresponding  pilot  model  quadruple 
data,  along  with  the  name  list  data,  are  given  in  Figure  47.  The  interconnection 
data  for  the  plant  are  given  in  Figure  48,  and  the  corresponding  plant 
quadruple  data,  along  with  the  name  list  data,  are  given  in  Figure  49.  The 


implicit  model  name  list  data  and  the  quadruple  data  are  given  in  Figure  50. 
The  interconnection  data  for  the  overall  system  are  given  in  Figure  51, 
and  the  corresponding  overall  system  quadruple  data,  along  with  the  name 
list  data,  are  given  in  Figure  52.  The  overall  system  design  model  name 
list  data  and  quadruple  data  are  given  in  Figure  53. 

KONPACT-2  output  data,  obtaining  optimal  state  feedback  gains  using 
DIAK  subprogram,  are  shown  in  Figure  54.  Even  though  the  corresponding 
input  data  given  in  Figure  38  contain  five  different  quadratic  weighting 
matrices,  the  output  data  shown  here  contain  the  iterations  for  the  fourth 
quadratic  weighting  matrix.  The  results  for  the  other  quadratic  weighting 
matrices  are  summarized  in  the  next  subsection. 

KONPACT-2  output  data,  for  obtaining  feedback  gains  on  specified 
measurements  using  the  FFOC  subprogram,  are  shown  in  Figure  55. 

KONPACT-2  output  data,  for  obtaining  time  response  of  closed  loop  system 
using  the  DIAK  subprogram,  are  shown  in  Figure  56. 

DISCUSSION  OF  RESULTS 

In  the  first  step  of  the  design  optimal  gains,  K^g,  Kj^,  Kge,  ant*  ^rp  are 
determined  by  varying  the  quadratic  weights  and  using  the  DIAK  subprogram 
until  the  handling  quality  criterion  is  met.  Figure  29  shows  the  variation 
in  gains  and  the  closed  loop  eigenvalues  when  the  weight  Qg  is  varying 
from  .001  to  10.0.  The  movement  of  the  closed  loop  eigenvalues  is 
also  shown  in  Figure  31, 


I 
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In  the  second  step  of  the  design,  the  number  of  feedback  gains  is  reduced 
by  using  the  FFOC  subprogram.  When  both  K^e  and  are  reduced,  the 
resulting  closed  loop  eigenvalues  do  not  meet  the  handling  quality  criterion. 
When  only  K^e  is  reduced,  the  handling  quality  criterion  is  met.  The 
variation  in  gains  and  the  closed  loop  eigenvalues  when  the  reduction  of 
gains  takes  place  are  shown  in  Figure  30.  The  movement  of  the  closed 
loop  eigenvalues  is  also  shown  in  Figure  31. 

The  step  response  for  q,  a,  and  &e  for  the  open  loop  system  is  shown  in 
Figure  32.  The  step  response  for  q,  a,  and  6  for  the  full  state  optimal 
closed  loop  system  is  shown  in  Figure  33.  The  step  response  for  q,  a, 
and  6^  for  the  reduced  (only  K5g  is  reduced)  optimal  closed  loop  system 
is  shown  in  Figure  34.  As  can  be  seen,  the  reduced  feedback  optimal 

system  develops  responses  similar  to  the  full  ~»cate  optimal  system. 
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SECTION  VI 


DEMONSTRATION  EXAMPLE  2 


INTRODUCTION 

The  example  chosen  in  this  section  for  demonstration  is  the  repeat  design 
of  the  Active  Lift  Distribution  Control  System  (ALDCS)  for  the  C-5A 
vehicle  (cruise  flight  condition)  per  Reference  10.  The  vehicle  model 
was  generated  from  the  Air  Force  Level  2.01.00  FLEXSTAB  Computer 
Programming  System  (Reference  1). 

The  problem  statement  and  the  method  of  solution  are  discussed  first. 

Next,  the  deck  set-up  for  solving  the  problem  using  KONPACT  is  explained. 
Finally,  the  computer  printout  and  the  final  design  results  are  discussed. 

PROBLEM  STATEMENT 

ALDCS  design  goals  are  shown  in  Table  17.  The  problem  is  to  design 
a  controller  for  the  C-5A  vehicle  (cruise  flight  condition)  to  meet  these 
goals. 

METHOD  OF  SOLUTION 

The  C-5A  ALDCS  controller  design  via  the  optimal  control  synthesis 
technique  is  carried  out  in  three  parts:  a)  design  model  generation, 
b)  controller  design,  and  c)  performance  evaluation. 
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Design  Model  Generation 

Figures  57a  and  57b  show  schematic  diagrams  of  the  design  model.  In  the 
following  paragraphs,  the  steps  for  obtaining  the  design  model  are  described. 


C-5A  Vehicle  for  Cruise  Flight  Condition--The  C-5A  vehicle  for  cruise 
flight  condition  (Residual  Elastic  Symmetric  Model)  is  provided  by  the 
FLEXSTAB  (LSA  program)  in  the  form  of  a  simulator  deck  data  written 
on  file  VDATA.  A  full  description  of  the  model  is  presented  in  References 
8  and  9  and  therefore  will  not  be  repeated  in  this  report.  This  model  is 
given  system  number  1  (S-l,  Figure  57a),  and  the  STAMK1  subprogram  in 
KONPACT-1  is  used  to  obtain  the  state  space  model  from  the  simulator 
deck  data.  For  a  definition  of  the  load  reference  axis  systems  used  to 
define  (Bi,  Ti  i  =  1,  ...15),  see  Figure  4  of  Reference  4.  The  system 
variables  are  the  following: 

States:  u,  w,  q,  8,  Tj^,  T]^,  ...  ,  ^5*  ^2*  *  *  *  *  ^15 

Outputs:  qg,  nftl,  n^,  SI,  Bl,  Tl,  S2,  B2,  T2,  ...  , 

S5,  B5,  T5 

Inputs:  6a.  6ei,  6^,  «a>  6^,  wgl<  wg2>  wg3>  *gl>  *g2>  wg3 


Table  18  describes  these  symbols  and  forms  the  basis  for  the  namelist  data. 
Note  that  a  redefinition  of  generalized  coordinates  has  been  made  from 


"UjjM  to  "Tl."  to  facilitate  later  discussion. 


Table  18.  Description  of  the  System  Variables  for  C-5A  Residual 
Elastic  Vehicle  (Cruise  Flight  Condition) 
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Table  18.  Description  of  the  System  Variables  for  C-5A  Residual 
Elastic  Vehicle  (Cruise  Flight  Condition)  (Continued) 
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Table  18.  Description  of  the  System  Variables  for  C-5A  Residual 
Elastic  Vehicle  (Cruise  Flight  Condition)  (Continued) 
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Table  18.  Description  of  the  System  Variables  for  C-5A  Residual 
Elastic  Vehicle  {Cruise  Flight  Condition)  (Concluded) 
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-C-5A  Converted  Vehicle  for  Cruise  Flight  Condition--It  is  assumed  that  the 

ALDCS  controller  design  requires  only  the  short  period  dynamics  (i.e.,  w 

and  q)  and  the  flexure  mode  states  (7]  ,  ...  ,  T)  ,  T|, ,  ...  ,  1,).  In  this 

1  151  15 

case  the  assumption  is  valid,  as  the  unaugmented  C-5A  has  substantial 
separation  between  its  short  period  and  plugoid  modes.  For  other  aircraft 
of  reduced  static  stability  that  do  not  possess  this  eigenvalue  separation, 
the  assumption  is  then,  of  course,  not  valid.  Also,  the  design  requires 
that  units  of  some  of  the  variables  need  to  be  changed  (see  Table  19). 

The  C-5A  converted  vehicle  is  obtained  by  truncating  the  states  u  and  6  and 


Table  19.  Change  of  Units  for  Some  of  the  System  Variables  of  C-5A 
Residual  Elastic  Vehicle  (Cruise  Flight  Condition) 


Variable 

Old  Unit 

New  Unit 

Scale  Factor 

q 

Radian/Sec 

Inch /Sec 

0 . 164789E  04 

w  . 

gl 

Inch /Sec 

Feet/Sec 

0. 833333E-01 

Wg2 

Inch/Sec 

Feet/Sec 

0. 833333E-01 

Wg3 

Inch /Sec 

Feet/Sec 

0. 833333E-01 

V 

Inch/ Sec 

Feet/ Sec 

0. 833333E-01 

Wg2 

Inch /Sec 

Feet/Sec 

0 . 833333E-01  ; 

**3 

Inch /Sec 

Feet/Sec 

0. 833333E-01 

nal 

Inch /Sec2 

1G 

0 . 258800E-02 

na2 

Inch /Sec2 

1G 

0 . 258800E-02 

na3 

j  Inch /Sec2 

1G 

0 . 258800E-02 
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changing  the  units  of  the  variables  specified.  System  number  1  is  retained  for 
the  converted  vehicle  (S-l,  Figure  57a)  and  is  obtained  using  the  CONDK 
subprogram. 

Actuator  Model--The  actuator  block  diagram  is  shown  in  Figure  58.  The 
actuator  is  given  a  system  number  2  (S-2,  Figure  57a)  and  is  modeled 
using  the  STAMK2  subprogram  in  KONPACT-1. 

Gust  Model — The  gust  model  block  diagram  is  shown  in  Figure  59,  and  the 
gust  model  (Reference  10)  coefficients  are  shown  in  Table  20.  The  gust 
model  is  given  a  system  number  3  (S-3,  Figure  57a)  and  is  modeled  using 
the  STAMK3  subprogram  in  KONPACT-1. 


Table  20.  Gust  Model  Coefficients 


Coefficient 

Description 

Value 

/V 

!  ai 

Wind  filter  coefficient 

0.42 

a 

o 

Wind  filter  coefficient 

0.0441 

Wind  filter  coefficient 

0.7937 

/V 

b 

o 

Wind  filter  coefficient 

-0.2371 

VNT 

Kussner  coefficient  for  nose  and  tail 

9.0284 

V 

w 

Kussner  coefficient  for  wing 

4.4596 

T 

w 

First  order  time  delay  filter  coefficient 

0.0744 

A 

ai 

Second  order  time  delay  filter  coefficient 

16.  0044 

A 

a 

o 

Second  order  time  delay  filter  coefficient 

96.0523 

A 

bl 

Second  order  time  delay  filter  coefficient 

-8.0022 

A 

b 

o 

Second  order  time  delay  filter  coefficient 

224.1219 

63 


Interconnection  of  Converted  Vehicle,  Actuator  and  Gust  Model  (Plant) — 

The  interconnections  for  the  plant  are  shown  in  Figure  57a.  The  non-null 
interconnection  matices  for  these  interconnections  are  obtained  as  in 
Section  V.  The  interconnection  data  are  used  by  subprogram  STAMK3  of 
KONPACT-1  to  model  the  plant.  Two  additional  responses  a  and  q ^  are 
defined  here  for  the  plant.  The  plant  is  given  a  system  number  4  (S-4, 

Figure  57a). 

Plant  Design  Model — The  load  rate  responses  (B^,  T^,  Bg.  Tg*  •••  *  ®5»  ^5) 
are  obtained  here  by  using  the  CONDK  subprogram  of  KONPACT-1.  The 
system  number  4  is  retained  for  the  plant  design  model  (S-4,  Figure  57a). 


Reduced  Order  Plant  Models — The  plant  obtained  in  the  preceding  step  con¬ 
tains  15  flexure  modes.  This  is  reduced  to  the  six  lowest  frequency  flexure 
modes  by  residualization  and  truncation  procedures  (Reference  4),  Three 
different  models  are  obtained.  The  procedure  for  obtaining  the  three  models 
is  shown  in  Table  21.  The  system  reference  specification  is  used  to  obtain 


Table  21.  Reduced  Order  Plant  Model  (Six  Modes)  Obtained  from 
Residual  Elastic  Plant  Design  Model  (15  Modes) 


Name 

Procedure 

F24RR 

Obtained  by  residualizing  states,  design 

responses,  and  sensor  measurements. 

j  F24RT 

Obtained  by  residualizing  states  and  design  \ 

responses  and  truncating  sensor  measure- 

ments. 

F24TT 


Obtained  by  truncating  states,  design 
responses  and  sensor  measurements. 


the  three  models  from  the  plant  design  model  (see  Figure  79).  The  F24RR 
model  is  used  for  the  ALDCS  controller  design. 


ALDCS  Controller  Model- -The  block  diagram  of  the  ALDCS  controller 
is  shown  in  Figure  60.  This  model  is  given  system  number  5  (S-5,  Figure  57b). 
The  STAMK4  subprogram  in  KONPACT-1  models  the  controller  dynamics 
using  the  user  written  subroutine  SIMX2.  The  program  listing  for  SIMK2  for 
the  ALDCS  controller  is  given  in  Figure  61.  The  system  variables  for 
the  ALDCS  controller  are  the  following: 

States:  (A21RL,  GLAF,  MLC1,  HP,  F3E,  MLC2,  P) 

Outputs:  (u  ,  u  u.  ) 

K  6a  6ei  6eo 

Inputs:  (u  u  ,  71  ,  6  ,  6  .,  A21R,  AFUS,  TFUS,  ACG). 

c  j.  oo  p  a  6i 

The  gains  KAF  and  KM1  for  the  ALDCS  controller  are  obtained  by  conducting 
three  steady  state  tests  (Reference  4)  on  the  corresponding  plant  model 
(F24RR).  The  steady  state  tests  made  are  given  in  Table  22.  These  tests 
are  described  more  fully  in  Reference  4. 

Reduced  ALDCS  Controller- -The  general  ALDCS  controller  modeled  in  the 
previous  step  is  reduced  by  truncating  the  states  (HP,  F3E,  MLC2).  The 
block  diagram  of  the  reduced  controller  is  shown  in  Figure  62.  The  system 
number  5  is  retained  for  the  reduced  ALDCS  controller  (S-5,  Figure  57) 
and  is  obtained  using  the  CONDK  subprogram. 

Interconnection  of  the  Reduced  Plant  and  the  Reduced  Controller  (Overall 
System)- -The  interconnections  for  the  overall  system  are  shown  in  Figure  57b. 
The  non-null  interconnection  matrices  for  these  interconnections  are 
obtained  as  in  Section  V.  The  interconnection  data  are  used  by  subprogram 
STAMK3  of  KONPACT-1  to  model  the  overall  system.  The  model  following 
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Table  22.  Steady  State  Responses 


error  rates  for  handling  quality  are  defined  as  additional  responses.  The 
overall  system  is  given  a  system  number  6  (S-6,  Figure  57b). 


% 


w 

& 


Overall  System  (Design  Model)- -The  design  response  specifications,  sensor 
specifications,  control  input  specifications,  and  disturbance  input  specifica¬ 
tions  are  used  by  the  CONDK  subprogram  of  KONPACT-1  to  obtain  the  design 
model.  System  number  6  is  retained  for  the  design  model  (S-6,  Figure  57b). 
The  specifications  for  this  model  are  listed  in  Table  23. 

Controller  Design 

Figure  63  represents  a  block  diagram  of  the  design  process.  The  design 
response  weights  are  shown  in  Figure  64.  In  the  following  paragraphs, 
the  steps  for  obtaining  the  ALDCS  controller  design  are  described. 

Optimal  State  Feedback  Gains  for  ALDCS  Controller  Design--Using  the 
overall  system  design  model  described  earlier,  the  optimal  state  feedback 
gains  are  computed  by  the  DLAK  subprogram  of  KONPACT-2  for  the 
specified  quadratic  weights.  The  quadratic  weights  are  varied,and  the 
corresponding  optimal  state  feedback  gains  are  computed  until  the  ALDCS 
design  goals  are  met. 

Reduced  Optimal  Feedback  Gains — Starting  from  the  optimal  state  feedback 
control  law,  reduced  optimal  feedback  gains  are  obtained  by  the  FFOC 

0 

subprogram.  The  reduced  control  law  is  given  by  the  following  equations: 

U*a  =  (K1DBLA)  'a  +  (K1A21RL>  A21RL  +  (K1GLAF>  GLAF  <19) 

u*ei  =  (K2A21R)  A21R  +  (K2AFUS)  +  <K2TFUS)  TFUS  +  (K2P)  P-  {20) 
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White  Noise  Input  to  Pilot  Flite 


The  gain  variations  required  to  produce  the  reduced  optimal  feedback  gains 
are  shown  in  Figure  43  and  44  of  Reference  4. 

The  Tuning  of  Reduced  Optimal  Feedback  Gains  for  Steady  State  Requirements-- 
The  gain  K2^  is  adjusted  to  meet  the  steady  state  requirements  in  the  elevator 
loop  (Reference  4). 

Performance  Evaluation 

The  time  response  and  covariance  response  of  the  final  closed  loop  system 

are  obtained  using  the  DIAK  subprogram.  Figures  66  through  74  demonstrate 

that  the  ALDCS  load  relief  and  the  handling  quality  design  goals  are  satisfied. 

The  resulting  eigenvalues  for  the  free,  SAS  and  ALDCS  aircrafts  are  shown 

in  Figure  66.  Figures  67  through  69  present  the  angle  of  attack  (a) time 

responses  due  to  a  pilot  elevator  command.  For  the  free,  SAS  and  ALDCS 

aircrafts.  Figures  70  through  72  present  the  pitch  rate  (q)  time  response. 

Figures  73  and  74  present  a  tabular  listing  of  the  steady  state  values  of 

response  variables  due  to  a  1-g  incremental  maneuver  and  of  the  covariance 

results  due  to  T)  and  T|  , 
g  P 

DECK  SET-UP 

Figure  75  shows  the  simulator  deck  data  which  are  calculated  by  the  FLEXSTAB/ 
LSA  program.  Using  this  data,  the  plant  design  model  is  obtained.  The 
precompiler  data  needed  for  this  step  are  shown  in  Figure  76,  and  the  ' 
corresponding  KOMPACT-1  input  data  are  shown  in  Figure  77.  In  the  second 
step,  the  reduced  order  plants  (F24RR,  F24RT  and  F24TT)  are  obtained.  The 
precompiler  data  for  this  step  are  shown  in  Figure  78,  and  the  corresponding 

KONPACT-1  input  data  are  shown  in  Figure  79.  In  the  third  step,  the 

* 
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reduced  order  plant  is  combined  with  the  reduced  ALDCS  controller  to  obtain 
the  overall  system.  The  precompiler  data  and  KONPACT-1  input  data  tor 
this  step  are  shown  in  Figures  80  and  81.  In  the  fourth  step,  the  overall 
syrtem  design  model  is  obtained.  The  precompiler  data  and  KONPACT-1  input 
data  for  this  step  are  shown  in  Figures  82  and  83.  The  KONPACT-2  input 
data  required  to  obtain  optimal  state  feedback  gains  (using  the  DIAK  sub¬ 
program)  are  shown  in  Figure  84. 


The  KONPACT-2  input  data  required  to  compute  the  feedback  gains  only  on 
specified  measurements  (using  the  FFOC  subprogram)  are  shown  in  Figure  85. 
The  KONPACT-2  input  data  required  to  obtain  time  response  (using  the  DIAK 
subprogram)  of  the  closed  loop  system  are  shown  in  Figure  86.  The  KONPACT-2 
input  data  required  to  obtain  covariance  response  (using  the  DIAK  subprogram) 
of  the  closed  loop  system  are  shown  in  Figure  87.  The  KONPACT-2  input 
data  to  prepare  frequency  domain  data  for  the  LSA  program  are  shown  in 
Figure  88,  and  the  corresponding  data  for  the  LSA  program  to  evaluate  power 
spectral  density  are  shown  in  Figure  89. 

OUTPUT  DESCRIPTION 

KONPACT-1  output  data  are  shown  in  Figures  90  and  91.  The  C-5A  vehicle 
quadruple  data  along  with  the  name  list  data  are  given  in  Figure  90.  This 
system  1  quadruple  data  is  the  state  space  representation  of  the  unaugmented 
FLEXSTAB,  residual  elastic  math  model  of  Figure  75  and  represents  the 
output  resulting  from  the  KONPACT-1  input  data  deck  shown  in  Figure  77. 

The  system  ALDCS  overall  design  model  is  shown  in  Figure  91  and 
represents  the  output  from  the  KONPACT-1  input  data  deck  of  Figure  83; 
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the  model  contains  the  F24RR  model,  control  surface  actuator  dynamics, 
the  gust  model  dynamic^  and  the  ALDCS  controller.  This  model  is  used 
in  DIAK  and  FFOC  to  produce  the  reduced  controller  gains.  The  KONPACT-2 
output  data  are  summarized  more  fully  in  the  next  section. 

DISCUSSION  OF  RESULTS 

The  optimal  reduced  feedback  gains  for  the  ALDCS  design  are  shown  in 
Figure  65. 

The  eigenvalue  comparison  for  open  loop  vehicle,  SAS  vehicle,  and  closed 
loop  (reduced  feedback  only)  vehicle  are  shown  in  Figure  66.  The  time 
response  of  the  open  loop  vehicle,  SAS  vehicle,  and  closed  loop  (using  the 
reduced  controller  of  Figure  62)  vehicle  are  shown  in  Figures  67,  68  and 
69  for  response  cr  and  in  Figures  70,  71  and  72  for  response  q.  It  can  be 
seen  that  the  handling  quality  requirements  are  met. 

The  covariance  response  due  to  gust  load  for  open  loop  vehicle  and  closed 
loop  (reduced  feedback  only)  vehicle  are  shown  in  Figure  73.  It  is  seen  that 
the  gust  load  alleviation  requirements  are  met. 

The  steady  state  responses  for  1  g  maneuver  for  open  loop  vehicle  and 
closed  1oop  (reduced  feedback  only)  vehicle  are  shown  in  Figure  74.  It  is 
seen  that  maneuver  load  control  requirements  are  met. 

The  torsional  moment  (Tl)  design  requirement  for  this  controller  is  not 
satisfied.  Tl  is  increased  80  percent  due  to  this  reduced  controller  of 
Figure  62.  This  controller  was  designed  by  a  single  run  through  the  FFOC 
program  as  described  in  Reference  4.  Due  to  the  time  and  funding  contraints 
of  this  study,  the  controller  was  not  refined  to  meet  this  design  requirement. 
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Figure  6.  Control  Card  Arrangement  to  Create  VDATA  File 
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Figure  7.  Control  Card  Arrangement  to  Execute  KONPACT-1  and 
Create  QDATA  and  NDATA  Files 
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Figure  8.  Control  Card  Arrangement  to  Execute  KONPACT-1  and 
Add  New  System  Data  on  QDATA  and  NDATA  Files 
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Figure  11.  KONPACT-1  Modeling  Data  System 
Arrangement 
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Figure  12,  KONPACT-1  System  Reference  Control 
Card  Arrangement 
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Figure  13.  System  Description  Data  (for  using  STAMK1) 
Card  Arrangement 
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Figure  17.  System  Description  Data  (for  using  CONDK) 
Card  Arrangement 


Figure  17.  System  Description  Data  (for  using  CONDK) 
Card  Arrangement  (Concluded) 
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Figure  18.  Data  Written  on  Files  QDATA  and  NDATA 
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Figure  21.  Control  Card  Arrangement  to  Execute  KONPACT-2  Program 
and  Create  Files  DDATA,  FDATA  and  SDSTP 
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Figure  22.  Typical  KONPACT-2  Input  Data  Subprogram  Execution 
Arrangement 
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Figure  24.  KONPACT-2  Input  Data  (for  using  FFOC) 
Card  Arrangement  (See  Reference  2  for 
Variable  Definition) 
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Figure  25,  KONPACT-2  Input  Data  (To  Compute  Frequency  Domain 
Data  for  LSA  Program)  Card  Arrangement 
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Model  Generation  for  Design  of  Handling  Quality  Controller 
(C-5A  Cruise  Flight  Condition) 
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Figure  28.  Design  Response  Weights  for  Handling  Quality  Controller 
Design  (C-5A  Cruise  Flight  Condition) 
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Figure  29.  Variation  of  Design  Performance  with  Change  in  Quadratic  Weight 
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Figure  32.  Response  of  C-5A  Open  Loop  Static  Elastic 
System  to  Elevator  Command 
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Figure  32.  Response  of  C-5A  Open  Loop  Static 
Elastic  System  to  Elevator  Command 
(Concluded) 
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Figure  33.  Response  of  C-5A  Full  State  Feedback  Static  Elastic 
System  to  Elevator  Command 
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Figure  33.  Response  of  C-5A  Full  State  Feedback  Static  Elastic 
System  to  Elevator  Command  (Concluded) 
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Figure  34. 


Response  of  C-5A  Reduced  Feedback 
(Kfte  =  0)  to  Elevator  Command 
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Figure  34.  Response  of  C-5A  Reduced  Feedback  Static  Elastic  System 
(K*  *  0)  to  Elevator  Command  (Concluded) 
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T/RI 

4 

1 

0. 

0. 

0.  0. 

T/OELSO 

4 

2 

0. 

0. 

0.  0. 

0. 

3. 

0. 

0. 

T/DELS1 

4 

2 

.00770793 

0. 

0.  0. 

0. 

0. 

0. 

0. 

T/NGO 

4 

3 

0. 

0. 

0.  0. 

0. 

0. 

o« 

0. 

0.  0. 

6. 

0. 

T/WGSO 

4 

3 

00 

0. 

0.  0. 

0. 

0. 

0* 

0. 

0.  0. 

0. 

0. 

•FINISHED*  0  0 


Figure  35.  FLEXSTAB/LSA  Static  Elastic  Simulator  Deck  Data 
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stamki 

5 1 AHK£ 
STAHK3 
COMO* 
NXH-  6 
NUH*  lb 
NRN=  10 
N*M-  10 
MSft=  3 
MTR--  1 


Figure  36.  Figure  37  Precompiler  Data  (KONPACT-1) 


C  INPUT  OAT*  FOR  DEMONSTRATION  EXAMPLE 
C  SPECIFY  PRINTING 
PRINT  OUTPUT  DATA 
PRINT  INPUT  DATA 
C  DEFINE  VEHICLE 

SYSTEM  NO  I  VEHICLE  I  STATIC  ELASTIC  SYMMETRIC  > 

SLSA  DATA 

STATIC-ELASTIC  SYMMETRIC 

ENO 

C  NAME  LIST  DATA 
STATE 


1 

XI  II 

U 

VELOCITY  ALONG  X  AXIS 

2 

X<  £1 

« 

VELOCITY  ALONG  Z  AXIS 

INCH/SEC 

3 

XI  31 

0 

PITCH  RATE 

RADIAN/SEC 

A 

-1 

OUTPUT 

XI  A 1 

THETA 

PITCH  ATTITUDE 

RADIAN 

1 

R I  II 

SAS6Y 

PITCH  RATE  GYRO 

RADI AN/SEC 

2 

R <  2i 

AZAP 

NORMAL  accelerometer 

INCH/SEC2 

3 

R<  31 

azfb 

NORMAL  ACCELEROMETER  FRQNTSPAR 

INCH/SEC2 

A 

-1 

INPUT 

Rl  A I 

AZRB 

NORMAL  ACCELEROMETER  BACKSPaR 

INCH/ SEC?. 

t 

U(  |  1 

9DAIL 

AILERON  DEELECTION 

RADIAN 

2 

U<  2) 

BOELV 

ELEVATOR  DEFLECTION 

RADIAN 

w 

Ul  3) 

BOA ILOOT  AILERON  DEFLECTION  RATE 

RADIAN/SEC 

A 

Ul  A) 

RDELVDOT  ELEVATOR  DEFLECTION  RATE 

RADIAN/SEC 

s 

Ul  SI 

UGt 

GUST  INPUT  AT  -1020  IN  FROM 

CG 

INCH/SEC 

6 

Ul  6) 

*G2 

GUST  INPUT  AT  0  IN  FROM 

CG 

INCH/SEC 

7 

Ul  7) 

W03 

GUST  INPUT  AT  1020  IN  FROM 

CG 

INCH/SEC 

8 

Ul  B> 

WGIOOT 

GUST  INPUT  RATE 

INCH/SEC2 

9 

Ul  91 

WG200T 

GUST  INPUT  RATE 

INCH/SEC2 

10 

U 1 10  > 

WG300T 

GUST  INPUT  RATE 

1NCH/SEC2 

11 

U 1 1 1 ) 

WGS1 

STEADY  GUST  INPJT 

INCM/SEC 

12 

Ul  121 

WGS2 

STEADY  OUST  INPUT 

INCH/SEC 

13 

Ul  131 

KGS  3 

STEADY  GUST  INPUT 

INCH/SEC 

1A 

U  1  1 A  1 

WGS100T 

STEADY  GUST  INPUT  RATE 

INCH/SEC2 

IS 

UI15I 

NGS700T 

STEADY  GUST  INPUT  RATE 

INCH/SEC2 

16 

-1 

ENO 

U 1 161 

VGS3DOT 

STEADY  OUST  INPUT  RATE 

1NCH/SEC2 

Figure  37.  KONPACT-1  Input  Data  for  Static  Elastic  Model 
(Resulting  KONPACT-1  Output  Shown  in 
Figures  41  through  53) 
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C  DEFINE  REDUCED  VEHICLE 

SYSTEM  MO  l  VEHICLE  (  STATIC  ELASTIC  SYMMETRIC  -  REDUCED  > 

SCONO I T I  ON  I  NO  DATA 
C  NO  SCALING  DATA 
ENO 

C  RESPONSE  SpECIEICATIONs 
SEIECT  OUTPUTS 
»  ( I )  .R  <2»  .X  (2)  .X  (  3>  , 

ENO 

C  REDUCTION  ANO  SHUEELING  DATA 
RETAIN  STATES 
X  (?) .X (31 . 

RETAIN  INPUTS 
U<2>  . 

ENO 

C  DEE INE  ACTUATOR 

SYSTEM  MO  2  ACTUATOR 

s transfer  function  data 

BLOCK  1 

1  2  .7SOOOOE  01  2  1  . 1 OOOOOE  01  2  a  .750000E  01 
-1 
END 
UI/U 

1  I  • I OOOOoE  01 
-I 

R/R I 

1  1  • 1 OOOOOE  01 
-1 
ENO 

C  NAME  LIST  DATA 


STATE 

1 

-1 

OUTPUT 

I 

X  ( 

1  1 

DELE 

ELEVATOR 

OEFLECTION 

RADIAN 

R  ( 

1  1 

DELE 

elevator 

deflection 

RADIAN 

INPUT 

1 

Ut 

1* 

DELEC 

ELEVATOR 

COMMAND 

RADIAN 

-I 

END 


Figure  37.  KONPACT-1  Input  Data  for  Static  Elastic 
Model  (Continued) 
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i§gs 


C  0EF1NE  PILOT  MOOCL 
SYSTEM  MO  3  PILOT  MOOCL 

STRANSFER  FUNCTION  DATA 
BLOCK  I 

1  2  .223610C-03  2  1  .IOOOOOE  *122  .IOOOOOE  00 
-I 
ENO 
UI/U 

I  1  •IOOOOOE  #1 

-I 

B/BI 

I  1  .IOOOOOE  01 
-I 
END 

C  NAME  LIST  OATA 


3IMIL 

1 

-1 

OUTPUT 

1 

-1 

X< 

1  > 

XP 

PILOT 

MODEL  STATE 

RADIAN 

R  ( 

I  1 

FP 

PILOT 

COMMAND 

RADIAN 

INPUT 

1 

U< 

I  » 

CTAP 

PILOT 

MODEL  INPUT 

RADIAN 

-1 

ENO 

C  DEFINE  PLANT 

SYSTEM  NO  A  PLANT  I  PILOT  MODEL  *  ACTUATOR  ♦  VEHICLE  » 

INTERCONNECTION  DATA 

un/u 

-I 

UI2/U 

I  2  .IOOOOOE  01 
-I 

UI3/U 

1  1  .IOOOOOE  01 
-I 

UI1/RI2 

I  I  .IOOOOOE  01 
-I 

R/RII 

1  I  .IOOOOOE  01  2  2  .IOOOOOE  01  3  3  .IOOOOOE  01  A  A  .IOOOOOE  01 
-I 

R/RI3 

S  1  .IOOOOOE  01 
-1 
END 

C  NO  NAME  LIST  DATA  IS  NEEDED 
ENO 

C  DEFINE  IMPLICIT  MOOCL 

SYSTEM  NO  9  IMPLICIT  MODEL 

SOUADRUPLE  DATA 

XOOT/X  22 

1  1-.6T6600E  00  I  2  .87*120E  OA  2  1-.756200E-03  2  2-.352I30E  01 
-1 

XDOT/U  2  I 

1  1-. 330790C  03  2  1-.160637E  0! 

-1 

R/x  2  2 

1  1  .IOOOOOE  01  2  2  .IOOOOOE  01 

-I 
END 


Figure  37.  KONPACT-1  Input  Data  for  Static  Elastic  Model 
(Continued) 


C  NAME 
STATE 

LIST 

OATA 

V  , 

I 

X( 

n 

MI 

IMP 

MODEL 

VELOCITY 

FEET/SEC 

2 

-l 

OUTPUT 

X< 

2) 

01 

IMP 

MODEL 

PITCH  RATE 

RADIANS/SEC 

1 

R< 

It 

MI 

IMP 

MODEL 

VELOCITY 

FEET/SEC 

2 

-1 

INPUT 

R  ( 

2 1 

01 

IMP 

MODEL 

PITCH  RATE 

RADI ANS/SEC 

1 

U( 

L  • 

DELEI 

I  HP 

MOOEL 

INPUT 

RADIANS 

-I 

END 

C  OEriNE  OVERALL  SYSTEM 

SYSTEM  NO  5  OVERALL  SYSTEM  <  PLANT  ♦  IMPLICIT  MODEL  > 

SINTER CONNECT  ION  DATA 
UIA/U 

1  1  .IOOOOIE  01  2  2  . IOOOOOE  01 

UI9/U 

-l 

UI9/RI4 

1  5  .IOOOOOE  01 
-1 

R/RI A 

1  1  .IOOOOOE- 01  2  2  .IOOOOOE  01  3  3  .U31«0£*03  S  3  .IOOOOOE  01  6  A  .IOOOOOE  01 

4  S  .IOOOOOE  01 
-1 

R/RI9 

5  1*. IOOOOOE  01  6  2-. IOOOOOE  01 
-I 

ENO 

C  NAME  LIST  DATA 
OUTPUT 

3  R (  31  ALPHA  ANGLE  Or  ATTACK  RADIAN 

‘1 
ENO 

C  OEFINE  OVERALL  SYSTEM  WITH  DESIGN  RESPONSE  SPECIFICATIONS 
SYSTEM  NO  5  OVERALL  SYSTEM  I  DESIGN  MODEL  ) 

SCONOI T ION ING  DATA 
C  NO  SCALIN6  OAT A 
END 

C  RESPONSE  SPECIFICATIONS 
SELECT  CONTROL  INPUTS 
U<2I. 

SELECT  COMMAND  INPUTS 
U<!>. 

CONSTRUCT  DESIGN  RESPONSES 

X(2t  »RI3>  »X<3> .XOOTIJ>.ROOT<5»-ROOT<4). 

SELECT  SENSOR  OUTPUTS 
R(|I«R(?>»X(3)«R(A). 

ERO 

C  NO  REDUCTION  AND  SHUFFLING  DATA 

END 

STOP 


Figure  37.  KONPACT-1  Input  Data  for  Static  Elastic  Model 
(Concluded) 
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C  DESIGN  USING  OIAK  FOR  THE  DEMONSTRATION  EXAMPLE 
READ  FOR  WHAT  PROGRAM  (  DIRK  «FF  OC  »l.S*  t  THE  DATA  IS 
40IAK  nATA 

C  READ  IF  0*TA  IS  ON  CAROS  ONLY  OR  ON  CAROS  AND  TAPE 
DATA  ON  CAROS  AND  TAPE 

C  IF  DATA  is  ON  CARDS  AND  TAPE  READ  THE  LABEL  TO  OBTAIN  OATA  ON  TAPE 
STSTEH  NO  5  OVERALL  SYSTEM  '  DESIGN  MODEL  > 

C  READ  DATE  AND  USER  ID 
AUG  17.  7S  J  K  HAMESH 
C  NOP  -  NO  OF  VARIABLES  being  PLOTTED 
C  REAO  NOP 
0 

C  GO  TO  100  IF  NOP.ER.O 

C  READ  (RLR ( I) • ITITl ( I) .UNIT ( I> .YMINI I > .VMAXt I> .SCAL  <I>  *l»I.N0P> 

C  READ  T.DT.ST.T1.T2 
C  100  CONTINUE 
C  READ  IMA*. ITER. 1TCRQ 
4030  A 

C  N0C0V=1  NO  COVARIANCE  ANALYSIS 
C  NOCOV*2  COVARIANCE  ANALYSIS 
C  NOCOVO  SNIP  CORRELATION  ANALYSIS 
C  NSTEP=0  NO  STEP  INPUTS 
C  NSTEPs 1  STEP  COMMANDS 
C  NSTEP=2  STEP  GUSTS 
C  NSTEPO  BOTH  (l  AND  2) 

C  NSTEP=M  NO  STEP  INPUTS  -  TRANSIENTS  ONLY 
C  NRAND*0  NO  RANDOM  INPUTS 
C  NpANO* I  GUSTS 

C  NPRIN=P  DO  NOT  PRINT  RESPONSES 
C.  NPR IN*  I  PRINT  RESPONSES 
C  NPLOT *0  NO  PLOTS 
C  NPLDT=1  C«LCOMP  plots 
C  NPLOT *2  LINE  PRINTER  PLOTS 
C  NPLOT* 3  POTH  (1  AND  ?) 

C  READ  NOCOV.NSTEP.NRANO.NPR1N. NPLOT 

3  0  0  0  0 

C  INDK*I  NEW  INPUT  GAINS 
C  INP<*?  NEW  STARTING  ROUTINE  GAINS 
C  INP«*3  USE  GAINS  IN  STORAGE 
C  INPM*A  USE  INPUT  GAINS  IN  STORAGE 
C  READ  INPK 
I 

C  NCONT *0  OONOT  COMPUTE  OPTIMAL  GAINS  -  USE  INPUT  GAINS  ANO  DATA  IN 
C  COVARIANCE  ANO  TIME  RESPONSE  ANALYSIS  ONLY 

C  NCONT  *  1  COMPUTE  OPTIMAL  GAINS 

C  NC0NT*2  COMPUTE  OPTIMAL  GAINS  WITH  AUTOMATIC  0  SELECTION  ON  CONTROL  RATES 
C  READ  NCONl 

1 

C  READ  FLIGHT  CONDITION  NUMBER 
45 

C  N*  -  NO  OF  STATES 
C  NR  -  NO  OF  RESPONSES 
C  NU  -  NO  OF  CONTROL  INPUTS 
C  NN  -  NO  OF  DISTURBANCE  INPUTS 
C  NF  -  NO  OF  FEEDBACK  STATES 
C  NG  -  NO  Of  GUST  INPUTS 

C  NCS  -  NO  OF  COMMAND  INPUTS  *  NO  OF  COMMAND  STATES 
C  NGLG  -  NO  OF  GUST  LIFT  GROWTH  STATES 

C  NSCRR  -  START  OF  CONTROL  RATE  RESPONSE  IN  THE  RESPONSE  VECTOR 

C  READ  NX.Nft.NU.NN.NF.NG.NCS'.MGLG.NfCRR 

4  6  114  0  0  0  7 


Figure  38.  KONPACT-2  Input  Data  for  Static  Elastic  Model 

(Employing  DIAK  to  Compute  Optimal  State  Feedbc  , 
Gains)  (Resulting  KONPACT-2  Output  Shown  in 
Figure  54) 
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) 


L  GO  TO  200  IF  INPK.GT.l 
C  RE AO  <NORQ<  n . 1  =  1. NX  > 

12  3a 

C  200  CONTINUE 

C  F  IS  STATE  TRANSITION  MATRIX 
READ  TAPE  FOR  MATRIX  F 
C  Gl  IS  CONTROL  INPUT  MATRIX 
REAO  TAPE  FOR  MATRIX  Gl 
C  G2  IS  DISTURBANCE  INPUT  MATRIX 
REAO  TAPE  FOR  MATRIX  02 
C  XI  IS  INITIAL  CONDITION  MATRIX 
READ  CARO  FOR  MATRIX  XI 

C  XLOXL  IS  STATE  LIMIT  -  RATE  LIMIT  MATRIX 
REAO  CARD  FOR  MATRIX  XLDXL 

1  I  . 1 00090E  20  2  I  . 100000E  20  3  1  .IOOOOOF  20  4  1  .100000E  20 

1  2  .IOOOOOF  20  2  2  .IGOOOOE  20  3  2  .100000E  20  4  2  .100000E  20 

C  CL  IS  command  level  matrix 
READ  CARD  FOR  MATRIX  CL 

C  M  IS  STATF  RESPONSE  MATRIX 
REAO  TAPE  FOL  MATRIX  H 
C  0  IS  CONTROL  RESPONSE  MATRIX 
REAO  TAPE  FOR  MATRIX  D 
C  AM  IS  MEASUREMENT  MATRIX 
READ  TAPE  FOR  MATRIX  AM 
C  RK  IS  INITIAL  FEEDBACK  GAIN  MATRIX 
READ  CARO  FOR  MATRIX  RK 

C  0  IS  QUADRATIC  WEIGHTS  MATRIX 
READ  CARO  FOR  MATRIX  Q 

3  3  .IOOOOqE  01  4  4  . 1 OOOOOE  01  6  6  .lOOOOOE-Ol 

C  I0UM«0  ANOTHER  run 
C  IDUM*1  NO  MORE  RUNS 
C  READ  IDUM 
0 

C  INpO*l  COMPLETELY  new  data 

C  IHP0*2  CHANGE  SELECTED  QUADRATIC  WEIGHTS  ONLY  -  USE  SOME  GAINS  IN  STORAGE 
C  INPO*3  CHANGE  SELECTED  QUADRATIC  WEIGHTS  ONLY  WITH  OPTION  FOR  NEW  GAINS 
C  INP0*4  CHANGE  SELECTED  data 

C  INPD*5  CHANGE  SELECTED  DATA  IN  MEASUREMENT  MATRIX  AND  QUADRATIC  WEIGHTS 
C  WITH  OPTION  FOR  NEW  GAINS 
C  READ  INPD.INPK 

2  3 

C  READ  NCONT 

1 

C  READ  t!OCOv.NSTEP«NRANO.NPRIN*NPLOT 

3  0  0  0  0 

REAO  CARD  FOP  MATRIX  0 

3  3  .10000  IF  01  4  4  .lUOOOOE  01  <S  4  .100000E-00 

C  READ  IDUM 
0 

C  READ  INPD.INPK 

2  3 

C  REAO  NCONT 

C  READ  NOCOV.NSTFP.NRAND.NPRIN.NPLOT 

3  0  0  0  0 

READ  CARO  FOP  M4TPJX  0 

3  3  .10000»£  01  4  4  .100000E  01  4  4  .IOOOOOE  01 

Figure  38.  KONPACT-2  Input  Data  for  Static  Elastic  Model  (Employing 
DLAK  to  Compute  Optimal  State  Feedback  Gains)  (Continued) 


C  BEAD  I OUM 
0 

C  READ  INpO.JNPK 

?  3 

C  READ  NCONT 
I 

C  READ  NOCOV.NSTEP.NRANO.NPRIN.NPlOT 
3  0  0  0  0 

READ  CARD  EO»  MATRIX  « 

3  3  . 1 OOOOOE  01  4  4  .100000E  01  &  0  .IOOOOOE  OJ! 

c  head  idum 

0 

C  BEAD  INPD.1NPK 

i  3 

C  BEAD  HfONT 
l 

C  BEAD  NOCOV.NSTePtNRAR0.NPRJN*NPLOT 
3  0  0  0  0 

BEAD  CARD  E^B  MATRIX  0 

3  3  .IOOOOOE  (/I  4  4  .IOOOOOE  01  6  6  .100000E  03 

C  PEAD  IDUM 

1 

END  OF  DIAK  DATA 
STOP 


Figure  38.  KONPACT-2  Input  Data  For  Static  Elastic  Model  (Employing 
DIAK  to  Compute  Optimal  State  Feedback  Gains)  (Concluded) 
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c  REAO  FOR  -AT  PROGRAM  I  01AK.rnC.L5A  )  T ME  0»TA  15 

cFprAr)°irAf.;T4  i~  r.<  cab..s  only  o-  on  laus  and  tape 

C*IF  r>AT  A^*!  '  ^ON^C  AflFI .  A  NT  T  APr  wf.AD  Tut  LAm£l  TO  OBTAIN  DATA  ON  TAPE 
SYSTEM  Ni  S  OVERALL  SYSTEM  <  DESIGN  MODEL  ) 

C  INAa  *  wA «  I  Miff  NO  OF  •  vPUNOV  SOlUI.jN  ITERATIONS 
C  N I T v  -  MAflMH;*  NO  OF  COST  CALCULATIONS 
C  NOPR=D  JS-  PROJECTED  GRADIENT 
C  NORSs I  DO-.oT  USE  PROJECTED  gBrDTEnT 
C  NOCOV^l  Ni  COVARIANCE  ANALYSIS 
c  NOCOVsP  C')NA(T  I  ANCt  ANALYSIS 

i  sssc.w'.grsa'jjw’tsj? « »«*.. 

C  NBEOINsO  ;0  test 

C  PEAlO  I«A  A  .NI  TM.NOPB. NOCOW. NMEGiN 

in  n  i 

C  NX  •  NO  OF  STATES 
C  NS  -  NO  OF  SF5P0NSFS 
C  NJ  -  NO  OF  CONTROL  INPUTS 
C  NN  -  SO  OF  DISTURBANCE  INPUTS 
C  NFF  -  NO  rFEO  FORWARD  STATES 
c  NF  -  NO  OF  FIXEO  GAINS 
C  READ  NX.N /.NU.NN.NFF.NF 
4  G  I  1  1  1 

c  BEAD  (NORmI I ) . I*t *NX> 

1  ?  0  A 

C  E’ST  •  INITIAL  STEP  S?7E 
C  REAP  E°SI 
• IOOOF  00 

C  AJSTAP  -  :.OWEbT  COAT  E»uErTED 
C  RE  AO  AJSTAB 
.9993E-0R 

C  CROC  -  DE‘IBED  RATIO  OF  COSTS 
C  READ  0«OC 
.HOOF  01 

C  A|_A>*  -  INTEGRATION  PARAMETER  -  LAmOA 
C  OELT  -  INTEGRATION  STFP  size 

C  AlAMO  -  L  'WF.m  ROUND  ON  LANDA  FOS  the  present  run 
c  RFAO  AL AM.O£L T . ALA“0 
.1000E  o l  .PJC.E  00  .OOSLE  00 
C  IF  -  fixed  GAIN  PON  INDEX 
c  Jf  -  FIXED  GAIN  COLNMN  INDEX 
C  RE  AO  (  IF  (  •  I  .  JF  (  I)  .  1=1  .'.'FI 


1  1 

l  2  1 

read 

T  A3E 

r  \q 

matrix 

F 

READ 

T  A°E 

f^r 

MATRIX 

Gl 

RF.A0 

ta3E 

T'iP 

MATRIX 

02 

READ 

TA=E 

r  >r 

matrix 

H 

READ 

T  A°E 

F'lP 

MATRIX 

f) 

READ 

TAOE 

FOR 

MATRIX 

AM 

9£AD 

CARO 

FOR 

MATRIX 

0 

3  3 

.1000)  E 

91  4  4 

.lOriMF  01 

G  b  .  1 

ftFAO 

C»RO 

F  »R 

MATRIX 

ak<wopt  iv*l 

SICCA  1 

1  l 

.5050S-F-LJ  I  P 

,2bS790E  H 

l  t-.: 

READ 

CARO 

F'»R 

MATRIX 

AK 1  (  I  )  ) 

READ 

CARO 

FOR 

MAT©! * 

RK (K?) 

C  If 

(ALAM#r,T 

..99)  GO  TO  100 

RE  AO 

CARO 

f  tR 

MATRIX 

AK  < K 1 (LAMBDA) ) 

REAO 

CARD 

F^R 

MATRIX 

DELK  <DFL*1 (LAMBDA 

C  100  CONTINUE 
FNO 


Figure  39. 


KON PACT-2  Input  Data  for  Static  Elastic-  Model  (Employing 
FFOC  to  Compute  Reduced  Feedback  Gains)  (Resulting 
KON  PACT -2  Output  Shown  in  Figure  55) 
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C  DESIGN  USING  i»a*  f  r-  .1,-  ..  1  ..  (.  jTijvi  f  m«"-. 

C  REAP  f(W  -AT  P-OS-’.*-'  I  ,L.-.4  )  '-It  .lit.  1 -• 

SOIAK  OATA 

C  °EAO  IF  OTA  IS  O')  <■.->••  s  I'll  V  1)3  CA-  |>  tvo  r apl 
DATA  ON  CAPE'S  ANO  Ta^t 

C  If  OATA  It  ON  f»»0'  A'lr*  TARf  -a‘l'  In-  L»H-L  I'  OOIaIN  IiATA  ON  TAPE 
SYSTEM  NO  s  ■jVti'GLl.  Y  S  T  L  .  1  uCSluS  MOll'L  ) 

C  PE  AO  OATS  AN1)  USE-  f 
AUG  17.  7S  J  *  <a«--s- 
C  NOP  -  NO  "  VAO;A-|l.fS  i.z  l  ■  r.  D.rtrT-l. 
c  orAO  NOP 
3 

C  r.o  TO  100  If  voo.fO.  • 

C  READ  IPLPi  II  .  ITiTi.  1 1 1  .111!  r  I  I  •  1  U  •(  I  I  .  1  Mt.  1 1 1  .SCAL'  l)  .Isl.'iOPi 
1  0 

2  ALPHA 

3  DEI.TAE 

C  READ  T.nT.ST.Tl  .T’ 

S.G  *).  'P  O.A  w  •  '*  ‘1.0 

C  100  CONTI  It  If 
C  PEAO  IMA«.ITf».IT.-:-j 
*030  4 

C  NOCOVs]  NO  COVA->!4>rf  ANALYSIS 
C  N3C0V*?  C'VAPIA-lCf  ANALYSIS 
C  NOCOVfl  S-IP  COPPCLATIVi  ANALYSTS 
c  NSTfPzO  N  .  STEP  I>|P||T>> 

C  NSTEP=1  STfP  COMMANDS 
C  NSTE°*P  STfP  GUSTS 
C  NSTEPfO  e-'Tn  (1  A.NO  ?i 

c  nste°*4  in  step  inputs  -  tk ans i rN r s  only 

C  N»4N0*0  NO  PA.-IOOM  INP.jYS 

c  naano*  1  gists 

C  NPRINfO  0-1  not  PRINT  RESPONSES 
C  NOQIN*l  PSTNT  RESPONSES 
C  N°L0T *0  NO  PLOTS 

C  n»lot*i  c *lco *°  plots 

C  N»LOTi»  LTNf  PPl'ir^u  PI.  OT; 

C  NPLOTfT  NOTH  (|  A'H  P> 

C  READ  NOCCi'.NSTEO.Npania. -iPJIN.nPlO' 

I  1  "  I  ? 

C  INPS  I  NF  ■*  INPUT  GAINS 
C  INPK*?  NF-  STARTING  SOUTINE  SAINS 
C  1 NPS * 3  USf  SAINS  1. 1  STOPAG-- 
C  INP<*4  USr  INPUT  GAINS  IN  STOJA'.t 
C  READ  IMPS 
I 

C  NC0NT*0  O'.NOT  COhPijTE  optimal  GaInS  -  USE  INPUT  GAINS  AND  DATA  IN 
C  COVARIANCE  AND  tint  RESPONSE  ANALTSIS  only 

c  ncont *  1  r-pouif  optima/  gains 

C  NCONT.A  compijt £  opt  tpai  gains  *iT-  AUT jAAl  !C  o  SFLfCTIOIl  3-v  CONTROL  PAltS 
C  PEAO  NCON* 
t* 

C  PEAO  FLIGHT  CONDITION  '■i|M«ca 

4* 

C  NX  -  NO  OF  STATES 
C  NO  -  VO  Of  RESPONSES 
C  NJ  -  NO  Of  CONTROL  »N»-iTS 
C-NN  -  NO  Of  DISTURBANCE  INPUTS 
C  NF  -  NO  Of  fFEORACY  Sfe'ES 
C  NG  -  NO  Of  GUST  INpuTr 

r.  DCS  -  .'If.  Of  c*.  ,r  -•  i>  :•  :C--«N>  STA'ES 

C  NGLG  -  *'0  Of  I,i *S T  I  ..  M.;.  -'1*.  ES 

C  NSCP"  -  ST'-T  Of  *T M I  L  -  If  P-'S-'>.iSE  IN  T-t  PESPiiNSF  VF.fTOP 
c  PFAO  •IX.NR.  l'I.M  ••IF.I.I..  .  .SLG.N'C  IP 

4*1100'  ’ 

C  go  to  p'f  :»  gi  .  ‘ 

r  Of  AO  I'.OPpl  1 » •  1  r\  .1 
I  >  3  > 

r  *s*»  co’itiwe 

Figure  40.  KONPACT-2  Input  Data  for  Static  Elastic  Model 
(Employing  DIAK  to  Evaluate  Time  Responses) 
(Resulting  KONPACT-2  Output  Shown  in  Figure  56) 
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Figure  40 


r»5r  :  • .  . 

r  -  r  .  ,t  .  , 

rra->  TA=*r  *•  •  -'■')•  •' 

r  r. ■>  T ^  f  •  T  t  •  *  ’  .  *  •  ■- 

cr  an  1  *  ••.!-** 

r  *  T  T4'  »*<»  l~.  '■  V"*'*  : 
ocan  r  a  ?•  r  *•::-!»  s  * 

r  »L«*l  ?  :  *T-.T'  •_  I  '*  -  •  .'  • 

erAT  c.r.-^  c  .J  ••  7.  .  n_  .>!_ 

I  »  #  t  >  •»  1  .  ] 

1  *.!*•"•  *  *  *  5  . 1 

r  f'L  TP  r*v«  V- ■)  I.  ■■••-. L  “  T-f' 

orao  r*or>  F  z  •  ■  •  t  ,  ]  <  .  , 

I  1  ,4W’KV  - 

r  M  TS  STAT  •  ,  ..  J.  4..  •, 

or  AT  T-or  e  ■ 

f  f>  TC  *'i':s-T  v.  J  “«.»-■’!  •  r - '  T  •. 

or  An  TA“*r  r  ;  ••  •.  7-  • 

r  am  15  '  »•-*"  *  r  r.  .  . 

prftn  T*or  ►  t-  .. 

r  °<  T\!*!AL  »  *  •V»'i  :MT.'|* 

opao  r  ao:>  r  >  ••  .1-  • . 

l  1  ,111^  r  I  ,  /  ..  -  .  i-.-n.i**iur  -.10  )  4  .l-*ov*0f  01 

r  0  ?«.  wv.  *.!  * :  -*ff ••-*?'*  c:  •;* 

orao  caon  r  ..  r 

C  ?  •.*  *•'  *  *  4 

r  rnijMs;  \r  ”V~ 

r  °Fto  1 

• 

TMn  0r  ^  f  /*  •*  •  •  A  T  a 

STHO 


L  »  •  r  I  •**  A  r  •'  i  A 

•  co  4  1  , 1  (jooooh  ^o 

1  *  *■  #U<‘J0Or  to  ^  t  . lOOOOCt  20 


KONPACT-2  input  Data  for  Static  Elastic  Model  (Employing 
D1AK  to  Evaluate  Time  Responses)  (Concluded) 


•KMCMMMttMM 


*  SYSTEM  NO  1  VEHICLE  (  STATIC  ELASTIC  SYMMETRIC  » 


LSA  -  FLEXSTAB  data  ••• 
STATIC-ELASTIC  SYMMETRIC 


VP/VPO  MATRIX  .  SIZE  *  3  X  3 


1 -COLUMN 


2-COLUMN  3-COLUMN 


1  -ROM 

2- ROW 

3- ROW 


-.3979210E-02 
-.1189251E*00 
. 1 920000E-05 


-.2045144E-01 

-.678679c:*00 

IB74000E-03 


-•33121 15E*03 
•874|202E»04 
-•1101 102E*0l 


VP/RO  MATRIX  •  SIZE  >3X| 


1 -COLUMN 


1- ROW  ».3858256E»C3 

2- ROW  14851B5E*02 

3- ROW  .4919200E-03 


VP/OELSO  MATRIX  «  SIZE  >3X2 
1 -COLUMN  Z-COLUMN 


1- ROW  -. 1 162480E-01  -.3486477E*01 

2- ROW  - , 1544287£»03  - . 3307943E»03 

3- ROW  -,2908141E*00  -.  160636<»E*01 


I -ROW 


IBANOING)  MATRIX  •  SIZE  > 
1 -COLUMN  2-COLUMJ4 

-. 141 7935E«04  0. 


3-COLUMN 

•101793SE*04 


Figure  41.  KONPACT-1  Output  for  C-5A  Static  Elastic  Model 

(Cruise  Flight  Condition)  (See  Figure  37  for  KONPACT-1 
Input) 
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mtf\ 


VP/wGO 


MATRIX  •  SIZE 


3X3 


I -COLUMN 


2-COLUMN 


3-COLUMN 


.  1 1 34790E-0? 
.  1#07768E*00 
.4263200E-03 


•l^ZbRSAE-Ol 
•578021 3E*00 
•19O2000E-03 


.471 0000E-04 
■  •  1 203000E-03 
■•3970000E-04 


VP/WG1 


MATRIX  .  SIZE  =  3  X 


1 -COLUMN 


2-COLUMn 


3-COLUMN 


■ .  6400000E-04  -  •  77 1  OOOflE-03 

■. I 746490E-02  -.1756977E-01 

■.2000Q00E-OS  . 3970000E-04 


•6610000E-OS 
•12A6630E-0? 
■  .4030000E-05 


VP/WGSO  MATRIX  .  SIZE  *  3  X  3 


1 -COLUMN 


2-COLUMN 


3-COLUMN 


•1I34790E-02  .19Z6954E-01  .47JOOOOE-04 
,1007788E»00  .S78021 3E*00  -. 1203000E-03 
.4263200E-03  -. 1992000E-03  -.3970000E-04 


VP/MGS  1  MATRIX  .  SIZE  *  3  X  3 


I -COLUMN 


2-COL UNN 


3-COLUMN 


1  -ROM 

2- ROW 

3- ROW 


•  . 640000  *t-04  - . 77 1 OOOOE-03  .66I0000E-05 
-.1746490E-02  -.1756977E-01  .1246M0E-0? 
- . 2000000E-05  .3970003E-04  -.4030000E-05 


R/VPO  MATRIX  .  SIZE  -  1  X  3 


1 -COLUMN 


2-COLUMN 


3-COLUMN 


1-ROW  -.0 


.IQOOOOOE’Ol 


Figure  41.  KONPACT-1  Output  for  C-5A  Static  Elastic  Model 

(Cruise  Flight  Condition)  (See  Figure  37  for  KONPACT-1 
Input)  (Continued) 


[ 


R/RO  MATRIX  .  SIZE  3  I  X  I 


1 -COLUMN 


l-ROW  - ,0 


T/VPO  MATRIX  .  SIZE  *  A  X  .1 


1 -COLUMN 


2-COLUMN 


1- ROW  0. 

2- ROW  0. 

3- ROW  0. 

4- ROW  0. 


3-COLUMN 


•  I OOOOOOE  *01 
■,8839912E*04 
•>8B39912E*04 
>.BB3991ZE*04 


MATRIX  •  SIZE  *  4  X  3 


I -COLUMN 


2-COLUMN 


3-COLUMN 


1- ROW  . 3670000E-06  .4320000005  .3448910E-0? 

2- ROW  0.  . lOOOOOOC’O 1  -.7242167E»03 

3- ROW  0.  • 1 OOOOOOOO 1  .24A8J33E»03 

4- ROW  0.  >1000000001  >3498233E*03 


MATRIX  *  SIZE  »  4  X  1 


(-COLUMN 


1- ROW  0. 

2- ROW  ,14S8824E»02 

3- ROW  .1458824E*02 

4- ROW  ,14S8824E*02 


MATRIX  *  SIZE  «  4  X  1 


l-COLUMN 


1- ROW  0. 

2- ROW  0. 

3- ROW  0. 

4- ROW  0. 


Figure  41.  KONPACT-1  Output  for  C-5A  Static  Elastic  Model 

(Cruise  Flight  Condition)  (See  Figure  37  for  KONPACT-1 
Input)  (Continued) 


3^283* 


\  m 


T'"W ■  — 


T /DO- SO  MATRIX  .  SIZE 


AX  2 


1-COLUMN  2-COLUMN 


1- ROW  0.  0. 

2- ROW  0.  0. 

3- ROW  9.  0. 

4- ROW  0.  0. 


T/DELSt  MATRIX  i  SIZE  =4X2 


1 -COLUMN  2-COLUMN 

1- ROW  .7707R39E-02  0. 

2- ROW  0.  0. 

3- ROW  0.  0. 

4- ROW  9,  0. 


T/*G0  MATRIX  •  S17E  =4X3 


1  -ROW 

2- ROW 

3- ROW 

4- ROW 


I -COLUMN 

2-C0LUMN  3-COLUMN 

/ 

f 

0. 

0.  /•• 

0. 

0.  /  0. 

0. 

0 .  y  ‘  0 . 

0. 

0.  y  0. 

T/H6S0  jHtfRlX  •  SIZE  =4X3 


\ -COLUMN  2-COLUMN 


3-COLUMN 


1- ROW  0.  0. 

2- ROW  0.  0. 

3- ROW  0.  0. 

4- ROW  0.  0. 


0. 

0. 

0. 

0. 


•FINISHED* 


Figure  41.  KONPACT-1  Output  for  C-5A  Static  Elastic  Model 

(Cruise  Flight  Condition)  (See  Figure  37  for  KONPACT-1 
Input)  (Concluded) 


SYSTEM  NO  1 


VEHICLE  (  STATIC  ELASTIC  SYMMETRIC  > 


NUMBER  OF  STATES  *  4 
NUMBER  OF  OUTPUTS*  4 

number  of  inputs  *16 


*•*  NAME  list  table 


VARIABLE  NANE 


description 


u  VELOCITY  ALONG  a  AMS 

W  VELOCITY  ALONG  Z  AMS 

G  PITCH  RATE 

THETA  PITCH  ATTITUDE 


INCH/SEC 

INCM/SEC 

RAOIAN/SEC 

RADIAN 


OUTPUT 


sasgy  pitch  rate  gyro 

AZAP  normal  accelerometer 

azfb  normal  accelerometer  FROnTSPAR 

azrb  normal  accelerometer  backspar 


RAOIAN/SEC 

IMCH/SEC2 

INCH/SEC2 

INCH/SEC2 


BOAIL  AILERON  oeflection 
roelv  elevator  deflection 
boailoot  aileron  oeflection  rate 
roelvdot  elevator  deflection  rate 
NGl  GUST  INPUT  AT  -1020  IN  f ROM  CG 

*f,2  GUST  INPUT  AT  0  IN  FROM  CG 

Nf,3  GUST  INPUT  AT  1020  IN  FROM  CG 

MGlOOT  GUST  INPUT  RATE 

WG2D0T  GUST  INPUT  RATE 

VG300T  GUST  INPUT  RATE 

WOS)  STEADY  GU5T  INPUT 
MGS2  STEADY  GUST  INPUT 
VGS3  STEADY  GUST  INPUT 
MG51O0T  STEADY  RUST  INPUT  RATE 
VGS200T  STEADY  GUST  INPUT  RATE 
VGSIl'OT  STEAOY  GUST  INPUT  RATE 


RADIAN 

RADIAN 

RADIAN/SEC 

RAOIAN/SEC 

INCH/SEC 

INCH/SEC 

INCH/SEC 

INCH/SEC2 

INCH/SEC2 

INCM/SEC2 

INCH/SEC 

INCH/SEC 

INCH/SEC 

1NCH/SEC2 

INCH/SEC2 

INCH/SEC2 


•••  QUADRUPLE  DATA 


Figure  42.  KONPACT-1  Output--C-5A  Static  Elastic  Vehicle  Name  List 
Table  and  Quadruple  Data  (Cruise  Flight 'Condition)  (See 
Figure  37  for  KONPACT-1  Input) 


I 


/ 


MATRIX 

A 

SIZE  «  <• 

X  4 

1 -COLUMN 

2-COLUMN 

3-COLUMN 

4-COLUMN 

1- ROW 

2- ROW 

3- ROW 

4- ROW 

-•3979?|*C-»? 
-.11 99251000 
. 1920000005 
0. 

-.2045144E-01 

-.6786798E*00 

-.1874008E-03 

0. 

-.3312115E-83 
. 8  74 1 202E  *04 
-.1101102E»0I 
•  1806800E-8 1 

- • 3858256E  *03 
-.1485185E-82 
.49 1 9288E-83 
0. 

1 

1 

MATRIX  0 

SIZE  *  4  X 

16 

1 -COLUMN 

2-COLUMN 

3-COLUMN 

4-COLUMN 

5-COLUMN 

6-COLUMn 

7-COLUMN 

1- ROW 

2- ROW 

3- ROW 

4- ROW 

-.11  6248000 1 
-  .  1 544287003 
-.?908i4iooo 
0. 

-.34B6477E*01 
-.3307943E-03 
- • I60636QE  *0 1 
0. 

0. 

». 

«. 

0. 

0. 

«. 

«. 

0. 

.1I34790E-02 
* 1007788E-00 
•42632Q8E-03 
8. 

•  1926954E-01 
•578821 3E-88 
-.  1992000E-03 
8. 

•4710000E-04 

-.1203000E-03 

-.3970000E-04 

0. 

8 

i 

a 

o- COLUMN 

9-COLUMn 

io-column 

11 -COLUMN 

12-COLUMN 

13-COLUMN 

14-COLUMN 

1 

1- ROW 

2- ROw 

3- ROW 

4- ROW 

-.640000000* 

174649000? 

o. 

-•7710000E-03 

-.I7S6977E-01 

.3970880E-84 

0. 

•6610U80E-85 
.  I246638E-82 
-.4030000E-05 
0. 

•  1 1 34790E-82 
. I807788E  *00 
.4263200E-03 
0. 

.1926954E-01 
.5780213E-00 
-, 1992000E-03 
0. 

.4710000E-04 
- • 1 203000E-03 
-.3970000E-04 
8. 

-.6400600E-04 

-.1746490E-02 

-•2000000E-05 

0. 

1 

I 

Is-COLUMN 

16-COLUMN 

1- ROW 

2- ROw 

3- ROw 

4 - ROW 

- .  7  71 0000003 
-.1  756977001 
.  3970000004 
0. 

•66I0068E-05 
.1246630E-02 
-  .4030000E-OS 
«. 

MATRIX  C 

SIZE  ■  4  X 

4 

1 -COLUMN 

2-COLUMN 

3-COLUMN 

4-COLUMN 

1- ROW 

2- ROw 

3- ROW 

4- RQW 

-.5085949006 
- ,  1 2031 56000 
-.1 1R4$I?E*00 
-.11 825 34£ • 00 

-.358S728E-85 

-.S429616E-88 

-•7249364E*00 

-.?442367E*60 

•  1 O33043E  *0 i 
.69S7265E'83 
-.3704996E*#! 
-.4839816E»83 

-.2I48614E-83 
-. 619868 7E-80 
-.142181 8E *00 
-.915I893E-6I 

Figure  42.  KONPACT-1  Output--C-5A  Static  Elastic  Vehicle  Name  List 
Table  and  Quadruple  Data  (Cruise  Flight  Condition)  (See 
Figure  37  for  KONPACT-1  Input)  (Continued) 
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MTRIX  0 


1  -ROW 

2- BOW 

3- BOW 


I -COLUMN 


size  >  u  i6 

2-COLUMN  3-column  4-COLUMN 


-.|670SSAe-02  -.6970532E-0?  .7707930E-02 

.561 8  375E  •  02  .8325667E *03  0. 

- .22621 1 3E  *03  -.727299SE *03  0. 


*-eow  -,2SM62?C*03  -.»9?71«J4E«03  0. 


1- l>Ow  -.  I«.<.*>6l4f-fl7 

2- 90M  -,?98CHe>6C-03 

3- ftOW  ??•**  1S7E-92 

4- PO*  37E-0* 


I  = -COLUMN 


l  -ROW  .M>73736F-©7 

2- ROW  -.4632U?E-0l 

3- iitOW  -./  770‘.R9E-©» 

*-R0W  -.36H17A6E-0? 


9. column 


,f>07373<*F-07 
-•«6jru  7E-oi 
-  ,  7  7T')4HQE-02 
-.36«l78*E-02 


-.«SU?fc*E-08 
.M6S223E-0? 
.2S18O19E-03 
-  «  I  b3 IS7ME-03 


10-COLUMN 


0. 

0. 

0. 


0. 

II -COLUMN 


5-COLUMN 


•1906I20E-05 
-•2079693E*00 
•  296OO80E  *00 


6-COLUMN 


•  18 1 7 1 0  IE-05 
•7222S5lE*00 

•  526052  IE  *00 


7-COLUMN 


-  •  1 37424 lE-06 
•20631 10E-01 
-•9919581E-02 


-.85I1240E-0M  .1906I20E-05 

.<*l652?3E-0?  -,20796«J3E*00 

.2Sl«919E-0t  .2060088E*Q0 

-.I63IS78F-03  ,2499I5SE*00 


•249915S£*00  .4»08336'iE  •  00  1400*2«E-0  l 

12-COLUMN  13-COLUMn  | 4-COLUMN 


•  1 B 1 7 1 0 1 £-05  -.137424)£-q*  1446614E-0 1 
,7?22853£*00  .2863U0E-O1  -.P9S0566E-0J 
•528BS? I E  *00  - .99 1 958 ) C-02  -•  2240 1 57E-02 
,S083365£*00  -.140082BE-01  -.24461 37E-02 


Figure  42,  KONPACT-1  Output- -C-5A  Static  Elastic  Vehicle  Name  List 
Table  and  Quadruple  Data  (Cruise  Flight  Condition)  (See 
Figure  37  for  KONPACT-1  Input)  ( Concluded) 


.  agSMMuyi few  'kmm 


SYS1EM  NO  J 


VEHICLE  (  STATIC  ELASTIC  SYMMETRIC  -  REOUCED  ) 


NUMBER  OF  STATES  *  2 
NUMBER  OF  OUTPUTS-  A 
NUMBER  OF  INPUTS  -  1 

NAME  LIST  TABLE  ••• 

VARIABLE  NAME  DESCRIPTION  UNIT 

STATE 


1 

X  ( 

1> 

w 

VELOCITY  ALONG  z  AXIS 

INCH/ SEC 

2 

X( 

21 

0 

pitch  rate 

RADIAN/SEC 

OUTPUT 

1 

R  ( 

I  > 

SASGY 

PITCH  RATE  GYRO 

RADIAN/SEC 

2 

R< 

2) 

AZAP 

NORMAL  ACCELEROMETER 

INCH/SEC2 

3 

R  ( 

3) 

W 

velocity  ALONG  Z  AXIS 

INCH/SEC 

A 

R( 

<.» 

0 

PITCH  RATE 

RADIAN/SEC 

INPUT 

I 

U( 

I) 

BDELV 

elevator  oeflection 

RADIAN 

•••  QUADRUPLE  ORTA 


MATRIX  A  SIZE  -2X2 


I -COLUMN  2-COLUMN 


1- ROW  -.6TBG798E-00  ,»74|202E*0* 

2- ROW  -.187AM0E-03  1 101 I02E-01 


MATRIX  B  SIZE  -2X1 


I -COLUMN 


I  -ROW  -,33079*3E*03 
2-ROW  - . 1606369E»01 


Figure  43.  KONPACT-1  Output--Reduced  C-5A  Static  Elastic  Vehicle  Name 
List  Table  and  Quadruple  Data  (Cruise  Flight  Condition) 

(See  Figure  37  for  KONPACT-1  Input) 


i 

t 


L 


m 


MATRIX  C 


SIZE  m  U  X  2 


1 -COLUMN 


2-COLUMn 


1 - ROW  -.3585728E-05  .1033843E»01 

2- ROW  -. 54296 16E‘ 00  .6987265t»03 

3- ROW  . 1 Q03000E  *01  0. 

4- ROW  3,  .  1 000000E«01 


MATRIX  0 


SIZE  •  4  X  I 


i -column 


1 - ROW  -.6973532E-02 

2- ROW  .8325647E-03 

3- ROW  0. 

4- ROW  0. 


Figure  43.  KONPACT-1  Output--Reduced  C-5A  Static  Elastic  Vehicle  Name 
List  Table  and  Quadruple  Data  (Cruise  Flight  Condition) 

(See  Figure  37  for  KONFACT-1  Input)  (Concluded) 


■■ytexKmmma 


t 


SYSTEM  NO  2 


ACTUATOR 


r 


•••  TRANSFER  function  data  for  blocks  ••• 


BLOCK  I 


S»1  TERM  S*»0  TERN 
NUMERATOR  B.  7.50000 

DENOMINATOR  1.0*000  7.50000 

**•  CONNECTION  data  for  blocks  *•* 


MATRIX  P 

I -COLUMN 
I -RON  •. 


SIZE  >1X1 


MATRIX  0 


SIZE  >1X1 


1 -COLUMN 

I -RON  , 1#0*0*0E»01 

MATRIX  R  SIZE  >  «  X  I 

1 -COLUMN 

I -RON  ,1MBMBC*B1 


Figure  44.  Actuator  Transfer  Function  Data 

(See  Figure  37  for  KONPACT-1  Input) 
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:<  ^  * 


MATRIX  S 


SIZE  -1X1 


I -COLUMN 


I -ROW  «. 


Figure  44.  Actuator  Transfer  Function  Data  (Concluded) 
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ACTUATOR 


SYSTEM  NO  2 


NUMBER  OF  STATES  *  1 
NUMBER  OE  OUTPUTS*  1 
NUMBER  OF  INPUTS  *  1 


VARIABLE  name 

STATE 

1  X<  l» 

OUTPUT 
1  R  1  1  ) 

INPUT 

1  U  <  It 


•••  NAME  LIST  TABLE  ••• 

DESCRIPTION  UNIT 


DELE 

ELEVATOR 

deflection 

RADIAN 

OELE 

ELEVATOR 

deflection 

RADIAN 

DELEC 

ELEVATOR 

COMMAND 

RADIAN 

•••  QUADRUPLE  DATA 


MATRIX  A  SIZE  *1X1 

I -COLUMN 

1-ROW  -,7SMMQE»01 


MATRIX  B  SIZE  *1X1 

1 -COLUMN 

1-ROW  .lMOBMC’Ol 


Figure  45.  Actuator  Name  List  Table  and  Quadruple  Data 
(See  Figure  37  for  KONPACT-1  Input) 
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•  - 


■Tr'^-vi' 


MATRIX  C 


SIZE  ■  t  X  1 


I -COLUMN 

I -ROM  .75tMMC*tl 


MATRIX  0 


SIZE  >1X1 


I  -COLUMN 


1-ROM  0. 


Figure  45.  Actuator  Name  List  Table  and  Quadruple  Data  (Concluded) 
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SYSTEM  NO  3 


PILOT  HOOEL 


•••  TPANSFER  ruNCTION  DATA  FOR  BLOCKS  ••• 
BLOCK  1 

S**l  term  S*»0  term 
NUMERATOR  0.  . 223610E-03 

DENOMINATOR  1.00000  .100000 


CONNECTION  data  for  blocks  ••• 


matrix  p 


SIZE  *  I  X  I 


1 -COLUMN 


1-ROW  0. 


MATRIX  0 


SIZE  -1X1 


\ -COLUMN 

1-ROW  .lOOOOOOE’Ol 


MATRIX  R 


SIZE  >1X1 


I -COLUMN 

t-ROW  , I OOOOOOE'OI 


MATRIX  S 


SIZE  >1X1 


1  -COLUMN 


1-ROW  0. 

Figure  46.  Pilot  Model  Transfer  Function  Data 

(See  Figure  37  for  KONPACT-1  Input) 
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SYSTEM  NO  3 


PILOT  MOOEL 


number  or  states  «  i 

NUMBER  or  OUTPUTS*  I 

number  Of  inputs  »  1 


•••  NAME  LIST  TABLE 
VARIABLE  NAME  DESCRIPTION 

STATE 

I  X<  n  XP  PILOT  MODEL  STATE 

OUTPUT 

I  R (  II  FP  PILOT  COMMAND 

INPUT 

I  U<  I)  ETAP  PILOT  MODEL  INPUT 

QUADRUPLE  DATA 

MATRIX  A  size  »  l  X  I 

1 -COLUMN 

1-ROW  - .  lOOOOOOE'OO 

MATRIX  B  SIZE  *  I  X  1 

I -COLUMN 

1-ROW  , I #OOCOOE*01 

Figure  47.  Pilot  Model  Name  List  Table  and  Quadruple  Data 
(See  Figure  37  for  KONPACT-1  Input) 
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UNIT 

RADIAN 

RADIAN 

RADIAN 


jmiwmiwm*.. 


. 


1 -COLUMN 


1-ROW  ,?2361ME-#3 


MATRIX  0 


1 -COLUMN 


1-ROW  0. 


MATRIX  C 


SIZE  -IX  1 


SIZf  *  1  x  I 


Figure  47.  Tilot  Model  Name  List  Table  and  Quadruple  Data  (Concluded) 


SYSTEM  NO  4 


PLANT  (  PILOT  MODEL  ♦  ACTUATOR  •  VEHICLE  > 


•••  INTERCONNECTION  DATA  ••• 

j 

/ 

i 

/ 

MATRIX  R  SIZE  =  3X6 


1 -COLUMN 

2-COLUMN 

3-COLUMN 

X1 

4-<iDLUMN 

i 

5-COLUMN 

6-COLUMn 

1-ROW 

0, 

0. 

0. 

f 

0,i 

. 1000000E'0I 

0. 

2-ROV 

0. 

0. 

«. 

0. 

ft. 

0. 

3-ROW 

9. 

0. 

0. 

0. 

0. 

0. 

/ 


MATRIX  O 

SI2E  *  3  X 

1 -COLUMN 

2-COLUMN 

/ 

f 

1 - RON 

2- RON 

3- ROM 

0. 

0. 

.iftosaooE'Oi 

0. 

.iooeeooE*oi 

0. 

f 

/ 

MATRIX  R 

SIZE  *  5  X 

6 

l -COLUMN 

2-COLUMn 

J-COLUWN 

4-COLUMN 

5-COLUMN 

6-COLUMN 

I  -RON 

2- RON 

3- ROW 

4- R0W 

5- R0W 

.  1  ftOuftOOE  »  01 
0. 

«. 

0. 

0. 

0. 

• 1000000E»01 
0. 

0. 

9. 

0. 

0. 

.lOOOOOOE’OI 

0. 

0. 

0. 

0. 

0. 

,I000090E«01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

. 1 OftOOOOE-O 1 

MATRIX  S 

SIZE  *  5  X  2 

1 -COLUMN 

2-COLUMN 

l-RON 

«. 

ft. 

2-RON 

0. 

ft. 

3-RON 

0. 

0. 

4-RCN 

0. 

0. 

Figure  48.  Plant  (Pilot  Model  and  Actuator  and  Static  Elastic 
Vehicle)  Interconnection  Data  (See  Figure  37  for 
KONPACT-1  Input) 
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SYSTEM  NO  4 


PLANT  (  PILOT  MODEL  ♦  ACTUATOR  ♦  VEHICLE  ) 


NUM8ER  or  STATES  -  4 
NUMBER  or  OUTPUTS*  5 
NUMBER  or  Inputs  *  2 


VARIABLE  name 


»**  NAME  LIST  TABLE 

description 


I 

x< 

I  ) 

w 

VELOCITY  ALONG  Z  AXIS 

INCM/SEC 

2 

XI 

2) 

0. 

PITCH  Rate 

RADI  AN/SEC 

3 

x< 

3> 

DELE 

ELEVATOR  OEELECTION 

RADIAN 

4 

x< 

4) 

XP 

PILOT  MOOEL  STATE 

RADIAN 

OUTPUT 

1 

R  < 

1  > 

3ASGY 

PITCH  RATE  GYRO 

RADIAN/SEC 

2 

R  < 

z> 

AZAP 

NORMAL  accelerometer 

INCH/SEC2 

3 

R  ( 

31 

W 

velocity  along  z  axis 

INCH/SEC 

4 

R  ( 

4) 

0 

p  iterate 

RADIAN/SEC 

5 

R  ( 

5) 

rp 

PILOT  COMMAND 

RADIAN 

INPUT 

1 

U< 

n 

ETAP 

'X 

PILOT  MOOEL  INPUT 

RADIAN 

2 

U( 

2) 

DELEC 

ELEVATOR  COMMAND 

RADIAN 

MATRIX  A 


•••  OUAORUPLE  DATA  ••• 

SIZE  *4X4 

H  2-COLUMN  3-COLUMN 


1 -COLUMN 


1- GOW  -,6786798E»*0  .87412i)2E»04  -,2480957E*04 

2- ROW  I874C08E-03  -.I1011*2E*0I  1204776E*02 

3- ROW  0.  0.  -.7500000E*0l 

4- ROW  0.  0.  0. 


4-COLUMN 


0. 

0. 

0. 

-.I000800E*00 


Figure  49.  Plant  Name  List  Table  and  Quadruple  Data 
(See  Figure  37  for  KONPACT-1  Input) 


MATRIX  8 


1 -COLUMN 


SIZE  «  AX  2 
2-COL'JMN 


1- ROW  0.  0. 

2- ROW  0.  0. 

0.  •  I  OOOOOOE’0 1 

‘•-Row  .  iooooooe*oi  o. 


MATRIX  C 

SIZE  >  5  X 

4 

1 -COLUMN 

2-COLUMN 

3-COLUMN 

4-C01.UMN 

1-R0W 

-.3S8S728E-05 

. 103384 3E»01 

-.5227899E-01 

I 

0. 

2-ROW 

-. 54296 16C* 00 

,6987265E*03 

.6244235E*04 

0. 

3-R0W 

.lOOOOOOEOl 

0. 

0. 

0. 

4-R0W 

0. 

.lOOOOOOE'Ol 

0. 

0  • 

5-ROW 

0. 

0. 

0. 

.2236100E-03 

MATRIX  0 

SIZE  «  5  X 

2 

i 

| 

I -COLUMN 

2-C0LUMN 

1- R0W 

2- ROW 

3- ROW 

4- R0W 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

1 

Figure  49.  Plant  Name  List  Table  and  Quadruple  Data  (See 
Figure  37  for  KONPACT-1  Input)  (Concluded) 


.o. 

V. 


r 


SYSTEM  N 0  <9 


IMPLICIT  MODEL 


NUMBER  OF  STATES  *  2 
NUMBER  OF  OUTPUTS*  2 
NUMBER  OF  INPUTS  *  1 


VARIABLE  NAME 

STATE 


*•  name  list  table  ••• 
description 


matrix  a 


•••  OUADRUPLE  DATA  ••• 

SIZE  *2X2 
1  2-COLUMN 


MATRIX  B 


I  -ROM 
2-ROW 


I -COLUMN 


•,6786h«OE‘«0  .  874 1 200E  *04 

•.7S62000E-03  - . 3S21 30»E*01 


SIZE  *  2  X  I 


I -COLUMN 


■.33«7<940E*03 
■ , I 606370E,0 1 


I 

XI 

l) 

WI 

IMP 

MODEL 

VELOCITY 

FEET/SEC 

2 

X( 

2) 

01 

IMP 

MODEL 

PITCH  RATE 

RADI ANS/SEC 

OUTPUT 

I 

R( 

D 

WI 

IMP 

MOOEL 

VELOCITY 

FEET/SCC 

2 

R  < 

2.) 

01 

IMP 

MODEL 

PITCH  RATE 

RADI ANS/SEC 

INPUT 

1 

U< 

1  > 

DELEI 

IMP 

MODEL 

INPUT 

RADIANS 

Figure  50.  Implicit  Model  Name  List  Table  and  Quadruple  Data 
(See  Figure  37  for  KONPACT-1  Input) 


MATRIX  C 


1 - ROM 

2- ROW  0 


MATRIX  0 


1- ROW  0 

2- ROW  0 


Figure  50. 


SIZE  *2X2 
1 -COLUMN  2-COLUMN 

.i»oo#m*oi  o. 

•  19MtME*tl 

SIZE  >2X1 

I -COLUMN 


Implicit  Model  Name  List  Table  and  Quadruple  Data 
(See  Figure  3?  for  KONPACT-1  Input)  (Concluded) 


149 


&&&&$$  •**>*?. 


•  <>  •«•«••'»< 

0« 

o  a  mm 

• 

• 

•  SYSTEM 

NO  S 

OVERALL  SYSTEM  ( 

»LANT  ♦  IMPLICIT  MODEL  > 

* 

<* 

« 

»«»» 

**•  r  te- 

'CO'J’iecI  ION  D«rs 

MATRIX  D 

SIZE  * 

3  X  7 

l 

-COLUMN 

Z-COllIMN 

J-COLUMN 

4-CDLUMN  5-COLUMN 

ft-COLUMN 

7-COLUMN 

1  -ROW 

0. 

0. 

0. 

0.  0. 

0. 

o. 

2-R  OW 

o. 

0* 

0. 

o.  o. 

0. 

0. 

3-R0W 

o. 

0. 

0. 

0.  .  1 00000QE  *01 

0. 

0. 

MATRIX  0 

?E  *  3  X  2 

J-COLUMN 

COLUMN 

1  -ROW 

. iooooooe*o» 

0. 

2-ROW 

o. 

. 100000*E*01 

3-ROW 

0. 

0. 

MATRIX  P 

SUE  =  6  x 

7 

I -COLUMN 

2-COLUMN 

3-COLUMN 

4-column 

5-COLUMN 

6-COLUMN 

7-COLUMn 

1-ROW 

o 1 OOOOOOE  *01 

■>. 

0. 

0. 

0. 

0. 

0. 

2-ROW 

o„ 

.iooooooc*oi 

0. 

0. 

0. 

0. 

0. 

3-ROW 

o. 

0. 

•  1 1 31000E-01 

0. 

0. 

0. 

o. 

4-ROW 

0. 

0* 

0. 

0. 

.1000000E«01 

0. 

o. 

5-ROW 

G. 

0. 

•  1 OOOOOOE *01 

0, 

0. 

-.1000000E.01 

0. 

6-ROW 

0. 

0. 

•  1 OOOOOOE  *0 1 

0, 

0. 

-#1000000E*01 

MATRIX  S 

SIZE  >  *.  »  ? 

’-COLUMN 

7-COLUMN 

1-ROW 

0. 

0* 

2-ROW 

0. 

0. 

3-ROW 

0, 

0. 

‘•-QOt  0,  c. 

5- «0.,  0.  0. 

6- »3W  0.  0. 

Figure  51.  Overall  System  (Plant  and  Implicit  Model) 

Interconnection  Data  (See  Figure  37  for  KONPACT-1 
Input) 
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•  SYSTEM  MO  5  OVERALL  SYSTEM  <  PLANT  .  IMPLICIT  MOOEL  I 


NUMBER  OF  STATES  *  6 
NUMBER  OF  OUTPUTS*  6 
NUMBER  OF  INPUTS  *  2 


•••  NAME  LIST  TABLE  ••• 


VARIABLE 

:  name 

DESCRIPTION 

UNIT 

STATE 

l 

XI 

1) 

M 

VELOCITY  ALONG  Z  AXIS 

INCH/SEC 

2 

x( 

?> 

0 

PITCH  RATE 

RADIAN/SEC 

3 

X( 

1> 

OELE 

ELEVATOR  DEFLECTION 

radian 

A 

X  ( 

41 

XP 

PILOT  MODEL  STATE 

RADIAN 

S 

XI 

SI 

MI 

IMP  MOOEL  VELOCITY 

FEET/SEC 

6 

XI 

M 

01 

IMP  MODEL  pitch  rate 

radians/sec 

OUTPUT 

1 

R  1 

i  • 

SASOY 

pitch  rate  GYRO 

RAOIAN/SEC 

2 

R  1 

2) 

AZAP 

normal  accelerometer 

INCH/SEC2 

3 

R{ 

31 

alpha 

ANGLE  OF  ATTACK 

RADIAN 

A 

R  I 

A 1 

FP 

PILOT  COMMAND 

RADIAN 

5 

Rl 

SI 

EMI 

MOOEL  FOLL  ERROR 

FEET/SEC 

6 

R 1 

hi 

EOI 

MOOEL  FOLL  ERROR 

radians/sec 

INPUT 

1 

U< 

X 1 

ETAP 

PILOT  MODEL  INPUT 

RAO!  AN 

2 

Ul 

?) 

DEL  EC 

ELEVATOR  COMMAND 

RADIAN 

•••  QUADRUPLE  DATA  ••• 


MATRIX  A 

Size  *  6  X 

6 

1 -COLUMN 

2-COLUMN 

3-COLUMN 

a-column 

5-COLUMN 

6-COLUMN 

1-ROM 

-.6786798£«#0 

•87A129?£»8A 

-.2AB8957E*0A 

e. 

0. 

». 

2-ROM 

-.167A808E-83 

- • 1 101 I0?E»01 

-.1Z0A776E»02 

e. 

0. 

0. 

3-ROM 

0. 

0. 

-,T50#OOOE‘Ol 

0. 

0. 

0. 

A-ROM 

0. 

0. 

0. 

-.1900080E»00 

0. 

0. 

5-ROM 

0. 

1. 

0. 

-.7396B85E-01 

-,6786880E*00 

,47A1200E'0A 

6-ROM 

«. 

0. 

8. 

-.359280AE-03 

-.7562080E-03 

- , 3521 300E  *01 

Figure  52,  Overall  System  Name  List  Table  and  Quadruple  Data 
{See  Figure  37  for  KONPACT-1  Input) 
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1 -COLUMN 


2-COLUMn 


1  -ROW 

0. 

0. 

2-BOW 

0. 

0. 

3-ROW 

0. 

. 100000oE*01 

6-ROW 

.  1900000E*01 

0. 

5-O0W 

0. 

0. 

6-»0W 

0. 

0. 

MATRIX  C 

SUE  «  6  K 

6 

i -column 

2-COLUMN 

3-COLUMN 

4-C0LUMN 

5-COLUMN 

6-COLUMN 

1  -ROW 

-.358S72BE-0S 

•  J  033863E*0 1 

-.5227B99E-0I 

0. 

0. 

0. 

2-BOW 

-.S629616E-00 

.698726SE*03 

.6246235E*#4 

0. 

0. 

0. 

3-B0W 

. 1 1 31 000E-03 

0. 

0. 

0  . 

0. 

0. 

6-BOW 

0. 

0. 

0. 

•2236I09E-03 

0. 

0. 

5-P0W 

, 1 009890E' 0 1 

0. 

0. 

0. 

-. 1000000E*01 

0. 

6-BOW 

0. 

. I080000E*OI 

0. 

0. 

8. 

- • i OOOOOOE 

MATRIX  0 

SUE  *  6  X  2 

1  -COLUMN 

2-COLUMs 

1  -BOW 

0, 

0. 

2-R0W 

e. 

0. 

3-B0W 

0. 

0. 

fc-pow 

0. 

0. 

S-BOW 

0. 

9. 

6-BOW 

0. 

0. 

Figure  52.  Overall  System  Name  List  Table  and  Quadruple  Data 
(See  Figure  37  for  KONPACT-1  Input)  (Concluded) 
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,  *  > 


»  system  no  5  overall  system  t  design  msdel  t 

o 


NUMBER  OF  ST&TES  *  A 
NUMBER  OE  OUTPUTS® 1 0 
NUMBER  OE  INPUTS  =  P 


«««  name  LIST  TABLE  ••• 


VARIABLE  name 

DESCRIPTION 

UNIT 

STATE 

I 

X  (  1) 

w 

VELOCITY  along  z  axis 

INCH/SEC 

2 

X  (  at 

a 

pitch  RATE 

RADIAN/SEC 

3 

x  (  3) 

DELE 

ELEVATOR  DEFLECTION 

RADIAN 

A 

X  (  H) 

x° 

PILOT  MODEL  STATE 

RADIAN 

DESIGN 

OUTPUT 

l 

r<  n 

0 

PITCH  RATE 

RADIAN/SEC 

2 

R  (  21 

ALPHA 

ANGLE  or  ATTACK 

RADIAN 

3 

R  (  31 

DELE 

ELEVATOR  deelection 

RADIAN 

A 

R<  A) 

O/DT  OE 

(  DELE  ELEVATOR  DEELECT 

) 

RADIAN 

/SEC 

5 

R  <  51 

D/OT  OE 

(  EMI  MODEL  FOLL 

ERROR 

) 

feet/sec 

/SEC 

s 

R  (  6) 

fl/OT  OE 

<  E01  MODEL  FOLL 

ERROR 

> 

RADIANS/SEC 

/SEC 

SENSOR 

OUTPUT 

7 

R  (  71 

SASGY 

PITCH  RATE  GYRO 

RADIAN/SEC 

9 

R  (  #1 

AZAP 

NORMAL  ACCELEROMETER 

INCH/5EC2 

9 

R  (  01 

dele 

ELEVATOR  DEFLECTION 

RADIAN 

10 

P<10> 

EP 

PILOT  COMMAND 

radian 

CONTROL 

INPUT 

1 

U(  II 

OELEC 

ELLVATOR  COMMAND 

RADIAN 

commanc 

1  INP11T 

2 

111  21 

ETAP 

PILOT  MODEL  INPUT 

RADIAN 

•••  OUADRUPLE  DATA  ••• 


matrix  a 

SIZE  *  AX  4 

1 -COLUMN 

2-COLUMN  3-COLUMN 

4-COLUMN 

Figure  53.  Overall  System  (Design  Model)  Name  List  Table  and 
Quadruple  Data  (See  Figure  37  for  KONPACT-1  Input) 
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MATRIX  8 


1  -ROW  -  .f>786798E»00  ,8T4120?£*04  -,24B0957E»O4  0. 

2- ROW  -.1874000E-03  -.1I0U0?E»01  «.  1 204776E‘02  0. 

3- ROW  0.  0.  -,7500000E*01  0. 

4- ROW  0.  0.  0.  -.1000600E*00 


I -COLUMN 


0. 

0. 

.  I900000E-0! 

0. 


MATRIX  C 


SIZE  =4X2 


2-COLUMN 


0. 

0. 

0. 

.  IOOOOOOE.01 


SIZE  ■  10  X  4 


I -COLUMN 


2-COLUMN 


3-COLUMN 


4-COLUMN 


1-R0W 

6. 

• 1660066E*61 

6. 

0. 

2-ROW 

, 1 1 31 600E-63 

0. 

0. 

6. 

3-ROW 

0. 

6. 

. 1000000E*0 I 

0. 

4-ROW 

0. 

0. 

- • 7500000E*01 

0. 

5-ROW 

. 2609666E-66 

.21 00606E-62 

-.2486957E*64 

•7396B85E-0I 

6-ROW 

. 5688660E-03 

•242619flE*0i 

-,1264776E*62 

.3592664E-03 

7-ROW 

- . 3565728E-65 

. 1033B43E»0I 

-.5227899E-61 

6. 

6-ROW 

-,54?9616E*60 

.6987265E«03 

.6244235E*04 

6. 

9-ROW 

0. 

0. 

•  1 OOOOOOE’Ol 

6. 

16 -ROW 

6. 

0. 

6. 

.2236160E-03 

IATRIX  0 

SIZE  «  10  X 

2 

1 -COLUMN 

2-COLUMN 

1-R0W 

6. 

0. 

2-ROW 

6. 

0. 

3-ROW 

0. 

0. 

4-R0W 

, I666600E»01 

6. 

5-ROW 

6. 

0. 

6-ROW 

0. 

0. 

7-ROW 

e. 

0. 

8-ROW 

0. 

0. 

9-ROW 

0. 

6. 

16 -ROW 

0. 

6. 

Figure  53.  Overall  System  (Design  Model)  Name  List  Table  and 
Quadruple  Data  (See  Figure  37  for  KONPACT-1  Input) 
(Concluded) 


TOOAyXS  OATE  AUG  I  7.  75 


IDENTIFICATION 


J  K  MAHESM 


MAX  DUMBER  OF  INNER-LOOP  ITERATIONS  40  MAX  NUMBER  OF  OUTER-LOOP  ITERATIONS  30 
MAX  NUMBER  OF  ITERATIONS  ON  ELIMINATING  CONTROL  SURFACE  FEEDBACKS  4 


MF..,  5  -CHL1. 

'JO r3V  “  1 

'ISTFP  *  1 

N=RI':  -  i 


»  I  T  !  •  -’•) 

ioi  n  r  !z 


'ICON  T  =  0 


FLK-mt  C3.NO! TICK. 


RUN  > 


OROE  •  OF  SYSTEM  -  * 

NllM-rO  Of  PFSOQXyS  =  t 
NUMBER  OF  cONTiOL*1  -  1 

NUMriPR  Or  DlST:JRPf.'jCc  INPUTS  =  ( 

NDPO.-O  OF  r«  is^,  STATES  r  3 
NUMBER  0C  GUST  HrtjTs  s 
NUMBER  Of  COM 'A NO  sTaT-S  =  i 

N’JMMTR  Or  GUST  rrT  OR  OUT*  STATES  *  e 

CONTROL  R«Tf  RpSOo'jSeS  STa:T  rlTr*  RES°0NQE  7 


STAT'S  are  owriswETT  AS 
1  2  3  4 


F  MATRIX 


ROW  I 


-.67B6BE*  iC 

.B74i?E*U4 

-,24«IOE*u« 

0. 

ROW  2 

-»IB740E-!’3 

-.1101 1 E  *  0 1 

- ,  1 2048F.D2 

0. 

ROW  3 

0. 

0. 

-, 75000E  *0 1 

0. 

ROW  4 

0. 

0. 

0. 

-.lOOOOE’OO 
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Figure  55.  KONPACT-2  Output  (Employing  FFOC  to  Compute  Reduced 

Feedback  Gains)  for  Static  Elastic  Design  Model  (See  Figure  39 
for  KONPACT-2  Input)  (Continued) 
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-.11037E-08  .35502E-09 

.39849E-04 

-. 16769E*01 

-.13807E-03  .39849E-0A 

.50000E*01 

R.M.s.  RESPONSES 

1  .6795501  IE-04 

2  .8831 13O4E-04 

3  , 1B842019E-04 

4  .6A980631E-04 

5  . 1 221 3247E*00 

6  .737255O7E-04 

TOTAL  R.M.S.  RESPONSES 

1  .6T955011E-04 

2  .883 1 1 304E-04 

3  .18S42019E-04 

4  .6498C631E-04 

5  .  I?2l3?47£«00 

6  .73725507E-84 


EIGENVALUES 

REAL  IMAGINARY 


0AMPIN3  RATIO  EREOUCNCV 


-.10000000 

-1  0.621:68367 

-2.24I880S5  2.16364058  -.71956724  3.1155956$ 


NEW  PREDICTOR 


ROW  1 

-. 16090E-0 1  .84593E-0S  0.  -.16729E-01 


KONPACT-2  Output  (Employing  FFOC  to  Compute  Reduced 
Feedback  Gains)  for  Static  Elastic  Design  Model  (See  Figure  39 
for  KONPACT-2  Input)  (Concluded) 
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TOOAYXS  DATE  AUG  I  7 *  75 


IDENTIFICATION 


J  K  MAHESH 


s 


TIME  RESPONSES  PLOTTING  TIME  *  S.000 

SAMPLE  INTERVAL  *  .2000E-01 

PLOTTING  SAMPLE  INTERVAL  =  .2000 

F1»ST  DELAY  TIME  a  0. 

SECOND  DELAY  TIME  =  0. 


PLOTTING  VARIABLES 


RESPONSE  NUMBER 


RESPONSE  VARIABLE  RESPONSE  OMITS  MIN  SCALE  MAX  SCALE 


scale  factor 


t  a 

2.  ALPHA 

0  DELTAE 


-.0 

-.0 

-.0 


-.0 

-.0 

-.0 


1.00 

1.00 

1.00 


MAX  NUMBER  of  INNER-LOOP  ITERATIONS  AO  RAX  NUMBER  OF  QUTER-LOOP  ITERATIONS  30 
MAX  NUMBER  OF  ITERATIONS  ON  ELIMINATING  CONTROL  SURFACE  FEEDBACKS  A 


NEW  PROBLEM  WITH  INPO  *  1  INPK  *  l  NCONT  •  1 

NOCOV  *  3 

NSTEP  *  0  NRANO  a  0 
NPRIN  *  0  NPLOT  *  0 


FLIGHT  CONOITION  AS  RUN  I 


ORDER  OF  SYSTEM  a  A 

NUKBtR  OF  RESPONSES  «  b 

NUMBER  OF  CONTROLS  »  1 

NUMBER  OF  DISTURBANCE  INPUTS  «  1 

NUMBER  OF  FEEDBACK  STATES  a  A 

NUMBER  OF  GUST  INPUTS  ■  0 

NUMBER  OF  COMMAND  STATES  ■  0 

NUMBER  OF  GUST  LIFT  GROWTH  STATES  «  0 

CONTROL  RATE  RESPONSES  START  WITH  RESPONSE  7 


STATES  ARE  ORDERED  AS 
I  2  3  A 


I 


Figure  56.  KONPACT-2  Output  (Employing  DLAK  to  Evaluate  Time  Responses 
to  Elevator  Command)  for  Static  Elastic  Design  Model  (See  Figure 
40  for  KONPACT-2  Input) 


F  MATRIX 


ROW  1 

*.67868E*00 

. 874 1 2E ♦ 04 

ROW  2 

- . I8740E-03 
ROW  3 

1 1 01 1 E*01 

0. 

ROW  4 

0. 

0. 

0. 

gi  matrix 


row  i 

o. 

ROW  2 

0. 

ROW  3 

, lOOOOE'O l 
ROW  4 
0. 


”.248I0e*04  0. 
-.1204BE«02  0. 

7S000E*0l  0. 

0.  1 OOOOE  *00 


32  MATRIX 


ROW  I 

0. 

ROW  2 

0. 

ROW  3 

0. 

ROW  4 

• tOOOOE’Ol 


INITIAL  CONDITION  MfttP-I* 


ROW  ] 

0. 

ROW  2 

0. 

ROW  3 

0. 

ROW  4 

0. 

I 

i  i 


Figure  56.  KONPACT-2  Output  (Employing  DLAK  to  Evaluate  Time  Responses 
to  Elevator  Command)  for  Static  Elastic  Design  Model  (See  Figure 
40  for  KONPACT-2  Input)  (Continued)  ’ 
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STATE  limit  -  RATE  1 

ROW  1 

.  1 0000E  *2i 
ROW  2 

. IOOOftt-20 

•  1 0000E*  20 
ROW  3 

.  1  OOQOE.*20 

•  1 OOOOE  *20 
ROW  <. 

. IOOOOE‘20 

. I0000E>20 

. 1 0000E  *20 

COMMAND  level  MATRIX 


ROW  I 

.44716E*  4 


h  matrix 


ROW  1 

0. 

. 1 OOOOE*0 t 

0. 

0. 

ROW  ? 

.U3IOE-0-1 

0. 

0. 

0. 

ROW  3 

0. 

0. 

.  10000E‘01 

0. 

ROW  4 

O. 

0. 

-,75000E»0I 

0. 

ROW  5 

.20000E-04 

•21000E-02 

-.2i*8l<»E.*04 

.73R6RE-0I 

ROW  6 

.S6880E-03 

,24202E«01 

-.U04Be»02 

.35920E-03 

o  matrix 


row  l 

o. 

ROW  2 

0. 

ROW  3 

0. 

ROW  4 

.  1 ©000E*01 

ROW  5 

0. 

ROW  6 

0. 


Figure  56.  KONPACT-2  Output  (Employing  DIAK  to  Evaluate  Time  Responses 
to  Elevator  Command)  for  Static  Elastic  Design  Model  (See  Figure 
40  for  KONPACT-2  Input)  (Continued) 
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imtmstsimmm 


M  «*?R1X 


ROW  I 

-.35637E-  5  .10338C*0I 

ROW  2 

-.542966*00  .69873f03 

ROW  3 

0.  0. 

ROW  <\ 

0.  0. 


-.52279E-01  0. 
.6£442E*04  0. 
.lOOOOfOl  0. 

0.  .22061E-03 


INPUT  (MINS  MATRIX 


ROW  l 

•  3 1 1 69£* 


-•95075E-03  -.901<j]g.00  .169396*01 


aUflORATfC  WEIGHTING  matrix 


ROW 

I 

0. 

ROW 

2 

0. 

0. 

0. 

«. 

O. 

0. 

ROW 

3 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

4 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

5 

0. 

0. 

0. 

0. 

0. 

0. 

ROW 

6 

Q. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

Figure  56.  KONPACT-2  Output  (Employing  DIAK  to  Evaluate  Time  Respons, 
to  Elevator  Command)  for  Static  Elastic  Design  Model  (See  Fieui 
40  for  KONPACT-2  Input)  (Continued)  gU1 
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AIRCRAFT  RESPONSES  WITH  PRESCRIBED  gains 


TIME  responses 

AUG  ,7’  75  j  *  *AH&r  ,  FLIGHT  CONDITION  45  RUN  1 

THERE  ARE  3  RESPONSES  TO  CQHPUTt 


TIME  RESPONSES  FOR  DISTURBANCE  1 


time  - 

0.000 

a 

=  0. 

ALPHA 

*  0. 

OELTaE 

X 

.abe-oi 

time  * 

.200 

a 

=  -.17 

ALPHA 

*  -.21E-01 

DELTAE 

X 

•81E-01 

TIME  * 

.4oo 

a 

=  -.27 

Alpha 

=  -.64E-01 

DELTAE 

m 

.S2E-0 1 

TIME  * 

.600 

a 

=  -.28 

Alpha 

=  -.11 

DELTAE 

X 

• 32E-0 1 

time  * 

.doo 

a 

*  -.24 

alpha 

*  -.15 

DELTAE 

■ 

.24E-01 

TIME  = 

1.000 

a 

=  -.20 

alpha 

*  -.17 

DELTAE 

s 

.23E-01 

TIME  = 

1  .200 

a 

=  -.16 

ALPHA 

=  -.18 

DELTAE 

X 

.25E-01 

TIME  = 

1  .400 

a 

=  -.13 

Alpha 

=  -.10 

DELTAE 

X 

.29E-01 

TIME  * 

1.600 

a 

=  -.11 

alpha 

=  -.10 

DELTAE 

X 

.32E-01 

TIME  * 

1.800 

a 

=  -.11 

ALPHA 

*  -.18 

DELTAE 

X 

.34E-01 

time  = 

2.000 

a 

=  -.11 

ALPHA 

*  -.18 

DELTAE 

X 

.35E-01 

time  = 

2.200 

a 

=  -.1  1 

alpha 

(  =  -.18 

OELTAE 

X 

.35L-01 

TIME  * 

2.400 

a 

=  -.1  1 

alpha 

»  -.18 

OELTaE 

X 

.35C-01 

Figure  56.  KONPACT-2  Output  (Employing  DIAK  to  Evaluate  Time  Responses 
to  Elevator  Command)  for  Static  Elastic  Design  Model  (See  Figure 
40  for  KONPACT-2  Input)  (Continued)  (See  Figure  33  for  On-Line 
Time  History  Plots) 
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TIME  = 

2.600 

Q 

=  -.11 

ALE-iA 

=  -.18 

DELTAE 

S 

.35E-Q1 

TINE  = 

a. eoo 

a 

=  -. 1 1 

ALPHA 

=  -.  18 

OELTaE 

S 

• 35E-0 1 

TIME  = 

3.600 

a 

=  -.11 

alpha 

*  -.16 

DELTAE 

a 

•35E-01 

time  * 

3.200 

a 

=  -.11 

alpha 

=  -.18 

OELTaE 

c 

.35E-01 

TIME  = 

3.4fft 

a. 

=  -.11 

ALPHA 

=  -.18 

DELTAS 

a 

.3SE-01 

time  * 

T.60C 

a. 

=  -.11 

ALPHA 

=  -.16 

DELTaE 

a 

•35E-01 

time  = 

3.000 

a 

=  -.11 

alpha 

»  -.18 

DELTAE 

9 

.3SE-01 

TIME  * 

4.000 

0 

=  -.11 

alpha 

*  -.18 

OELTaE 

s 

•35E-01 

time  * 

4.200 

0 

=  -.11 

alpha 

*  -.18 

DELTAE 

a 

•35E-01 

time  = 

4.400 

a 

=  -.11 

alpha 

*  -.18 

DELTaE 

a 

. 35E-01 

time  * 

*.600 

a 

=  -.11 

ALPHA 

«  -.18 

OELTAE 

a 

• 38E-01 

TIME  * 

4.600 

a 

=  -.11 

ALPHA 

=  -.18 

DELTAE 

a 

•35E-01 

TIME  = 

5.000 

0 

=  -.11 

ALPHA 

=  -.10 

OELTAE 

a 

•3SE-01 

Figure  56.  KONPACT-2  Output  (Employing  DIAK  to  Evaluate  Time  Responses 
to  Elevator  Command)  for  Static  Elastic  Design  Model  (See  Figure 
40  for  KONPACT-2  Input)  ( Continued)  (See  Figure  33  for  On-Line 
Time  History  Plots) 
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Figure  59.  Gust  Model  Block  Diagram 


SUBROUTINE  SIHK2IXOOT. 

Y.X,U»XDOTL» YLtRL*NX.NY*NRtNU»iNIT.T> 

SINK2 

2 

c 

SIMK2 

3 

c 

PURPOSE  -  TO 

implement 

SIMULATION  EQUATIONS  FOR  CSA  CONTROLLER 

SIMK2 

4 

c 

ANALISIS  -  A 

F  KONAR  / 

J  K  MAMESH  -  THE  HONETMELL  INC 

SIMK2 

5 

c 

OATE  WRITTEN 

-  197S 

SIMK2 

6 

c 

SIMK2 

7 

c 

ARGUMENTS  LIST 

SIMK2 

8 

c 

XDoT 

ARRAY  FOR  STATE  DERIVATIVES 

SIMK2 

9 

c 

Y 

ARRAY  FOR  Y  EQUATIONS 

SIMK2 

10 

c 

X 

ARRAY  FOR  STATES 

SIMK2 

11 

c 

U 

ARRAY  FOR  EXTERNAL  INPUTS 

SIMK2 

12 

c 

XOOTL 

OUTPUT 

ARRAY  FOR  DERIVATIVE  OF  STATE 

SIMK2 

13 

c 

YL 

OUTPUT 

ARRAY  for  y  equation  variables 

SIMK2 

16 

c 

KL 

OUTPUT 

ARRAY  FOR  EXTERNAL  RESPONSE  VARIABLES 

SIMK2 

15 

c 

NX 

OUTPUT 

NUMBER  OF  STATES 

SIMK2 

16 

c 

NY 

OUTPUT 

NUMBER  OF  Y  EQUATIONS 

SIMK2 

17 

c 

NR 

OUTPUT 

NUM8ER  OF  OUTPUTS 

SIMK2 

18 

c 

NU 

OUTPUT 

NUMBER  of  inputs 

SINK  2 

19 

c 

INI  T 

INPUT 

INITIAL  MODE  FLAG 

SIMK2 

20 

c 

T 

OUTPUT 

sample  time 

SIMK2 

21 

c 

SIMK2 

22 

DIHENs ION  XDOT(NX) .Y(NYt iXCNXt  «U <NU 1 1 XDOTL <NX I , YL (NY > .RL (NR) 

SIMK2 

23 

CONNON  /INOUT/  IR. IWt IPRINT %  INSERT 'LOCATE «NULL«MARK (20) . JN» JO. JS 

SIMK2 

26 

DIMENSION  CARO (20) 

SIMK2 

25 

REAL  WNI.KN2 

.KAFtKO.KP 

♦KNF,MLCI*MLC2 

SIMK2 

26 

OATA  hEND8»HVITC»HAINB/<.MEN0  .4H*ITC« 4HA1N  f 

DATA  hMLC1.HMLC?.HSASR*HAL0C/4HMLC1*<*NNLC2*4HSaS  .ahaldc/ 

DATA  mKMIB«MKM2B»HK AF0 » HK0BR/4HKM 1  i4HKW2  tAHKAF  «  4HKQ  / 

DATA  HKPBB.HKNFB/4HKP  .4MKNF  / 

CHECK  IF  INITIALIZATION  MOOE 

irilNlT.NE.*)  DO  TO  100 

SET  FILTER  GAINS 

AP—.l  S  BP-.22361E-03 

ANF--4.0 

AF«-.f)2 

AMl«-.01 

AM2— .01 

AHF«-1.  >  BHFa-1. 

ATF—4.0  t  BTFmA.O 

SET  CONTROLLER  SNITCHES 

SAS-0.0  l  ALDCS-0.0  *  MLC1>0.0  §  MLC2-0.0 

SET  CONTROLLER  GAINS 

KMl«l,0/0.26 

KM2-I.0/#. #5«9l 

KAF»36.0#0.26 

KO«O.S 

KP»0«  7068 

KNFa-0.09 

PEAD  CONTROLLER  SNITCHES  ON  AND  CONTROLLER  GAIN  VALUES 

CONTINUE 
READ! IR»28) CARD 
EORHAT (20AAI 

IF(CAPOU)  .EO.HENOBtGO  TO  8c 

Figure  61.  Subroutine  SIMK2  Program  Listing 


SIMK2  27 
SIMK2  28 
SINK2  29 
SIMK2  3« 
SIMK2  31 
SIMK2  32 
SINK2  33 
SINK2  3A 
SINK2  35 
SIMM  36 
SINK 2  37 
SIMM  38 
SIMM  39 
SINK2  A6 
SIMK2  A I 
SINK2  A2 
SINK2  A3 
SINK2  66 
SIMM  65 
SINK2  66 
SINK2  67 
SIMM  66 
SINK2  69 
SINK2  59 
SINK2  51 
SINK2  S2 
SIMM  53 
SINK 2  56 
SINK2  55 
SINK 2  56- 
SIMM  57 
SINKS  56 
SINKS  59 
SINKS  66 
SINKS  *1 
SINKS  M 
SINKS  6} 
SINKS  96 


f 


IF (CARO <41 .NE.H* ITC)GO  TO  4o 

SIMK2  65 

c 

SIMK2  66 

c 

READ  CONTROLLER  SWITCHES  ON 

SIMK2  67 

c 

SIMK2  68 

30 

CONTINUE 

SIMK2  69 

READ  <  JR*2C 1  CARO 

SIMK2  70 

trtCARom .eq.hendrigo  to  io 

SIMK2  71 

IF<CA»0<1» .EQ.HMLCl»MLCl«l.i 

SIMK2  72 

IF<CAPO(l» .EO.HHLC1IOO  TO  3c 

SIMK2  73 

IF (CARO I 1  * •EO»HNLC?)MLC2*l . i 

SIMK2  7<* 

IF<CAR0<1 J .EO.HMLCPIGO  TO  30 

SIMK2  75 

IF (CAPO! 1 > .EO.HSASR)SAS«1.0 

SIMK2  76 

!F(CAPO<l>*EQ.HSASR)GO  TO  3fi 

SIMK2  77 

IF (CAPO < 1 ) .EO.HALDC) ALDCS»1 .0 

SIMK2  78 

IF (CAPO ( I ) .EQ.MALOC) GO  TO  3o 

SIMK2  79 

STOP  III 

SIMK2  60 

c 

SIMK2  81 

c 

RE AO  CONTROLLER  GAIN  VALUES 

SIMK2  82 

c 

SIMK2  83 

40 

CONTINUE 

SJMK2  64 

IF (CARO (4 >  .NE.HAINRlSTOP  111 

SIMK2  85 

so 

CONTINUE 

SIMK2  86 

REAO<IR'20)CARO 

SIMK2  87 

IF  (CAPO  ( 1 )  .EO.HENOR)  GO  TO  i{. 

SIMK2  88 

IF(CApO(1» .EO.HKH|R)PfAO(IR.60)KHl 

SIMK2  89 

60 

FORMAT <E12.6» 

SIMK2  90 

IF (CARO ( 1 ) .EQ.HKM 1R) GO  TO  5o 

SIMK2  91 

IF (CAPO ( 1 1 .EO.HKM2r)ReAO(IR«60>KM? 

SIMK2  92 

IF(CAPOU)  .E0.HKM2BJ00  TO  5j 

SIMK2  93 

IFICAPOm  .EQ.HKAFR)REAOIIR«60>KAF 

SIMK2  94 

IFICApDU  )  .EO.HKAFRIGO  TO  So 

SIMK2  95 

IF (CAPO (1) .EO.HKQBR>REAO(IR*60VKO 

SIMK2  96 

IFICAPOIl ) .eo.hkrbrigo  TO  So 

SIMK2  97 

IF(CAROd)  .EO.MKPBR)READ»IR»60>KP 

SIMK2  98 

IFICAPOm.EO.HKPBBIGO  TO  5ft 

SIMK2  99 

IFICAROin  .EQ.HKNFR|REAOUR«60»KNF 

SIMK2100 

IF(CAQO(l».EO.HKNFR}GO  TO  so 

SIMK2101 

STOP  111 

SIMK2102 

80 

CONTINUE 

SIMK2103 

c 

SIMK2104 

c 

SET  DIMENSIONS  OF  SYSTEM 

SIMK2105 

c 

SIMK2106 

NX»7  S  NR*3  t  NU*9  S  NY»5 

SIMK2107 

c 

SIMK2108 

c 

return 

SIMK2109 

c 

SIMK21 10 

return 

SIMK2111 

c 

SIMK2112 

c 

SIMULATION  EQUATIONS 

SIMK2113 

c 

SIMK2114 

100 

CONTINUE 

SIMK2115 

c 

SIMK2116 

c 

differential  equations 

SIMK2117 

c 

SIHK21 18 

XDOTL  ( 1  >  »AP*X  <  1 )  ♦BP*U'(  3 > 

SIMK21 19 

XOOTL <  2) >ANF«X  <21 ♦ ALOCS*U (61 

SIMK2120 

XDOTL (3) *AM1*X 131 *MLC1#YI2) 

SIMK2121 

XOOTLT4) >AF*X (4 1  *ALDCS*Y  1 1 1 

S1MK2122 

XDOTL <5*«ATF»X<5»*RTF*Y<3» 

SIMK2123 

XD0TL(6)»AHF*X<6)»BHF*U<2> 

SIMK2124 

XDOTL ( T ) ■ AM2«X  <  7 1 *MLC2» Y  <4 ) 

SIMK2125 

c 

SIMK2126 

c 

SUMMING  POINT  EQUATIONS 

SIMK2127 

c 

SIMK2128 

YL<1>»KAF*X<2$*ANF*UU> 

SIMK2129 

YL<2l9BN»»U.(AV«!«!?» 

SIMK2130 

Figure  61.  Subroutine  SIMK2  Program  Listing 

(Continued) 
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no  o  non 


VL(3)»KP*X(1)*X(6) *U<2) *KNF  *U ( 7 ) 

SIMK2131 

YL<4)»KM2*U<5)-U<9t 

SIMK2132 

YL  <5) =*LOCS*X (5) *X ( 1 > «SAS«K0*U<8) 

SIMK2133 

SIMK2134 

RESPONSE  equations 

S IMK2 l 35 

SIMK2136 

RL(l)cUU) 

SIMK2137 

RL(2)«Y(5) 

SIMK2138 

RL (3>»X  ( I ) 

SIMK2139 

SIMK2140 

RETURN 

S1MK214I 

SIMK2142 

return 

SIMK2143 

END 

SIMK2144 

Figure  61.  Subroutine  SIMK2  Program  Listing  (Concluded) 
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Figure  62.  Reduced  ALDCS  Controller  Block  Diagram 


Response 


Weight 


Value 


MLC1 

0.  800E+01 

B1 

«2 

0. 100E-01 

T1 

% 

0. 100E-08 

qs 

% 

0 

B2 

<*5 

0 

T2 

0  ! 

6 

a 

«*7 

0.  500E+04 

B3 

S 

0 

T 

3 

S 

0 

6  . 
ei 

Q10 

0.  600E+06 

B4 

«U 

0 

T4 

Q12 

0 

6 

a 

Q13 

0 

B5 

Q14 

0 

T5 

Q15 

0 

6ei 

Q16 

0 

B1 

Q17 

0.  750E-13 

*1 

Q18 

0. 100E-10 

Q19 

0 

Figure  64.  Design  Response  Weights  for  ALDCS  Controller  Design 
(C-5A  Cruise  Flight  Condition) 
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Figure  64.  Design  Response  Weights  for  ALDCS  Controller  Design 
(C-5A  Cruiss  Flight  Condition)  (Concluded) 
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Figure  65.  Reduced  Feedback  Gains  for  ALDCS  Controller  Design 
(C-5A  Cruise  Flight  Condition) 
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Figure  67.  Alpha  Response  of  C-5A  Open  Loop 
F24RR  Model  to  Elevator  Command 
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Figure  70.  Pitch  Rate  Response  of  C-5A  Open  Loop  F24RR 
Model  to  Elevator  Command 
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Figure  71.  Pitch  Rate  Response  of  C-5A  SAS  F24RR  Model 
to  Elevator  Command 
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Figure  72.  Pitch  Rate  Response  of  C-5A  ALDCS  F24RR  Model 
to  Elevator  Command 
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Figure  73,  Covariance  (RMS)  Comparison  of  Open  Loop  and  Closed  Loop  ALDCS  (Reduced 
Feedback)  Residual  Elastic  Vehicle  (C-5A  Cruise  Flight  Condition) 
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-,53S7E*04- 
-.24737*04- 
-.24817*03 
.71497*03- 
-.22987*01 
-.39447*01- 
.50777*04 
.76257*03- 
-.16987*04 
-.30267*03 
-.13177*43 
.6B03E*0) 
.64807*00- 
-.4524E«01 
-.24017*03 
-.36837*04- 
.22627*02 
-.53177*02 
.32567*03- 
-.95767*00 
-.63617*01- 
.21617*04 
.53887*03- 
,38337*02 
-,1S65E*03 
•?5?IE*03- 
.4827E*01 
.13137*01- 
.70147*00 
.39597*0? 
-.14657*04- 
. B2SIE*02 
-.66977*02 
-.74477*02- 
L/UE2 
.62467*01 
-.75687*00- 
-.29277*01- 
-.85767*03 
-.3985E‘«3 
-.31647*64 
.14647*63- 
-.56677*02 
-• 1922E*0I 
-,79867*00 
•66527*04 
-.18477*83- 
.3I72E*B3 
-.15367*03 


.72107 
,??Mt 
,952  IE 
.56887 
.1B25E 
, 1 5B9E 

.3802E 
.1 35CE 
.61737 
.50967 
,96107 
15 

. » «9lt 
.29047 
,54337 
.43767 
,  321  BE 
■  1 0457 
.7455E 
,?3*1£ 
.II02E 
.5B55F 
»  3665E 
, 36B5E 
.974BE 
.12507 
.21507 
.6370E 
.3I45E 
. 1360E 
.  I 399E 
, 62907 
,69697 
•  I B99£ 
.4267E 

Ml 

,31837 
,2423 1 
.4A90E 
.7861E 
.754 3E 
.1156E 
. 5662E 
.6497E 
.I273E 
. 3168E 
,2983E 
•2903E 
.4974E 
15 

.4944E 
, 2887E 
, 1993E 
, 1 883E 
.17487 
, 2266E 
, 1582E 
.25497 
,31627 
.14567 
.23397 
.82637 
.45827 
.65797 


•  1  >'f*t  ,  06- 
.97377.04 

-  . 14  76E  *05 

-  .  1 5767*04- 
-.21667.03 

.28367*03- 

.3?4?E*C» 

-.38657.66 

-.8678E.PS 

-.70417*04 

.88407*04- 

-.65127.03 

.15497*02. 
.18197.01- 
.67287*00 
.10977*06 
.11357*05- 
.75147*03- 
-.11167*04 
.26787*03 
.14247*02 
.16197*08- 
. 33027*04- 
.57637*04- 
.34457*03- 
.17997*04- 
.16907*04 
.51667*01- 
.52637*01- 
.67417*00 
.44607  *0* 
.37817.04- 
.12957*04- 
-.19627*03- 
-.11427*0! 
.90?9E*01 
.  1 195E.92- 
.61267*0.3- 
.38047*04- 
-.55977*02 

.10447*04- 

.63457*0) 

.32477*00 

.9230E.01- 

-.22927*00 

.14697*04 

-.24467*03 

.56547*03- 

-.54097*03 

.34857*02 

.16327*01- 

.I7T7E.8I 

.64497*00 

-.57127*03- 

.41557*03- 

.58327*03 

.45487*83- 

.14987*03 

-.28187*01- 

.21697-01- 

.22237*03- 

-.22187*03- 

-.15917*03- 

.70477*03 


.12067*06- 

. 1 65 17 ♦  04  ■ 

,744tC*03- 

.1)037*04- 

.*«54E*02- 

,3.-.337*03- 

,62S7£*04- 

,99777*04- 

.36667*04 

,  lr-567  *04- 

. I  9767*04 


,298'tE*Ul 

.3767E.IM 

.57287*04 

.70247*03 

.232liE*04 

,240oE*0>' 

.13517*04 

.56447*01 

•  I  3877*00* 

.36817*01 

.74267*03 

.4  778E  *04' 

.34547*03 

,399)7*02' 

.10507*03' 

.13767*01 

. 503*7*01' 

.33997*04 

.64627*03' 

.63547*03 

.20497*03 

.63257*02 

.64347*01 

.92997*00 

.16247*0) 

.36997*02 

.26797*04 

,67937*02 

.97497*0? 

.68147*0?' 

.18877*00 

.64127*01' 

.10447*04 

.32567*03' 

.14707*03 

.97387*0? 

.24047*03' 


.17767*01 

.43657*00- 

.79507*04 

.36677*03- 

.31987*03 

.11207*03 

.19187*03 

.90607*01 

,7?99E*00- 

.27557*01- 

.43547*03 

.26*11*83 

.42607*03 

.28847*0?* 


•  340 17*95' 
.660  •  7*04' 
.43237*04 
.10367*05' 
.17527*03- 
.10367*03' 
.24927*05' 
.39667*05 
.149?e*06 
,63757*04. 
.48517*04' 


*  1 544E*0? 
•.76987.01 
.686*7*04 
'.40887*04 
.10017*05 
.35067*03- 
.12607*04- 
.12977*02 
.24747*01 
.64477*01 
.35647*04 
.20407*04 
.8787E*03 
.11717*03 
.27007*03 
.57.337*01 
.36087*01 
.11977*04 
.30957*04 
.7)657*04 
.97757*03 
.54357*02 
.40227*01 
.44767*01 
.11667*02 
.58197*03 
.28407*03 
.56627*03 
.35097*03- 
.52807*03 
.44227*00 
.52427*00 
.82297*02 
.14667*04 
.27817*04- 
.46527*03 
.85697*62 


.36957*81- 

.56217*60 

.13737*04 

.70497*03 

.■4807*03- 

.12267*03 

.58677*62 

.44857*61 

.26467*01 

.11637*01- 

.99047*03- 

.16937*03 

.63757*83 

.36297*63 


.14007*06 

.10427*05 

.66407*04 

.11627*04 

.17577*03 

.26637*03 

.46237*05 

.66207*05 

.21257*04 

•,96937*04 

•.75167*03 


• 7097E*0? 
.2033E*02 
.6?73E*04 
.69737*04 
.66767*03 
.72307*03 
.12727*03 
.10517*02 
.24727*01 
.19247*01 
.73447*04 
.74847*04 
.18927*04 
.19937*03 
.29147*02 
.91737*01 
.13947*02 
.13287*04 
.49767*04 
.31497*01 
.14367*04 
.9745E*03 
.41)27*01 
.26897*01 
.13177*00 
.29877*04 
.14177*04 
•950OE*O3 
.62657*03 
.50317*0? 
.57417*01 
.42377*01 
•  1657E*03 
.22047*04 
.66567*0? 
.69697*03 
.42397*03 


.25847*01 

.22277*00 

.78157*03 

,?640E*03 

.14847*03 

.99897*03 

.42137*83 

.28517*81 

.11867*01 

.18377*88 

.24167*83 

.28237*83 

.78837*83 

.39867*83 


Figure  75.  FLEXSTAB/LSA  Residual  Elastic  Simulator  Deck  Data 
(Continued) 


■'v'-  '&S£Sfr«&'k 


-.27356*03 

•  9326C ‘0 1 
-.65146*00 

•  1 920C  *0 1 
-.74626*02 
-.6495E.02 
-.12|4£.02 
-.3934E.02 

.343RE*03 

-.26026*00 

.19706.00 

.26516*04 

.46546.02 

■7726E.02 

-,2l25r.02 

•  1  7486  ♦0Ji 


1O2?C«03 
I8S1E*01 
5«5Mf ♦00' 
19426.01 
I377£*03 
8680E  *0? 
40836*03 
ROSbE  *0  ?■ 
tISOf .03 
?119£.00‘ 
I 0056  *00 
7065E  *03' 
|46?E«03 
5078E  *0?' 
2I40£.©3 
6814£*02. 


.7222E.01- 
.95466-01- 
.433IE*00- 
.51206.02* 
-.34046*02 
-.5208E*02 
-.23T6C«02 
-.I503E.03 
L/OELSO 
.I228E*C6- 
-.6024C»00 
.I169C.09- 
.13696*06- 
-.70R6E*07 
,579TE«08- 
■  1 3156*06- 
-.70956*07 
t/OELSI 
-•97O0E»O4 
•I510C.07- 
-.3O90E*O7 
-.56866*04 
-•2TI9C.06- 
-.I2I7C.07 
-.2340E«04 
-.2I92E.06- 
L/0GO 
•2092E.02- 
-.9044C.04 
•  1 755C.05* 
.2527E.02* 
-.3963C*03- 
.6162C.04- 
.|646C*02- 
-. 33346*03- 
L/Httl 
.6OI5C00 
••*I7«.«' 
• I509£»03 
.22466*00 
.IJ78C*82 
• 3S2SE-02 
.36O6E-0I 
.6461C*01 
•f INISHCO* 


.2003E.01' 
•I829E.OO- 
. 1 I26£»0) ■ 
,7731£*C2 
.U12E-02- 
,I263£»03 
.I286E.03- 
.2634E.0? 
15  3 

•  IC27C  *06* 
.3188C.0R 
.6II8E<«6- 
*681 1£*05' 
.14696.07 

,2145E<0R' 
.5608E.05- 
. I |40E*07 
15  3 

.I027E-04 
,6836E*0»- 

•  1 174E*07 
.I602E*04 
.6I87C.05' 
.4434C>06 
,I336E»04 
.4S33C.05' 

15  3 

.28236-01 
.34086.05 
,832CE*05' 
.  1231E»83' 
.2981E.0* 
.2127£*05’ 
,4I84E«|2' 
.20O6E.O4 
IS  3 

.I400C-02 
. I64SE*04 

•  5757E.8* 
.84246*01 
.3026E-03 
•10116*0* 
.3755C.il 
. 1 36IE*03 

•  0 


.38jTE.03-.|27lE»03 

•  13t.3E.lH-,  M67£. 01 
.Iu7ot.0l-.9194E.00 
•1533E-01  .*32kE.04 
.34S8E.U3-.5S19E.02 
. 39? 1 E .03  .278»E«03 
•57fc9£.n3-.9?3iE«0? 
•6940E.fl>-.8333E»0? 
,73h8£.fl?  ,845 1 E * 0  I • 

•  4I03£.P8-.2(.24E.DO 
.17166.01  • 1 3B4E.O0 

•  ?5b6£.n3  ,5l?s£.02' 
.4442E«07-.R99a£«Ol 
.5344£.0’-.n5«E.03 

.4S4OE.O1-.1747£.07 

.2990E.02-.920 >E.O? 
•70?6E.08  .1573£»00 
.41S6E-01  .3704E-0I 
•2574E-01  .I1RRE.04 
•4616E.03  .567hE.O? 
.67466.02  .*?35E»0l 
•2308E«OJ-.288lE*02 
.30 1 ?E «03  .1987E.03 

•  4525C  *01 

.43476.05  ,13l4£»09 
•1356E.06  .17816.06 
•2604E.08-.1558E.88 
•37R0E.05  .R329E.08 
■6343E .06  .|4»m£.U6 
•9142E.07-.7460E.07 
•2435E.05  . 31 3?E«0R 
•4916E.06 


.3M7E.03-.5495E.07 

.2O98E.06-.7R23E.O4 

.359IE.O6-.lS23C.06' 

.5I43C.03-.2I53C.O7 

.1937C.O5-.4R66E.04 

•  I l60E.06-,2o#RE.06' 
.4118E.03-.5o9nE.06 
.1417E.05 

•  I700E.01  .22I6C0S' 
•364RE.03  .28246*02 
.8710E.03-.1431E.04' 
.I259C.0I  .lr.64E.0S' 
.239RC.02  .2124C.02' 
.3134C.03-.5504C.03' 
.7974E.00  .261SE.04' 
.1625E.02 

•17S8C.00  .2607C«03 
.1 325E.02  .457RE.0O 
•4175E.02  .1747E.02 
■6540C-OI  .6672E.02 
.594  IE .01  .I469E.00 

•  5340C. 01  .12*5C»02 
. 1494C-01  .1S39C.08 
.13RlE.il 


•123RE.C3-.104fcE*02 
.2OOnE.01-.lRiOE.Oi 
•6»6OE.OO-.|20OE»O0 
.5327E.U3-.1260E.03 
•6R72E.02-.1009E.03 
•3606E.03-.9994E.02 
•4533E.03  .6023E«03 
.6I27E.O2-.8940E.02 
.R643E.OO-.6036E.OO 
.7216E.00-.II05E.01 
.1770E.0I-.O590E-0I 
.9495E.02  .6207E»03 
•2375E.02  . 10556*03 
. I  34 1 E  > 0 3  . ?698E *03 
.113TE.03-.9683E.03 
.3264E.02  ,?894E*02 
•2S46E.00  .|62IE<01 
.R394E-01  .P20E-0I 
.4l73E.03-.|395E‘C3 
. 1 847E«03  . 1 383E«03 
•676SE.02-, 1504E*O? 
•1932E.03  .3927E.03 
.799RE .02-. 42 IRE* 02 


.7733E.08-.3287E.08 
•990RE.05-.4205E.05 
•52S3E.07  .?260E»07 
.3700E.08-.1576E.00 
•740SE.05-.7156E.05 
.1C6IE.07  . T954E*06 
.910RtO7-.3R82C.07 


. I 340E*07  .4O7RE.06 
.1457E.04  .43996*03 
.I633E.06-.S074C.05 
.7597E.06  .»330E«06 
. 14RRE.04  .4562E.03 
.7956E.05-.94R9E.05 
.18506.06  .5671E*05 


.I270E.06-.1II9E.04 
•2023E.03-.I447C.OI 
. I902E.04  .6633E‘02 
•4234E.05-.S329E.03 

•  9208E  *02-. 1 104E»0I 

•  381 1 E«04  • 1 997E*02 
.9762E*04-. 1 322E»03 


.81256.04  .77086*02 
.12036*02  • 1 409C*00 
.50566*03  .12556*02 
.317RE.04  .15296*02 
.68136.01  .26546-01 
.26476*03  . 1 0966*01 
.77396*03  .17136*01 


Figure  75«  FLEXSTAB/LSA  Residual  Elastic  Simulator  Deck  Data 
(Concluded) 
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ST,\MK| 

STAMK2 

STAMK3 

CPND* 
'iXM®  H2 
NR"=  31 
NUM=  12 
NYM=  *9 
HS8=  3 
MTS  -  3 


Figure  76.  Figure  77  Precompiler  Data  (KONPACT-1) 

•••  INPUT  OAT*  CAPOS  ••• 

C  SPECIFY  PRINTING 
PRINT  OUTPUT  OAT* 

PRINT  INPUT  DATA 
C  OCEINE  VEHICLE 

SYSTEM  NO  I  fLEASTAB  CSA  A/C  FJ*  <  RESIDUAL  ELASTIC  SYMMETRIC  -  RES  I 


SLSA  DATA 

RESIDUAL  ELASTIC 

SYMM  MOOEL 

FOR  C-5A  -  CRUISE  FC  *  SIMULATOR 

•LOADS 

END 

C  NAME 
STATE 

LIST  DATA 

1 

X<  1  t 

U 

VELOCITY  ALONG  X  AXIS 

INCH/SEC 

2 

X(  71 

M 

VELOCITY  ALONG  2  AXIS 

INCH/SEC 

3 

X(  71 

0 

PITCH  BATE 

RAOIAN/SEC 

A 

X  (  A  1 

THETA 

PITCH  ATTITUDE 

RADIAN 

5 

X<  SI 

UCl 

MENDING  MODE  DISPLACEMENT 

INCH 

A 

XI  91 

UE2 

BENDING  MOOE  DISPLACEMENT 

INCH 

7 

X(  71 

UE3 

MENDING  HOOF  DISPLACEMENT 

INCH 

B 

X<  «> 

•»t  i 

BENOING  mode  displacement 

INCH 

9 

X<  91 

UES 

benoing  mode  displacement 

INCH 

It 

XII,.  I 

UE6 

BENOING  mode  DISPLACEMENT 

INCH 

11 

XIII  1 

UE7 

benoing  MODE  DISPLACEMENT 

INCH 

12 

XI1PI 

UES 

benoing  node  displacement 

INCH 

13 

xi  m 

UE9 

benoing  mode  DISPLACEMENT 

INCH 

1* 

X  1 1 A  1 

UElt 

MENDING,  mode  DISPLACEMENT 

INCH 

IS 

X(|M| 

UEII 

BENDING  MOOE  DISPLACEMENT 

INCH 

It 

XI1M 

UEI2 

MENDING  MOOE  DISPLACEMENT 

INCH 

17 

X  1 1  7  1 

UE13 

MENDING  MODE  DISPLACEMENT 

INCH 

IB 

X  1  1  A  1 

UEIA 

MENDING  MODE  DISPLACEMENT 

INCH 

19 

*1191 

UE1S 

MENOING  MODF  OISRLACEMENT 

INCH 

2# 

XI2,M 

HF1DOT 

sending  mode  ratf 

INCH/SEC 

21 

X  121  1 

UE200T 

SENDING  MODE  rate 

INCH/SEC 

22 

X  12*1 

UE3D0T 

BENDING  MODE  rate 

INCH/SEC 

23 

*1231 

UEAOOT 

MENOING  MODE  RATE 

INCH/SEC 

?* 

X  12*1 

UE500T 

MENDING  NOOF  RATE 

INCH/SEC 

25 

XI2SI 

UE900T 

MENOING  NODE  RATE 

INCH/ SEC 

29 

X  1 2*  1 

UE700T 

BENDING  MOOE  RATE 

INCH/SEC 

27 

X  1 27 1 

UEBDOT 

BENDING  MODE  RATE 

INCH/SEC 

2B 

X  (291 

UE9DOT 

BENOING  MODE  RATE 

INCH/SEC 

29 

X  (29) 

UEIOflOT 

MENOING  NODE  RATE 

INCH/SEC 

39 

X 1 3'}  1 

UEllOOT 

BENOING  MODE  RATE 

INCH/SEC 

31 

XI3II 

UEI2O0T 

BENDING  MODE  RATE 

INCH/SEC 

32 

X  I  37  I 

UE 1300T 

MENOING  NODE  RATE 

INCH/SEC 

33 

X  (371 

UEIADOT 

MENDING  MOOE  RATE 

1NCH/5EC 

3* 

-1 

OUTPUT 

XOAI 

UEISnOT 

MENDING  MODE  RATE 

INCH/SEC 

1 

R(  |  1 

SASGT 

PITCH  RATE  GYRO 

RAOIAN/SEC 

2 

R(  >1 

A2AP 

NORMAL  ACCELEROMETER 

INCH/SCC2 

3 

R  f  71 

A7EM 

NORMAL  ACCELEROMETER  ERONTS»AR 

INCH/SEC? 

A 

R<  A 1 

AZRS 

NORMAL  ACCELEROMETER  BACKSPAR 

INCH/SEC? 

s 

R I  SI 

SI 

SHEAR  EoRCE  1120.61 

LB 

9 

R I  91 

«l 

PENDING  MOMENT  1 120.91 

INCM-LB 

7 

R I  71 

T 1 

TORSION  MOMENT  (129.91 

INCH-LB 

B 

R I  71 

S2 

SHEAR  FORCE  (32S.2I 

LB 

9 

R I  91 

M2 

MENDING  MOMENT  I 32B.2I 

INCH-LB 

It 

R  I  1  7  | 

T2 

TORSION  NOMENT  I32B.2I 

INCH-LB 

11 

Rllll 

S3 

SHEAR  FORCE  (5T5.I1 

L8 

12 

RIIPI 

M3 

MENOING  NOMENT  (STS. 11 

1NCH-L8 

13 

RII7I 

T3 

TORSION  MOMENT  (57S.II 

INCH-LB 

Figure  77.  KONPACT-1  Input  Data  to  Produce  F42D  Model 
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19 

R(l4) 

S4 

SHEAR  FORCE  (ThG.OI 

LB 

15 

RUS* 

R4 

RENDING  MOHFNT  (799.01 

INCH-LB 

JO 

RIIM 

T  4 

TORSION  MOMENT  (TAG. 01 

INCH-LB 

IT 

RUT) 

S5 

SHEAR  FORCE  (920.01 

LB 

IB 

RIIRI 

RS 

RENDING  MOMENT  (920.01 

INCH-LB 

19 

RU9) 

T5 

TORSION  MOMENT  (929.0* 

INCH-LB 

INPUT 

1 

U(  \) 

poail 

aileron  deflection 

RADIAN 

2 

U(  ?) 

ROE  I 

INBOARD  ELEVATOR  OErLECTjON 

RADIAN 

3 

U<  H 

ROEO 

outboard  elevator  deflection 

RADIAN 

A 

U(  4) 

R0AIL00T  aileron  deflection  rate 

RAOI AN/SEC 

s 

U<  SI 

ROE  1  dot 

inboard  elevator  deflection  rate 

RAOt AN/SEC 

6 

U<  M 

RDCOOOT 

outboard  ELEVATOR  DEFLECTION  rate 

RAOIAN/SEC 

7 

U(  71 

NGI 

GUST  INPUT  AT  -1020  IN  FROM  CG 

INCH/SEC 

S 

u<  m 

MG2 

GUST  INPUT  AT  9  IN  FpON  CG 

1NCM/SEC 

9 

U(  91 

VG3 

GUST  INPUT  AT  1920  IN  FROM  CG 

INCH/SEC 

10 

UUM 

MG| DOT 

GUST  INPUT  RATE 

INCH/SEC2 

II 

Udll 

KG200T 

GUST  INPUT  RATE 

INCH/SEC2 

12 

-1 

ENO 

U(|?l 

MG300T 

GUST  INPUT  RATE 

INCH/SEC2 

C  DEFINE  CONVERTED  VEHICLE 

SYSTEM  NO  I  CONVERTED  FLFXSTA8  C5A  a/C  F32  (  FLEXSTAR-RES  TO  Hl/GELAC  » 
•CONDITIONING  DATA 
C  SCALING  DATA 
SCALE  THE  VARIABLES 


X(3I 

. I697B9E  99 

RAOIAN/SEC 

INCH/SEC 

U(T» 

.837333E-0I 

INCH/SEC 

FEFT/SEC 

U(BI 

.833J33E-BI 

INCh/SEC 

FECT/SEC 

U  (91 

•H33333E-0I 

INCh/SEC 

FEET/SEC 

Udtl 

.831333E-0I 

INCh/SEC 

FEFT/SEC 

U(!!> 

.83’333E-01 

INCh/SEC 

feet/sec 

U<  121 

.8V3333E-0I 

INCh/SEC 

feet/sec 

R  (21 

.2S88B0F-02 

INCH/SEC2 

IG 

R(3I 

•2S8880E-02 

INCM/SEC2 

IG 

RIAI 

.2S8809E-02 

INCH/SEC2 

IG 

ENO 

C  RESPONSE  SPECIFICATIONS 
SELECT  OUTPUTS 
R(H-R(|9I,*(2).«(3I. 

END 

C  REDUCTION  AND  SHUFFLING  OATA 
RETAIN  STATFS 

X(2l.X(3:  *X|2U-XI3AI  .A (51 -A <191. 

ENO 

C  DEFINE  ACTUATOR 

STSTEN  NO  2  ACTUATOR  (  FIRST  ORDER  ) 

•TRANSFER  Function  OATA 
BLOCK  1 

I  2  .  10000'E  01  2  1  •  1 6A66TE  00  2  »  .IOOOOOE  01 
-l 

BLOCK  2 

I  2  >10000  'E  01  2  1  .J33333E  00  2  2  .  IOOOOOE  01 

-I 

BLOCK  3 

I  2  .10009  >F  01  2  l  , 1 31333E  90  2  >  .IOOOOOE  01 
*1 
ENO 

C  CONNECTION  OATA 
UI/U 

I  I  ,10000'if  01  2  2  .IOOOOOE  01  3  1  .IOOOOOE  01 
-I 

UI/RI 

-I 

R/RI 


Figure  77.  KONPACT-1  Input  Data  to  Produce  F42D  Model 
(Contini/  d) 
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I  (  .lOOOOlE  01  2  >  .I0OOO0F  01  3  1  .100000F  01 
A  I-.600000E  •»  5  2-.75C000E  01  0,  3-.750000E  01 
-1 
R/U 

0  1  .600000E  01  5  ?  .750000E  01  6  3  .7SOOOOE  01 
-1 
ENO 

C  NAME  LIST  OAT A 
STATE 


1 

x< 

KA 

actoator  state 

RADIAN 

? 

x< 

«El 

actoatop  state 

RADIAN 

3 

x< 

XEO 

actuator  state 

RADIAN 

•  1 

OUT»UT 

1 

R( 

n 

DEL  A 

aileron  position 

RADIAN 

I 

R< 

DELE  I 

INBOARD  ELEVATOR  POSITION 

RADIAN 

3 

R( 

i) 

DELEO 

OUTBOARD  ELEVATOR  POSITION 

RADIAN 

A 

R  C 

4  » 

OELADOT 

AILERON  VELOCITY 

RADIAN/SCC 

5 

R  < 

oeleiootinboaro  elevator  VELOCITY 

RADIAN/SCC 

6 

-1 

INPOT 

R  < 

M 

deleodotoutboard  elevator  VELOCITY 

NAOIAN/SEC 

1 

U< 

)  t 

OOELA 

AILERON  CONTROL  INPUT 

RADIAN 

2 

U< 

HDELEI 

INBOARD  ELEVATOR  CONTROL  IN»UT 

RADIAN 

3 

U1 

UOELEO 

OUTBOARD  ELEVATOR  CONTROL  INPUT 

RADIAN 

-1 

ETO 

C  OEEINE  COST  HOOEL 

SYSTEM  MO  .1  GOST  MODEL  (  WITH  FIRST  ORDER  KUSSNER  I 

(QUADRUPLE  DATA 
X00T/X  7  7 


»  i- 

.IBA93’£ 

02 

1 

7 

. 18a93AE 

02 

? 

2-.I3A270E 

02 

2 

5  • 1 3A270E 

02 

3 

1-.800330E 

01 

3  A 

.10000'IF 

01 

A 

1 

'22AlftfiE 

03 

4 

3-.90C7A0E 

0? 

4 

A-.160070E 

02 

5 

S-.9I555AE 

01 

5  7 
_  1 

.91555  IE 

01 

5 

6- 

.A2O0OOE 

00 

6 

7-.AAIO00E- 

■01 

7 

6  • 1O0OO0E 

01 

XROT/U 

7 

1 

6  1- 
-1 

R/X 

.237IOnE 

00 

7 

1 

.79370OE 

00 

6 

7 

1  1 

.lOOOOlE 

01 

2 

2 

•1O00OOE 

01 

3 

3  .1AOOOOE 

01 

4 

1-.IBA93CE 

02 

A 

7  .18A930E 

02 

5  2- 

.1 3527 sE 

02 

5 

5 

.13A270E 

02 

6 

3-.0OO330E 

01 

4 

A  » 1 OOOOOE 

01 

-I 

R/U  6  I 

-1 
END 

C  NAME  LIST  DATA 
STATE 


1 

X( 

01 

KUSSNER  STATE  (  NT  1 

FEET/5EC 

2 

X( 

>\ 

R? 

TRANSPORT  0EL*V  STATE  1 

w  > 

►EET/SCC 

3 

X  ( 

03 

TRANSPORT  OELAT  6TATE  ( 

T  ) 

keet/sec 

A 

X( 

41 

PA 

TRANSPORT  DELAY  STATE  ( 

T  ) 

5 

x< 

05 

KUSSNER  STATE  t  V 

1 

6 

x< 

4) 

06 

VINO  FILTER  STATE 

7 

-1 

OUTPUT 

x< 

7) 

ur 

VINO  DUST  STATE 

1 

R  i 

l  \ 

VGN 

VINO  DUST  TO  NOSE 

FEET/SEC 

2 

R  < 

7) 

WDM 

VINO  DUST  TO  VINO 

FEET/SEC 

3 

«< 

"»> 

vr.T 

VINO  RUST  TO  TAIL 

FEET/SEC 

A 

R  ( 

4) 

VONDOT 

VINO  r.UST  RATE  TO 

NOSE 

FEET/SEC2 

5 

R  ( 

VRVOnT 

vino  dust  rate  to 

VI  NR 

FEET/SEC2 

0 

-l 

INPUT 

R< 

M 

WGTDDT 

vind  rust  rate  to 

mil 

FEET/SEC2 

1 

U< 

«  > 

ETXO 

unite  noise  ;npjt 

to  gust  model 

FEET/SEC 

-I 


Figure  77.  KONPACT-1  Input  Data  to  Produce  F42D  Model 
(Continued) 


FND 

C  DEFINE  PL  1*  r 

SrSTEM  NO  4  PLAmT-fa*  (CiUU/E-UED  FL£xSTA«  CS«  A/C  < ►» L «  *  »  ACTUATOR  •  OUST  MODEL! 

ilNTEPCONNEfTlON  DATA 

Ull/u 

-I 

OU/RI2 

i  i  ,10000  t  fi  ?  t  . n  oii'h-:  iii  i  1  , iooooof  oi 
-1 

un/nn 

7  3  ,10000  F  01  8  ?  .lOl.liOfiE  01  a  1  .IOOOOOC  01 

-1 

U12/U 

1  2  ,10000  E  01  ’  3  ,10f>(lrt£  III  3  4  ,1100(1'./  01 
-1 

UI3/U 

1  1  .10000  -F  01 
-l 

P/P  11 

1  1  ,10000  F  01  2  2  .  |  'I n  0  0  0  E  Cl  3  1  ,  lOOOOOE  01  4  4  , 1 OOOOOE  01  5  8  ,  l  OOCOOE  01 

6  8  ,10000  F  01  7  7  , 1  v.OOOE  (I 1  «  A  .100000E  01  0  P  .10C000E  011010  .100900E  01 

111)  ,10000-  F  01121*  .liM'/dOE  011311  .HOOOJf.  0 1 1 4 1  *•  •  1 COOOOE  U1151S  .100001E  01 

1616  ,10000' F  011717  .IC'ilOOf  Oil'll*  .lOOOOOE  ililPlft  .100000E  012020  .1I337HE-03 

2121  .60660- F-032S20  ,IO'"iOOE  01 
-1 

ft/P  1 2 

22  1  .10000  F  0123  *  .U"'OOOE  0I.?4  3  .100DC0E  U1 

-1 

P/P  I  3 

20  2  , 1  1  337AE-C 3 

-1 

END 

C  NAME  LIST  DATA 


OUTPUT 

24  R12I 

ALPHA 

ANGLt' 

OF  ATTACK 

RADIAN 

21  R(2'l 

,  | 

OCG 

PITCH 

RATE  AT  CG 

MADIAN/SEC 

END 

C  OEFINf  PL  ".NT 

DESIGN  MODEL 

SYSTEM  NO  4  PLANT-F420  (  COnvFWTED  FLf  AST  A®  CSA  A/C  ♦  ACTUATOR  •  GUST  MODEL  1 
SC0N0ITI0N1NG  OATA 
C  NO  SCALING  OATA 
FNO 

C  RESPONSE  SPECIFICATION.; 

SELECT  CONTROL  INPUTS 
U12I-U14I . 

SELECT  GUST  INPUTS 
U(1 1  . 

CONSTRUCT  O.SIGN  PESDONSf S 

R(6),R(7l.fi  (PI, ft(101.P(l»l«P(13l.P<ISI.P  1161. PURI. PllPI. 

ROOT  161  .ROOT  (71  .ROOT  1PI  .MOOT  (10)  .BOOT  11?)  .BOOT  113)  .ROOT  115)  ,«OOT  ( 16)  « 

ROOT  (181  .ROOT  UP)  .R1201-M2S1  .R00T12SI . 

SELECT  SENSOR  OUTPUTS 
RI31.R(4I.R(2).R(t>, 

END 

C  NO  REDUCTION  AND  SHUFFLING  DATA 

END 

STO“ 


Figure  77.  KONPACT-1  Input  Data  to  Produce  F42D  Model 
(Concluded) 
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it 


1 


I 


I 


tim. 


CONt)* 

R2 

NRH  =  31 
NUH -  12 


Figure  78.  Figure  79  Precompiler  Data  (KONPACT-1) 


-4 ' 


I 

It 


1 


1 


input  d^ts  cards  •  ••> 


PRINT  INPUT  DAT* 

print  output  data 
CONTINUATION  run 

C  OtriNE  F14  plant  MODE!  -  RESIDUAL  IZINC,  STATES  .  RESPONSE  AND  MEASUREMENTS 
ST  STEM  NO  a  PLANT-  FJ4PR  (F|_C*STAi  CSA  A/C  *  ACTUATOR  •  OUST  MODEL  I 
SCOND I T ION IN6  DATA 
C  NO  SCALING  DATA 

end 

C  NO  RESPONSE  SPEC  I E I  CAT  J  ON  DATA 
END 

C  REDUCTION  AND  SHUFFLING  DATA 
RETAIN  STATES 

*  (I I-X18T  .X  (181-X (23)  .X  (33)  -X  (4£l  . 

RESIOUALIZE  states 
A  I9I-X  1 1  7>  t « (241-x  < 32.1  . 

RESIOUALIZE  STATES  IN  OUTPUTS 
ROl-POtl  . 

ENO 

REFERENCE 

SYSTEM  NO  <t  PLANT-M2D  <  CONVERTED  ECEXSTAg  CSa  A/C  .  ACTUATOR  ♦  GUST  MOUEL  ) 
ENO 

C  DEFINE  F24  PLANT  MODEL  -  RESIOUALlZING  STATES  AND  RESPONSES 
C  AND  TRUNCATING  MEASUREMENTS 

SYSTEM  NO  4  PLANT-  F £4R T  (FLEXSTA*  CSA  A/C  »  ACTUATOR  •  GUST  MODEL) 

sconoitionino  data 
C  NO  scaling  data 
ENO 

C  NO  RESPONSE  SPECIFICATION  DATA 
ENO 

C  REOUCT  I  ON  AND  SHUFFLING  DATA 

retain  STATES 

«  (I I -XI* I .XI |»>-A (231 .*( JJI-X <4E) . 

RESIOUALIZE  STATES 
X  (9>-X  (171 .« (241-x (32) . 

RESIOUALIZE  STATES  IN  OUTPUTS 
RIJ)-R(ZT) . 

TRUNCATE  states  IN  OUTPUTS 
R(ZRI-MH)  . 

END 

REFERENCE 

SYSTEM  NO  4  PLANT-FA2D  (  CONVERTED  FLE*STAB  C5*  A/C  ♦  ACTUATOR  .  GUST  MODEL  ) 
ENO 

C  DEFINE  Fas  PLANT  MODEL  -  TRUNCATING  STATES  •  RESPONSE  AND  MEASUREMENTS 
STSTEM  NO  4  PLANT-  F24TT  (FlEXSTAB  CSA  */C  •  ACTUATOR  •  GUST  MODEL) 

ICONOIT IONIhG  DATA 
C  NO  SCALING  DATA 
FNO 

C  NO  RESPONSE  SPECIFICATION  DATA 
f  HO 

C  REDUCTION  ANO  SHUFFLING  DATA 
RETAIN  STATES 

*ni-«(8i.xii#)-«izii.xmi.xi*z), 

END 

STOP 


Figure  79.  KONPACT-1  Input  Data  to  Produce  F24RR, 
F24RT,and  F24TT  Models 
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CjtlifcTatv  r  f  ■  ■ ; fc/*’  'j  rY.1.  |,i.  <r- 1 4S-~“ 


5TAHK3 
STAHK4 
NXM  =  28 
NRH=  31 
NOH=  9 
NfM=  47 
N5B=  2 


Figure  80.  Figure  81  Precompiler  Data  (KONPACT-1) 


I 


I 


•  input  mat  a  cards 


PRINT  INPUT  DATA 
print  OUTPU'  OaU 
CONTINUATION  RUN 

C  DEFINE  CONTROLLER  FOR  CSA  ALDCS  DESIGN 
SYSTEM  NO  5  CONTROLLER  FOR  C5a  ALDCS  DESIGN 
^SIMULATION  DATA 

C  SET  which  CONTROLLER  SWITCHES  SHOULD  BE  DM 

CONTROLLER  switches  ON 

ALDCS 

SAS 

H'.Cl 

MLC2 

END 

C  READ  ANY  CHANGE  IN  CONTROLLER  CAINS 
CONTROLLER  GAIN  VALUES 
FAF 
9  9 
WMI 

-1403. 7B 
KO 
0.0 
KNF 

o.o 

K  P 
0.0 
ENO 
END 

C  NAME  LIST  DATA 
STATE 


1 

*  ( 

1 1 

p 

PILOT  FILIEP 

2 

X  < 

A  P 1  HL 

L AGGEo  NORMAL  ACCELERATION 

1C 

3 

X( 

M 

MLC 1 

FULL  ST^TE  MLC  FOR  AILERON 

4 

X  ( 

u) 

C.LAF 

GUST  LOAO  ALLEVIATION  filter 

5 

X  { 

m 

F3E 

3RD  RfNOING  MOOE  FILTER  ON  ELEVATOR 

6 

X  ( 

*1 

HP 

high  pass  filter  on  elevator 

7 

X  ( 

n 

MLC? 

FULL  STATE  MLC  for  ELEVATOR 

-  * 

OUTPUT 

i 

p( 

I > 

HOELA 

AILERON  COHmawd  OUTPUT 

RADIAN 

2 

k  ( 

ITDELFI 

INBOARD  ELEVATOR  COMMAND  OUTPUT 

radian 

3 

R( 

i  \ 

IIDELFO 

OUTBOARD  ELEVATOR  COMMAND  OUTPUT 

RADIAN 

1 

INPUT 

i 

U( 

\ > 

U1 

AILERON  OPTIMAL  CONTROL  INPUT 

2 

U< 

■>» 

U2 

INBOARD  ELEV  OPTIMAL  CONTROL  INRUT 

3 

U( 

11 

ETAP 

MMITE  NOISE  INPUT  TO  PILOT  FILTER 

4 

U< 

U) 

DELA 

AILERON  POSITION 

Radian 

5 

U< 

*,) 

OELEI 

INBOARD  ELEVATOR  POSITION 

radian 

6 

U( 

M 

A21R 

ACCELFROMETEJJ  OUTPUT  (21) 

1G 

7 

U< 

7) 

AFUS 

FUSELAGE  ACCELEROMETER  OUTPUT 

1G 

S 

U( 

R» 

TEUS 

PITCH  RATE  GYRO  OUTPUT 

HADIAN/SEC 

9 

U( 

91 

ACC. 

ACCELERATION  AT  cg 

INCM/SEC2 

-I 

END 

C  DEFINE  REMICEO  CONTROLLER  FOR  CSA  ALDCS  DESIGN  REPFAT 
SYSTEM  NO  5  CONTROLLER  FOR  CSA  ALOCv.  OESIGN  IREDUCEDI 
SCONDITIONING  DATA 
C  NO  SCALINO  DATA 
ENO 

Figure  81.  KONPACT-1  Input  Data  to  Produce  F24TT  Plus 
Controller  Model 


I 
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*«%• 


J 


C  NO  3fSt>ON-,r  SPFCIFtCATtON 
FND 

C  REDUCTION  DATA 
RETAIN  ST  a TrS 
XID-XIA). 

RESIOUALIZE  states 
XIS)  . 

RESIDUAUZF  STATFS  Hi  OUTPUTS 
R(1I-R(3I . 

END 

REFERENCE 

SYSTEM  NO  A  PLANT-  F?AwR  IFLEXSTAs  CSA  a/C  »  ACTUATOR  •  OUST  MODEL I 
FND 

c  define  overall  system 

system  NO  6  OVERALL  SYSTEM  II  epi.oc  ♦  RFOJCED  CONTROLLER! 

INTERCONNECTION  DATA 

II14/IJ 

A  A  .10000  E  01 
-1 

UI5/IJ 


1  1 

.10000  E 

ul 

?  ? 

•  lOniHUJr 

•n  i  1 

. 1  no  0  0  of 

0) 

■  I 

UIA/RIS 

1  1 

.10000  E 

01 

->  ? 

•  1  MOOOF. 

01  1  1 

. ) Ooooof 

01 

1 

U15/QI4 

A23 

.10030  F 

01 

. 1 o^OOnr 

Of  f-2'i 

.1000006 

ul  730 

•lOCOOOfc 

01  S3 1  •  luOOOOE 

01 

927- 

.10000  F 

01 

<*?/> 

.M^OOpF 

04 

R/RI4 

1  1 

.10000  F 

01 

7  7 

.  1  wf’UUOP 

in  i  ) 

. looooor 

01  4  A 

, 100000E 

01  5  4  •  1000 OOF 

01 

6  6 

.10000  F 

01 

7  7 

,U«-',n,nF 

*  !* 

.ioooeor 

01  9  9 

.IQDGOOL 

011013  .100000F 

01 

mi 

.10000  F 

nn?iR 

.|.J  .OO  lF 

i-tnn 

. i nnoooF 

01  1  A  l  A 

.  1  ii  0  COOL 

on  si  s  •  iooooof 

01 

1616 

.  1  3000  'F 

011717 

•  K-OOdF 

i'U^H 

•  1 0  0  0  0  0  E 

01 1919 

.1.130002 

012020  .IOOOOOF 

01 

2121 

,10000  E 

01???? 

•  1  Ii*i4HQF 

.lonooOE 

0 1 2  4?9 

. 1030902 

012530  . IOoOOOE 

01 

2631 

.10000  F 

>>1 

272? 

.11485  IE 

03?7?3- 

.2297O0F 

J32VS- 

,4  7 1 200E 

0? 

2622 

•13?«S' f 

.601tfO')F 

03^2  4- 

.2667006 

0  6  2  4  2  S  - 

,61 ?000t 

07 

-l 

P/RIS 

27  3  .2765?  F  0  3?.^.J«  .3?*-90'>E  OA 


29  1  .10000  F  o I  Jo  ?  .lO'  OOOE  01 

-I 

End 

C  name  list  DATA 
DUT0UT 

?7  «<?M  FwllOT  IMP  M)IHL  ERROR  pate  FOR  a 

?8  m<?mI  EOIIOT  IMP  MODEL  ERROR  RATE  1  OR  0 

?9  Wi?i>  ill  A ILEPDN  OPTIMAL  CONTROL  INPUT 

3.1  RI3  I  U?  INH0A90  ELFV  OPTIMAL  CONTROL  IN»UT 

-t 

FND 

REFERENCE 

SYSTEM  NO  A  PLANT-  fua&i  (FI.EXSTtp  CSA  A/C  <  ACTUATOR  ♦  DUST  MODEL) 
FND 

C  OEPINE  OVtRALL  SYSTEM 

SYSTEM  NO  6  OVERALL  StSTfm  (I  F?4RT  .  REDUCED  C3NTROILE9I 
S INTERCONNECT  ION  DATA 

0]  A/O 

A  4  .10000  E  01 
-I 

UI5/U 

1  I  .10000  F  SI  ?  ?  .  IU'00»F  01  i  *i  .  lCOOOIIF  01 
-1 


.iroooop  oi 


Figure  81.  KONPACT-1  Input  Data  to  Produce  F24TT  Plus 
Controller  Model  (Continued) 


uia/pia 


t  l 

,  1  loot,  F 

1  2  2 

.  )  -j 

<M  I  ^  .i^000?F 

1 1 

”  1 

UIS/9I* 

,1100  F 

r  1  S2A 

•  <* 

u  ^2^  .icr*oonF 

Ul  730  •  1 OOOOOE 

01  831 

.100000E 

01 

Q?7- 

.1000  F 

Cl  1?2 

St'* 

~  1 

»/i»h 

\  1 

. luoor  t 

01  2  2 

.  1  F.^OlML 

'M  t  5  .  1  nrwrju*-: 

ul  A  A  .lOOUOUt 

01  5  S 

.1000002 

01 

b  b 

.ii/oo:  f. 

0  1  7  7 

J 1  M  fi  . lftOOOGt 

ii  »  )  .iioooae 

011010 

•  1 0000  IE 

01 

nil 

.  1  fUltO  F 

L 11 2 | 2 

omn  •  i oo«»oor 

01 1  a  1 A  .  1000COE 

0I151S 

. 1 OOOOOE 

01 

t  6 1  f> 

.1100.1  F 

11  717 

imhi-1  .innootiF 

Ul 1010  .1000001 

012020 

• 100000E 

01 

?\?\ 

.101/10  F 

0  1  2  2  2  2 

.  1 vnOOOr 

. l nnoooF 

012420  .lOCOOOt 

012530 

•  1 OOOOOE 

01 

Pf5.11 

.10100  E 

.-.1 

?7?2 

.llASSIf 

.132  721- 

W--.^P^700f 

012725-.M7I200E 

02 

?*?? 

.  I  VMS--F 

5A3A/1- 

u3?  »?*-. 2^S?r»UF 

04222S-.6I700UE 

03 

“  t 

*/»!  3 

P7  3 

.27182  F 

'.120  3 

,  3^9011? 

0* 

-  \ 

a/u 

?9  1 

.10000  F 

'll  Vi  2 

.  1  O^rtOf'F 

i>l 

-i 

F'in 

C  NAME  U^T  nATA 
ntJt°UT 

2  7  8(271  F,.bf/T  [MP  P'/liFi.  F»8./A  PATF  F0'(  8 

28  P(?..|  r,)ilOT  IMP  Mil]  L  FPPJt/  Mtlr  FOP  0 

?9  P(2'U  'll  AILERON  (IWIM4L  C0NT8CM.  [NPJT 

3,  M3  I  M2  INHOAii/  FLt'w  OPTIMAL  CONTROL  IN2(M 

-I 

FNO 

R£FE8fmCF 

SYSTEM  MO  A  PL  A, IT  -  fjlTT  ( FlExST  Am  FSa  A/C  ♦  ACTUATOR  *  OUST  HU1EL> 
FMO 

C  DEFINE  0V--PALL  SYSTEM 

SYSTEM  MO  <5  DVE8AU  SYST/M  (I  F24TT  •  kfO./CED  COMTMOllFP) 
MNTEBOHNEi'TIOU  OAT  A 
0 1  A  /II 

A  A  .11/0,13  f  1 

-1 

UIS/IJ 

I  1  . 10000  .F  11  2  ?  .li'.i’O'F  III  V  1  .l.iO'IOLF  ul 

-l 

UlA/PIS 

I  1  ,  1  0091)  F  II  2  2  . 1 U  uloOf  01  3  3  .lOOOOl’F  Ul 
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UIS/SIa 
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, 1 OOOu 
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,88’u  'I  2 

1 A 
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1  1 
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.1  /Mil  It 

11  3  3  . 1  iiilO.lE 

01  A  a  ,  1 50000E 
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•  1 OOOOOE 

01 

b  ft 

,11000 

F 

jl  7  7 

.to  Ifo-IE 

01  A  A  .l"3001'r 

v 1  09  .150000L 

01 1010 

. 1 OOOOOE 

01 

mi 

,1100'! 
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:n?H 

.11  10'  F 
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w 1 1 A l A  ,  1  OQOuOE 

01151ft 

. 1 OOOOOE 

01 
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OH  717 
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01181-  .110000" 
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012020 

. 1 OOOOOE 

01 

2121 

,11001 
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U12A20  •  1 OOOOOE 

012530 
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Figure  81.  KONPACT-1  Input  Data  to  Produce  F24TT  Plus 
Controller  Model  (Continued) 
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-I 

PND 

C  NAME  LIST  DATA 

output 

27  X<z?>  cMOOT  IMP  model  ERROR  PATE  POP  V 

2*  P(2«l  EOOOT  IMP  MOO  PL  ERROR  RATE  FOR  0 

29  R<2«>  til  AILERON  OPTIMAL  CONTROL  INPUT 

Pilot  U2  INSOARO  ELEV  OPTIMAL  CONTROL  INPUT 

ENO 

STOP 


K 


Figure  81,  KONPACT-1  Input  Data  to  Produce  F24TT  Plus 
Controller  Model  (Concluded) 
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CONDK 
NAM-  11 
NNM=  '0 
NUM=  4 


Figure  82.  Figure  83  Precompiler  Data  (KONPACT-1) 


(KONPACT-1  Output  is  Shown  in  Figure  90) 

•••  INPUT  DaTA  cards  ••• 

P^INT  INPUT  PAT* 

PRINT  OUTPUT  OATA 
CONTINUATION  RUN 
((TTCRtNCC 

SYSTEM  NO  6  OVERALL  SYSTEM  !<  ER..3H  ♦  REDUCED  CONTROLLER* 

END 

C  DEFINE  OVERALL  SYSTEM  DESIGN  MODEL 

SYSTEM  NO  6  OVERALL  SYSTEM  (<  ERROR  ♦  REDUCED  CONTROLLER*  DESIGN  MODEL 
SCONOITIONInG  DATA 
C  NO  SCALING  DATA 
ENO 

C  RESPONSE  SPEC1EICATI0NS 
SELECT  CONTROL  INPUTS 

0  O  I  .U (?)  a 

SELECT  GUST  INPUTS 
U(*l  .JO*  . 

CONS'TUCT  DESIGN  RESPONSES 

S  (?7I  .01 1 »  ««(?>  .ROM  .AIT*  .R<*»  aAr.OTIlS*  .R(S*  .RIG*  .XOOT  Oft*  •*< » XI  «R  <8*  •*  ( IS)  • 
R<<n.RU0t.<(IM. Rill). R<l?i«*ni. RIO*. Rmi.X<«I.R*15).RtlNI.A<Sl. ROT). RUB*. 

X  (ft)  .RU9I  . I XCO  .*(  El  .A  <81  .HIRE)  .3  (RBI  ,S<2|»  .R<?9*  .ROO)  . 

SELECT  SENS  IP  OUTPUTS 

X  os*  .X  081  -X  <?*)  .Ml  I  .X(OI  .X  os*  -XliEI  .A  <J*-«  (S)  aR| RAI  aX<ftl-X  (Ml  .Ml  25*  . 

R!?6I  .X  OOI  .XO?l  .XOAI  .X  <2ftl-X<2A>  , 

END 

C  REDUCTION  ANO  SHUFFLING  DATA 
RETAIN  STAT*  S 

X  (I  I -I  0  6) ,  «  <26* -*0  281  .X  OX*  .MRS*  .«  OR* -A  <R*l  . 

ENO 

REFERENCE 

STSTEM  NO  ft  OVERALL  SYSTEM  O  f»*RT  .  REDUCED  CONTROLLER* 

END 

C  DEFINE  OVERALL  SYSTEM  nESlGN  “ODEL 

SYSTEM  NO  6  OVERALL  SYSTEM  O  ERAST  •  REOUCEO  CONTROLLER*  DESIGN  MODEL 
SCONOITIONInG  DATA 
C  NO  SCALING  DATA 
ENO 

C  RESPONSE  ePECIEICATIONS 
SELECT  CONT  >0L  INPUTS 
UO  I  .U<?l  . 

SELECT  GUST  INPUTS 
UIAI .01 1* . 

construct  o.sign  responses 

XOTI.Rdl.MRI.HOM.RniaPIAl.xnOlOSt.-'tSiaRIGI  aAODT  0P*»M«El.R(8».X05). 

R<9*  .ROO*  .1(161  .RO  1 1  .RO?)  illll  .RO  J>  .HO*t  .*1*1  «ROS)  .ROfc*  .XIS1  .RTIEI.R0  8)  • 

X  (ft*  ,RI  I  R*  .SIRR)  .A*T*  .AIM*  ,UI?E*  .RIR8I  .  R  <  21  I  .R(RR)  .Wf  TO)  . 

SELECT  SENSOR  OUTPUTS 

XIZ5I  .XORI-XI2AI  .XII I  .AIR)  .AOS)  aXIIE)  .X  Ill-X  (SI  ,R  ( RAI  .  X  ( 6 1 -X 1 8 1  .RI25*  . 

R<2ft*  »X  O0«  .111X1  .XII  A)  .  <  |?M-X(?*t . 

ENO 

C  REDUCTION  ANO  SHUEEL 1NG  OATA 
RETAIN  STATES 

XO  t-XUGI  .*  (?8*-X(R8)  .X  OX)  aX  (  >5*  «>  08*-*  (RA)  a 
END 

REFERENCE 

SYSTEM  NO  ft  OVEMALL  SYST*m  (<  EratT  •  RE'^JCEO  CONTROLLER* 

ENO 

C  DEFINE  OVERALL  SYSTEM  DESIGN  MODEL 

SYSTEM  NO  ft  OVERALL  SYSTEM  O  ER*TT  .  HFOUCiO  CONTROLLER)  DESIGN  MODEL 
SCONOITIONING  DATA 

Figure  83.  KONPACT-1  Input  Data  to  Produce  F24RR  Plus 
Reduced  Controller  Model 
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C  NO  SCALING  data 
END 

C  RESPONSE  SPECIFICATIONS 
SELECT  CONTROL  INPUTS 
U<|1 .U(2I . 

SELECT  GUST  INPUTS 

uui  tunt . 

CONSTRUCT  DESIGN  RESPONSES 

R(2T)»R(|)«R(2)*RI2fc>«RU>*R<A)«ADOT<l5)»R(5l»R(6) .AOOT (161 »R(T) »R(A>  « A (IS) • 

R(RI  .R( 10). XI 161  .R(ll)  .fill?)  .AO)  ,Rll J)  .RUAI  .A  (A)  .Rtl5).R(l»>  .A  (5)  .RtITt  .R(l«>. 
A(6I •R(I«I *R (20) ,A(7I .AIR) .«(27l .9(28) .R(21) .R(2R) .ROO) . 

SELECT  SENSOR  OUTPUTS 

A(25).X(IAI-X<2*).AUI.»(RI.XU5)-A(l7l.xO|.X(5I.R<JA).At6l-A(*).R(25l. 
RI26UAII0)  .A  (12)  •A(1*I.A(2GI-A(28)  . 

ENO 

C  REDUCTION  AND  SMUFFL INO  OATA 
RETAIN  STATES 

A<1)*A(I6I.  A(26)*A<28)*A(|7)«AI2Sl«A!i8)-A(2AI. 

END 

STOP 


Figure  83.  KONPACT-1  Input  Data  to  Produce  F24RR  Plus 
Reduced  Controller  Model  (Concluded) 
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C  DESIGN  USING  DU*  FOR  THE  OfMiyisTfiAT  ION  E»«MPLE 
C  READ  FOR  WhAT  PROGRAM  I  OU*«FFoC,LS»  I  THE  DATA  IG 
SOU*  OATA 

C  READ  IF  OATA  is  ON  CARDS  ONLT  0*  ON  CARDS  anu  TAPE 
OATA  ON  CAR'lS  AND  TAPE 

C  IF  OATA  Is  ON  CAROS  ANn  TiPr  R£S|)  THt  LAhEL  TO  OBTAIN  OATA  ON  TAPE 
SYSTEM  NO  6  OVERALL  SYSTEM  <1  F£4pR  ♦  RtOUCfO  CONTROLLER*  DESIGN  MODEL 
C  REAO  DATE  AND  USER  10 
JAN  10*  TS  J  *  MAHFSH 
C  NOP  -  NO  nr  VaRURLFS  cFIng  PLOTIf u 
C  READ  NOP 
9 

C  GO  TO  100  IF  NOP.fO.O 

C  REAO  (PLRl  1 1  »ITITL  1 1 1  .UNIT  (II  •  YMINU  )  •YMAAM  I  .SCAL  I 1  >  .1*1. NOP  I 
C  REAO  T.DT.ST.TI.T2 
C  100  CONTINUE 
C  REAO  IMAX.ITER.ITFPO 
?»1S  0 

C  NOCOV.I  NO  COVARIANCE  ANALYSIS 
C  NOCOV*?  COVARIANCE  ANALYSIS 
C  N0C0V*3  SNIP  CORRELATION  ANALYSIS 
C  NSTER*0  NO  STEP  INPUTS 
C  NSTE°*I  STEP  COMMANDS 
C  NSTEP*?  STEP  GUSTS 
t  NSTER*1  BOTH  (1  ANO  21 

*  NSTEP*<*  NO  STEP  INPUTS  -  TRANSIENTS  ONLT 
C  NSANO*0  NO  PANOOM  INPUTS 
C  NRAND* 1  GUSTS 

C  NPRIN*0  00  NOT  PRINT  RESPONSES 
C  NPRIN*I  PRINT  RESPONSES 
C  N»L0T*i)  NO  PLOTS 

c  n»lot*i  calcomp  plots 
C  NPLOT*'  Line  PRINTER  PLOTS 
C  N»L0T.1  BOTH  (I  ANO  ?> 

C  READ  NOCOv.NSTCP.NRANO.NPRIN.NPlOT 

3(000 

C  IMPK«I  NEW  INPUT  GAINS 
C  INPW  NEW  STARTING  ROUTINE  CAINS 
C  INP*«3  USE  GAINS  IN  STORAGE 
C  IMP**fc  USE  INPUT  GAINS  IN  STORAGE 
C  READ  INP* 

I 

C  NCONT  *0  OONOT  CON»UTE  OPTINAL  GAINS  -  USE  INPUT  GAT  NS  ANO  DATA  IN 
C  COVARIANCE  AND  TINC  RESPONSE  ANALYSIS  ONLY 

C  NCONT*!  COMPUTE  OPTIMAL  GAINS 

C  NCONT*?  COMPUTE  OPTIMAL  GAINS  NTTM  AUTOMATIC  ft  SELECTION  ON  CONTROL  PAT 
C  REAO  NCONT 
I 

C  REAO  FLIGHT  CONOITION  NUMBER 

*12)01 

f,  NX  -  NO  OF  STATES 
C  NR  -  NO  OF  RESPONSES 
C  NU  -  NO  or  CONTPQI  INPUTS 
C  NN  -  NO  OF  OISTUNNANCE  INPUTS 
C  NF  -  MO  Or  FEEDBACK  STATES 
C  NO  -  NO  Or  OUST  INPUTS 

C  NCS  -  NO  OF  COMMAND  INPUTS  *  NO  OF  COMMAND  STATE5 
C  NGLG  -  NO  OF  GUST  LIFT  GROWTH  STATES 

C  NSCRP  -  START  OF  CONTROL  PATE  RESPONSE  IN  THE  RESPONSE  VECTOR 
C  READ  NX.NR.NU.NN.MF.NG.NCS.NGLG.NSCPP 
MIT  I  (M  0  0  013 
C  GO  TO  ?#•  If  INPM.GT.I 


Figure  84.  KONPACT-2  Input  Data  (Employing  DIAK  to 
Compute  Optimal  State  Feedback  Gains) 


C  RE AO  (NORnllWI'l'NX) 

1  n«  S  >  It  91011121  JU1SI61TH1920212223262S26>T2S 
C  200  CONTINUE 

C  f  IS  5TATC  TRANSITION  MATRIX 
READ  TAPE  FOR  MATRIX  F 
C  Gl  IS  CONTROL  INPUT  MATRIX 
READ  TAPE  FOR  MATRIX  fil 
C  G2  IS  DISTURBANCE  INPUT  MATRIX 
READ  TAPE  FOP  MATRIX  G2 
C  XI  IS  INITIAL  CONDITION  MATRIX 
READ  CARD  FOR  MATRIX  XT 

C  XLOXL  IS  STATE  LIMIT  -  RATE  LIMIT  MATRIX 
READ  CARD  FOR  MATRIX  XLO«L 

C  CL  IS  COMmANO  LEVEL  MATRIX 
READ  CARO  FOR  MATRIX  CL 

C  M  IS  STATr  RESPONSE  MATRIX 
READ  TAPE  FOB  MATRIX  H 
C  D  IS  CONTROL  RESPONSE  MATRIX 
READ  tape  FOR  MATRIX  D 
C  AM  IS  MEASUREMENT  MATRIX 
READ  TAPE  FoR  MATRIX  AM 
C  RK  IS  INITIAL  FEEDBACK  GAIN  MATRIX 
READ  CARO  FOR  MATRIX  NK 

III  .266I0<£  OR  126  .lAAlOOE  31  12*  .13S6Q0F  01 

2  I  .SORAflof  00  217  .OflOftTsE  00  22I-.962000F-OI  222  .1I9000E  01 

C  0  IS  QUADRATIC  HEIGHTS  MATRIX 
READ  CARO  FiR  MATRIX  0 

1  1  .MOOOO  F  01  2  2  .IdVlOLE'lv  3  7  .  1  njOOOE-OA  T  T  .SOOOOOE  061010  .600000C  06 
1717  .7S009-E-I31B1B  .100500E-lo202o  . 100000E-1 32121  . 1 OOOOOE-I 12J23  .200000E-13 
2626  .20000  F-1 12626  .6UPOOCE-1 12*27  .10000VF-I02929  .2000ME-123030  .2M0MC*I6 
3333  .ICOOOjE  613636  .lOMIOnE  01 

C  IOUM.0  ANOTHER  RUN 
C  I0UM.1  NO  MORE  RUNS 
C  REAO  IDUM 
1 

END  OF  01  AM  DATA 
SSTOR 


Figure  84.  KONPACT-2  Input  Data  (Employing  DIAKto  Compute 
Optimal  State  Feedback  Gains)  (Concluded) 
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C  PFAO  FOP  ^AT  PROGRAM  i  IiiAK.Ef'iC.I.SA  I  I*f.  UuIA  IS 
♦  FFOC  DATA 

c  PfAO  IF  DMA  IS  ON  CAP'S  ONLY  II)  "i  CMOS  AMO  TAPE 
n/TA  ON  CAP'S  AND  TAPS 

C  IF  DATA  1  .  ON  CAPnc  A’Jij  TA°F  ’O.AU  T  -if.  LAA'L  TO  OPT  MM  D«TA  ON  TAPE 
SYSTEM  NO  6  OVERALL  SYST1'1'  l(  F^A-P  .  RmJCin  CONTROILER)  DESIGN  MODEL 
C  1  PAX  -  MA  <  1  MUM  NO  OF  LVPUMiV  SOLOMON  ITERATIONS 

C  NITM  -  MA«MU»  Of)  OF  COST  CALCOLAMONs 

C  N0p9=0  US'  PROJECTS!)  GpaDTF'NT 
C  NOPRM  00  :0T  USf  PPr)jr/-Tn>  r,P'.ij|F  1 
C  NOCOVtl  N"  COVARIANCE  ANALYSIS 
C  NOCOVS?  C’VAPIANCF  ANAlYSIS 
C  NDC0V*3  SMP  CORRELATION  AVAL  Y*  IS 

C  NREGIN.GT.G  TEST  fop  L  «RGr  o-riiM  ON  F  IV  S  T  IN'-RfME>iT  OF  LAMDA 

C  NBfO!N  =  '1  "0  TFST 

C.  PEAO  IMAA.NITM.NOPp.NOoov.nvFoIv 

?.'  S  0  3  1 

C  NX  -  NO  Of  STATFS 

C  NR  -  NO  Of  RESPONSES 

C  MJ  -  NO  Of  CONTROL  INPUTS 

C  NN  -  NO  Of  DISTURBANCE  ) MOOTS 

C  NFF  -  NO  frro  FORWARD  STATFS 

C  NF  -  NO  Of  FIXED  GAINS 

C  REAO  NX.NJ,NU.NN,NFF.Nr 

2B37  2  2  «  t 

C  PFAO  (N0Ro<  1 1  .  1*1  ,*!XI 

I  2  3  A  S  *.  7  8  9101 1  |?i  AlA|S|F,|  >|  •l9?r,Rl?P?S2A2S?fr>72A 
C  E°SI  -  INITIAL  STFP  Sf/F- 
C  REAO  EPSI 
.S000E  00 

C  AJSTAa  -  LOWEST  COST  F«PECTro 
C  REAO  AJSTaR 
.BSflAE  0? 

C  DROC  -  DESIRED  RATIO  Or  COSTS 
C  READ  DROC 

•  1100E  01 

C  AlAM  -  INTEGRATION  PARAMETER  -  LA'DA 
C  BELT  -  INTEGRATION  STEP  SIZE 

C  ALAND  -  Lower  ROUND  on  LAMDA  fop  the  PRESENT  RUN 
C  READ  ALAM.OELT. ALAMO 

•  1 000E  01  ,?Q0-)E  00  .Od  ooE  00 
C  IF  -  FIXE''  GAIN  ROW  INDEX 

C  JE  -  EIXFD  GAIN  COLUMN  INDEX 
C  READ  ( !E< M iJFIII . I»l .ME) 

111  126  12-  2  1  217  2?1  PP? 

READ  TAPE  FOP  MATRIX  F 
REAO  TAPE  EOR  MATRIX  Gl 
REAO  TAPE  E''R  MATRIX  G2 
REAO  TAPE  FOR  MATRIX  H 
REAO  TAPE  FOR  MATRIX  D 
READ  TAPE  FOR  MATRIX  AM 
READ  CARO  FOR  MATRIX  0 

I  1  .flOOOO.  E  01  2  2  .10'M00E-li)  3  3.  .  10000UE-OB  7  7  .S00001E  041010  .600000E  06 
1717  ,75000  JEM  31 R 1  ft  .  1 0  "(iOt)E- 102  T2F  ,  lOUOOOE- 1 32  12  l  .  100000E-1 12323  .200000E-13 
2424  .200003E-I126P6  .SOlCCCE-l 32727  . 100000E-I 02929  . 200000D- 1 ?3030  .200000E-I0 
3333  ,10000'E  013434  .lOoOOOE  oi 

READ  TAPE  FOP  MATRIX  A*G(OPT IMAL  RICCATI  GAINS) 

GAINS  MATRIk  FOR  CASE  I 
READ  CARD  FOP  MATRIX  AKlKlUI) 

R£A0  CARO  FOR  MATRIX  HK(k?) 

C  IE  <ALAM.GT..99)  GO  TO  10o 
C  100  CONTINUE 
END  OF  FrOC  OATA 
SSTOP 


Figure  85, 


KONPACT-2  Input  Data  (Employing  FFOC  to 
Compute  Reduced  Feedback  Gains) 
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C  DESIGN  l)ST MO  DIAK  FOR  THE  pFMONSl  -  AT  ION  EXAMPLE 
C  REAP  FOR  ''HAT  PROGRAM  (  PI  As  .rr  oC.LSA  )  ME  DAlA  IS 
60IAK  PATA 

C  READ  IF  PITA  IS  ON  r  API'S  ONLY  op  VJ  C/.VOS  AMu  IAR£ 

OATA  ON  PAROS  AND  T A°F 

C  IF  OATA  I f  ON  CARPS  AND  tape  READ  The  LAHEL  10  OHTaIN  PATA  ON  TAPE 
SYSTEM  NO  6  OVERALL  SYSTrM  ((  F?-jk  *  ftrOJCEO  CONTROi.LE  R )  DESIGN  MODEL 
C  READ  DATE  ANO  USE0  IP 
JAN  10.  7S  J  K  MARfsh 
C  NOR  -  NO  ‘  'F  VARIABLES  Rf  I  NO  RLOHO 
C  REAP  NOP 

ia 

C  GO  TO  100  IF  NOP.ro. 1 

C  REAP  (PLPl  1 )  » ITITL  <  I  >  .UNIT  1 1 )  .Y*l  ;( II  .VMAU  It  .SCAL  U  1  •  1=1 <NOP> 

2  01 
3  T 1 

7  DELAOOT 
10  DFLEIOOT 
13  OELA 
16  PELE1 
35  ALPHA 
39  Q 
66  ET  A 1 
47  ETA3 
50  ETA6 
57  OELEO 

C  READ  T.PT.ST.T! .T» 

4.0  0.01  A . '  4  .1,0  -i.O 

C  100  CONTI  HIE 
C  READ  IMAX. ITER. ITERP 
0  0  0 

C  NOCOV  =  1  N'l  COVARIANCE  AMAI  YSIs 
C  N3C0V«?  COVARIANCE  ANALYSIS 
C  N0C0Vo3  S< IP  CORRELATION  ANALYSIS 
C  NSTEPfO  N  STEP  INPUTS 
C  NSTER=1  STEP  COMMANPS 
C  NSTEP*?  STEP  GUSTS 
C  NSTEP=3  DOTH  (|  ANN  PI 

C  NSTEP=4  No  STEP  INRlJTS  -  TRAI.SIrNTS  OnlT 

C  NRANO*0  No  ranpom  INPUTS 
C  NRANOfI  GmSTS 

C  NPRIN*0  O'O  NOT  PRINT  RESPONSES 
C  N°R IN= 1  PRINT  RESPONSES 
C  NDLOT *0  NO  PLOTS 
C  NPLOT=I  C f.LCOMP  PLOTS 
C  NpL0T s?  LINE  PHINTEP  RLPTs 
C  NPL0T*3  HoTh  (1  AND  21 
C  READ  NOCO  /,NSTFR,NRAND.NPRIN.N“i  l!T 
110  1? 

C  I NPKa 1  NE  INPUT  GAINS 
C  INPSiP  NE  1  STARTIN'--  poiiTjnf  (la  i  os 
C  INP<*3  US'  GAINS  I”  S r, -.5A.il 
C  INP<=4  USr  INPUT  GAINS  IN  STOhAGE 
C  READ  INPK 
1 

C  NCONT *0  OPNOT  COMPUTE  n°T|MAL  GAl'iS  -  USE  INP  JT  Gains  AND  DATA  IN 
C  COVARIANCE  «ND  tun.  RF»»P  :Sf  analysis  ONLY 

C  NCONT o |  C'MPUTE  ORTTMA|  Gains 

C  NCONT*?  Compute  OPTIMA!,  GAINS  i I T h  AUTOMATIC  0  EELFCTION  ON  CONTROL  HATES 
C  READ  NCONT 
6 

C  RF AO  FLIG.<T  CONDITION  MIMSES 


Figure  86. 


KONPACT-2  Input  Data  (Employing  DIAK  to 
Evaluate  Time  Responses) 
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r4 


I 

! 

i 


1 

I 


4123‘Jl 
C  MX  - 

MO 

Off 

STATES 

c 

M3  - 

MO 

0»* 

WtSPONSE S 

r- 

MU  - 

MO 

Of 

CONTROL  INPUTS 

c 

MM  - 

MO 

OP 

OISTtlttfUNCf  INPUTS 

c 

NF  - 

MO 

Or 

fef.dsack  statfs 

c 

NG  • 

MO 

Op 

GUST  INPUTS 

C  NCS  -  NO  >>F  COMMAND  INPUT*  t  "0  O'  COmmanJ  STATES 
C  NC-LG  -  No  OF  OUST  LIFT  C.ROtfTM  STATES 

C  NSCRR  -  START  OF  CONTROL  RATE  RERUNS*  IN  ThL  RFSRGNsE  VECTOR 
C  RfAO  NX. N>. Nil. NN«NF,N'i. NCR. NGLG. NSCRR 
2837  2  2P7  I  J  7 
C  00  TO  700  -f  INPK.OT.I 
C  READ  (NON  ii:v.I»l<NX> 

I  2  .1  A  S  -  7  rt  ‘9I0UlP7?HlSH>in-lR2:?2232<,2?2‘'2?P821 
C  20u  CONTI  "JE 

C  F  IS  ST  AT 1  TRANSITION  •■ATolX 
READ  TARE  F  'R  MATRIX  F 
C  G]  IS  CONTROL  INPUT  NATCH 
READ  TAPE  F-9  MATRIX  M 
C  G2  IS  OlSIURflANCL  INPUT  MATRIX 
REAO  TAPE  FtP  MATRIX  F-P 
C  XI  IS  INITIAL  CONDITIO'.  “'THa 
REAO  CARO  F''0  MATRIX  XI 

C  XLOXL  IS  STATE  LIMIT  -  LI'CI  MaIR|* 

REAO  CARO  F'  P  MATRIX  XLIJ»L 

C  CL  IS  COM -AND  LFVfl  MaTRIx 
REAO  CARD  F'R  MATRIX  CL 

1  I-.V8RS7  K  a-. 

C  H  IS  STAT •'  RESPONSE  matrI- 
READ  TAPE  F  ip  MATRIX  H 
C  0  IS  CONT-OL  HFSPQNSF  -ATCJX 
READ  TAPE  F'R  MATRIX  I) 

C  AN  IS  MEASUREMENT  matrix 

READ  Tape  F'R  matrix  am 

C  H<  IS  INITIAL  FEEDPAC.R  GAIN  MATH ' 

READ  CAPO  F  >U  MATRIX  rk 

111-, 78129  F  Cl  IPX,  .I|oAti'£-  OP  t?"  .  lRhRiloE  M\ 

2  1  ,01 780  E  SI  217  PRASE  «?  >2l-,(>4(i|6or-j|  22?  .UAROMuE  Cl 

C  0  IS  DUM)  ATIC  XEK-MTS  matrix 
REAO  CARO  F'R  MATRIX  N 

C  1DUM*0  ANOTHER  run 
C  10UM*1  NO  MURE  RUNS 
C  REAO  IOUM 
I 

FND  OF  01  AM  OATA 
SSTOP 


i 


> 


Figure  86.  KONPACT-2  Input  Data  (Employing  DIAK  to 
Evaluate  Time  Responses)  (Concluded) 
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C  DESIGN  USING  DIA«  FOP  THE  nfMOHsU'AT  ION  EaamplE 
C  READ  FOR  'HAT  PROGRAM  (  0  I  A* ,  rnc.LSA  I  THF  um  Ig 
SOIAK  DATA 

C  READ  IF  OaTA  IS  ON  CARi-.S  ONLY  0-  ,iN  C/.-N> S  Am,j  TAPE 
OATA  ON  CAP  IS  AND  TAPE 

C  IF  DATA  Is  ON  CARP'S  AND  Tape  ofll  T  ME  LAm'.l  IP  OBTAIN  DAI  A  ON  1  APE 
SYSTEM  NO  6  OVERALL  SYSTEM  ((  r-ji.jH  .  REDUCED  CONTROLLER!  DESIGN  MODEL 


C  READ  DATE  AND  USEP  ID 
JAN  10.  7S  J  A  MAHFSH 
C  NOP  -  NO  OF  VARIABLES  »FIl|F.  PLOTT.U 
C  READ  NOP 


GO  TO  100  IF  NOP.Fp.o 
READ  (PLPl  I  >  .ITITL  1 1)  •'PUT  ( I )  .  Y 
REAO  T.DT.ST.Tl .T? 

100  CONTINUE 

REAO  IMA*.  ITER.  ITFP() 

?ulS  0 

C  NOCOV* 1  N"  COVARIANCE  ANALYSIS 
C  NOCOV*?  COVARIANCE  ANALYSIS 
C  NOCOV*3  S'lP  CORNEL  AT I^N  ANALYSIS 
C  NST£P=0  NO  STEP  INPUTS 
C  NST£P*1  STEP  COMMANDS 
r  NSTEP*?  STFP  GUSTS 
NSTEP*3  BoTH  ( I  ANn  ?! 

No  STEP  INPUTS  -  TRANSIENTS  ONLT 


IM  I  I  .YMAAH  I  .SCAL  <  I  1  .  I  si  .NOP  I 


NSTEP** 

NRAN0*0 

NRANO-1 

N»RIN=0 


Nl  RANDOM  INPUTS 
GNSTS 

. . .  DO  NOT  PPIMT  Hf SPOUSES 

C  NRPIN*I  PRINT  RESPONSES 
C  N°L0T  =  0  NO  PLOTS 
C  NDLOT*l  CALCOMO  PLOTS 
C  N=10T*?  LINE  PRINTER  P|  OTi 
C  NDL0T*T  BOTH  (I  AND  ? I 
C  READ  NOCOu,NSTEP.NRAN!),NPPt\.MPLOT 
3  0  0  0  0 

C  INPK*1  N£ ;i  INPUT  GAINS 
C  INPK*2  NE.'  STARTING  ROUTINE  GAI'>S 
C  INPK.3  USo  GAINS  IN  STORAGE 
C  1NPS»*  USE  INPUT  GAINS  IN  STORAGE 
C  READ  INPK 
I 


C  NCONT *0  DONOT  COMDMTF  OPTIMAL  r,cl’:S  -  uSr  ImP'JT  GAINS  AND  DATA  IN 
C  COVARIANCE  AND  TIME  RESPONSE  ANALYSIS  ONLY 

C  NCONT*  1  COMPUTE  OPTIMAL  GAINS 

C  NCONT*?  COMPUTE  OPTIMAL  GAINS  «IT«  AUTOMATIC  0  SELFCY 1  ON  ON  CONTROL  RAT 
C  REAO  NCONT 
1 

C  REAO  FLIGHT  CONDITION  MUMMER 

*12301 

C  N*  -  NO  OF  STATES 
C  NR  -  NO  OF  RESPONSES 
C  NU  -  NO  Op  CONTROL  INPUTS 
C  NN  -  NO  OF  DISTURBANCE  INPUTS 
C  Nr  -  NO  OF  FEEDBACK  STATES 
C  NG  *  NO  Or  GUST  INPUTS 

C  NCS  -  NO  OF  COMMAND  INPUTS  *  NO  OF  COMMAND  STATES 
C  NGLG  -  NO  OF  GUST  LIFT  GROWTH  STATES 

C  NSCRR  -  START  OF  CONTROL  PATE  RESPONSE  IN  THE  RESPONSE  VECTOR 
C  REAO  NX.NPfNU.NN.NF.NG.NCStNGLG. NSCRR 
2B37  2  ?28  !i  t  <M3 
C  GO  TO  201 >  TF  tNPK.r.r.l 


Figure  87.  KONPACT-2  Input  Data  (Employing  DIAK  to 
Evaluate  Covariance  Responses) 
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C  BEAD  <NORnm«I*l.NX» 

1  !HU?I  VIOIM?!  JU1S161TIOI  WI??;>32(.?S76?72« 

£  Zlt  CONTINUE 

c  r  is  state  transition  matrix 
READ  TARE  FOR  MATRIX  F 
C  61  IS  CONTROL  INRUT  MATRIX 
READ  TAPE  E OR  MATRIX  «1 
C  6?  IS  DISTURBANCE  INPUT  MATRIX 
.READ  TAPE  FOR  MATRIX.  0? 

C  XI  IS  INITIAL  CONDITION  MATRIX 
READ  CARO  EnR  MATRIX  XI 

C  XLOXL  IS  STATE  LIMIT  -  PATE  LIMIT  MATRIX 
READ  CARO  EUR  MATRIX  XLOXL 

C  CL  IS  COMMAND  LEVEL  MATRIX 
READ  CARO  EAR  MATRIX  CL 

C  H  IS  STATE  RESPONSE  MATRIX 
REAO  TAPE  ro»  MATRIX  M 
C-O  iS  CONTROL  -RESPONSE  MATRIX 
KAO  TARE  EOR  MATRIX  0 
C  AM  IS  MEASUREMENT  MATRIX 
MgM  TARE  FOR  MATRIX  AN 
C  MR  IS  INITIAL  FEEDBACK  RAIN  MAToIx 
REAO  CARO  FOR  MATRIX  RK 

C  a  IS  QUADRATIC  ME  I  AMTS  MATRIX 
REAO  CARO  FOR  MATRIX  0 

..11  •QIMQaC  II  1  ?  . Lif atOE-IO  3  3  .100088E-8S  T  7  .S08888E  OMOlO  .600000E  06 
1T1T  .TJMIOE-lllIH  .lOPaOOE-lOTfEC  .  1 00008E-I VI 21  .IC88C0E-I 1*323  .20090VE-13 
2*2*  .2888811 1- 112*2*  .ROf  6«CE- ! 32T27  .IR0008E-I  02A2R  .200880E-1238J0  .280080E-10 
_  031i  ORIAA8E  8U*3*  ^IDCOODE  VI  . 

C  IDUM*t  ANOTHER  RUN 
-£  I  DUN*  L  NO  HOWL  RUNS 
C  REAO  IOUM 
I 

MTOR 


i 


Figure  87.  KONPACT-2  Input  Data  (Employing  DIAK  to 
Evaluate  Covariance  Responses  (Concluded) 
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iPfl? 


NWXA^VRiUMStWA?'.'''''  i  Vi 


C  ISA  INPUT  DATA 
♦  LSA  data 

C  READ  If  DATA  is  ON  CAROS  ONLY  OR  ON  CARDS  AMO  TAPE 
OAT  A  ON  CARDS  AND  TAPE 

C  If  DATA  IS  ON  CAROS  A  Up  TAPE  READ  TnE  LABEL  TO  OBTAIN  DATA  ON  TAPE 
SYSTEM  NO  4  PLANT-E420  <  CONVERTED  PLEXST  A8  CSA  A/C  .  ACTUATOR  ♦  GUST  MODEL  I 
C  INSERT  DATA  ON  CAPOS  ir  ANY 
f  NO 


C  NANE 
NAME 

V 

LIST  DATA 

0 

fOR  LSA 

UCIQOT 

UC3D0T 

UEA30I 

uesdot 

UESOOT 

Ut TOOT 

UEAOOT 

UEROOT 

UE1DOOT 

UEI IOOT 

UEIZDOT 

UE1300T 

uciadot 

UEI SOOT 

UEI SOOT 

UCI 

lie  3 

U£4 

IlCS 

UE6 

UET 

uea 

HER 

UCIO 

util 

mu 

UEI  3 

U£|  A 

UE15 

UEI  6 

XA 

*E! 

XfO 

PI 

P> 

P3 

Pa 

PS 

P» 

VO 

•1 

T  l 

B2 

T? 

63 

T3 

BA 

14 

as 

TS 

B1D3T 

TIOOT 

B2DOT 

T2D0T 

B303T 

TJOOT 

OAOOT 

TAOOT 

BSDOr 

T500T 

ALPHA 

OCG 

DEL  A 

Dilei 

OELEO 

■R 

v  DOT 

A2lf 

A2IR 

AfuS 

Trus 

END  Of 
♦  STOP 

UpfLA 

LSA  DATA 

UOFLf I 

UOCLEC 

etas 

Figure  88.  KONPACT-2  Input  Data  (to  Prepare  Frequency 
Domain  Data  for  LSA  Program) 
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OATA  fOR  LSA  (sroaTrM=?fiff  n 
1  Of  NT  |F  1C  AT  10KMFCSA 
RfAO  MSI 

POINT  STSTfM  MATPT< 

COMPUTE  POLrS 
compute  /Tt»os(oi/F»ACi 
compute  pso 

PLOT  osniUNEAP-n.l.U-.rt.'si.i  .t.-JE  I C> » 
COMPUTf  ?EOnS(Tl/‘T»f<> 

COMPUTE  pso 

PLOT  >SO  (LINE  »n  •0.1.1'  ,9.t.  .v.jE  It,  I 
COMPUTE  ZEPOScn/'ll'ELAI 
COMPlJ  i  F  PSO 

PLOT  oSO<LINFAO.;.l.!".f.  17| 

COMPUTF  ZEPOSlTI/.Jf'ELU 
COMPUTE  PSO 

PLOT  osoiLTNF A°.3.»  .)  .  i.oE  111 

COMPU'E  ZEBOS(«I/U|>fLF  II 
COMPU’E  PSO 

PLOT  PSOILINCAO.U.I.I '.Cit..  i.-.OE  ITj 
COMPUTE  ZEP0S(Tl/UOFLf I> 

COMPU’E  PSO 

°LOT  PSO  (LINE A3. il.  1  .] ir •  1 .  )E  lsl 
COMPUTE  ZEROS (Q/ETf 01 
COMPUTE  PSI1 

PLOT  'SOIL  tl!EA“ .0. 1  .  1  ii.P.i' .  ..7.(11 
COMPUTF  ZFPOS (O/UOFL A } 
compute  psn 

PLOT  -SO(LINFAP,.1.t.lL.C.;.  ,.|..)E  tiS) 
COMPUTE  ZEROS  (O/UOfl.Ett 
COMPUTE  PSO 

PLOT  PSIHLINEAP.O.t.lTT.iV.fl.'MA.OE  OS) 
COMPU'E  ZE»OS  (ilf  3/fTAF.i 
COMPU'E  PSO 

PLOT  ssihliNEAO.a.,  .n.o.o.  :.-»,0E-01* 
COMPUTE  ZEPOSUIEO/UOFlai 
COMPU’E  PSO 

PLOT  OSOILINEAB. 3.1. 1(1.9..*.,  .I.CE  •}<.* 
COMPU'E  ZEROS (UE3/CCELE 1 1 
COMPUTE  PSO 

DLOT  "SO(L1NEAP.0.1 .  |A. ,)•(...  .a.cE  0»l 
ENO 


Figure  89.  LSA  Input  Data  (to  Evaluate  Power  Spectral  Density) 
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{Cruise  Flight  Condition) 
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Figure  90.  C-5A  Vehicle  Name  List  Table  and  Quadruple  Data 
(Cruise  Flight  Condition)  (Continued) 
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Figure  90.  C-5A  Vehicle  Name  List  Table  and  Quadruple  Data 
(Cruise  Flight  Condition)  (Continued) 
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Figure  90.  C-5A  Vehicle  Name  List.  Table  and  Quadruple  Data 
(Cruise  Flight  Condition)  (Continued) 
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Figure  90.  C-5A  Vehicle  Name  List  Table  and  Quadruple  Data 
(Cruise  Flight  Condition)  (Continued) 
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Figure  90.  C-5A  Vehicle  Name  List  Table  and  Quadruple  Data 
(Cruise  Flight  Condition)  (Continued) 
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Figure  90.  C-5A  Vehicle  Name  List  Table  and  Quadruple  Data 
(Cruise  Flight  Condition)  (Concluded) 
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Figure  91.  Overall  System  Design  Model  Name  List  Table  and  Quadruple  Data 
(See  Figure  83  for  KONPACT-1  Input)  (Continued) 
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Figure  91.  Overall  System  Design  Model  Name  List  Table  and  Quadruple  Data 
(See  Figure  83  for  KONPACT-1  Input)  (Continued) 
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Figure  91.  Overall  System  Design  Model  Name  List  Table  and  Quadruple  Data 
(See  Figure  83  for  KONPACT-1  Input)  (Continued) 


Figure  91.  Overall  System  Design  Model  Name  List  Table  and  Quadruple  Data 
(See  Figure  83  for  KONPACT-1  Input)  (Concluded) 
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SECTION  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  volume  provides  two  demonstration  examples  to  aid  the  user  of 
KONPACT  in  preparing  input  decks  and  executing  the  programs. 

The  preliminary  emphasis  in  this  study  has  been  interface  software  develop¬ 
ment  and  the  demonstration  of  its  use.  These  have  been  achieved.  For 
future  work  it  is  recommended  that  the  ALDCS  design  using  the  FLEXSTAB/ 
LSA  data  be  reinvestigated  and  refined  to  meet  the  torsional  moment  con¬ 
straint  and  phase /gain  margin  requirements. 
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